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PREFACE 


This little book, like so many of its kind, had its roots in a course of 
lectures. I gave the lectures, by request, to an audience consisting mainly 
of graduates in physics, soil mechanics and the various branches of civil, 
mechanical and electrical engineering. The listeners had begun to find 
that progress in their researches and further studies was being impeded 
by their inability to handle Bessel functions; it was my task to remedy 
this defect in their equipment as best I could. 


Somewhat to my surprise the audience turned up in considerable 
strength; to my even greater surprise a high percentage stayed the course 
without hope of ulterior reward beyond the acquirement of wisdom. I 
was flattered into believing with venial weakness that I might have done 
worse; and maybe I lent all too readily “ a credent ear ” to the suggestion 
that the lectures be put into book form. However, now that the step has 
been taken, I am under a moral obligation to answer two questions. 

The book is definitely mathematical and is addressed to technicians. 
. ‘“mediately poses the question, what should be the mathematical 
equipment of the man who hopes to benefit by reading it? The lectures 
were designed to fall within the circle of ideas of those who had taken the 
mathematical course that usually-gogs with-4 degree course in physics or 
engineering. More specifically, no use was to be made of contour integration 
or of the complex variable in the analytic sense. On the other hand I 
assumed a knowledge of ordinary differential equations with constant 

further’ ^ OCCm "} the fcheor y of beams or of simple circuits. I 
urther presumed my audience to know that such partial differential 

equations as normally arise can be solved as the product of functions of 

the independent variables and might involve Fourier series 

*r-.M ked °? the IeCtUreS W ' th a g00d deal ^trepidation, for reasons 

mativeti, i„ g IZTitT^ IT °' “ ' UnCt r U *“ kast ^ 

bc <:ven in A° p° ! 10 
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much cogitation I decided that the special circumstances of the task 
justified the heterodoxy of making the recurrence formulae the starting 
point. I had no reason to regret my decision, for as pointed out in the 
bibliographical note, I found I had the moral support of the argute E. B. 
Wilson. 

I inserted the initial chapter, treating the ancillary functions, on 
pedagogical grounds. It is better to do a little preliminary spadework to 
make sure that the soil is ready for sowing, rather than to interrupt one’s 
discourse by parenthetic paragraphs and distracting footnotes. There is 
the additional defence that the Gamma function is in itself sufficiently 
interesting to justify the expenditure of a few hours on mastering its salient 
properties. 

Coming to the main obstacle, how to demonstrate the existence and 
nature of the zeros of the Bessel functions, I naturally inquired how other 
writers had approached it. In three current books I found that one em¬ 
ployed Bessel’s original method that a beginner would hardly relish; 
the second advocated the uninspired method of plotting from tabulated 
values; the third went on the principle ignotum per ignotius by borrowing 
the unproven asymptotic values. That removed all qualms I might have 
had about using the oscillation theory. The majority of English mathe¬ 
maticians are indebted to the standard treatise by G. N. Watson rather 
than to the continental school, and I think that if its author had felt more 
drawn to the Sturmian theory, other writers would have adopted it. In 
any case, mathematical physics has in recent years given the theory a 
renascence, and my audience was pleased to learn that a differential 
equation can be made to impart information by methods other than solving 
it. The slight advance on what is usually taught concerning differential 
equations is given in the second chapter, and the unbiased can hardly 
fail to recognize in it an instrument admirably suited to its purpose. 

The year of publication marks a centenary; F. W. Bessel died at 
Koenigsberg on 17th March, 1846. 

Concerning tables, it is confidently to be expected that the British 
Association Tables Committee will shortly issue a second volume on Bessel 
functions in continuation of their previous work. The new volume will 
tabulate the four functions J, Y, I and K, the argument running from 0 
to 25 and the integral order from 2 to 20. The entries will each contain 8 
reliable figures and there will be provision for interpolation. 

F. E. RELTON. 

Watford, 

July, 194G. 
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CHATTER I 


The Error Function ; 
Beta and Gamma Functions 


11. The study of functions. 

The functions that one encounters early in one’s studies are usually 
termed “ the elementary functions ”. They include the trigonometrical 
functions, the hyperbolic functions, the exponential and logarithmic 
function. Functions other than these are termed “ transcendental 
functions ” and for the most part they fall into one or other of two 
categories; either they are defined by an integral or they satisfy a 
differential equation. It occasionally happens that a function falls 
into both categories. 

A mathematician, faced with the study of a new function, lias 
very much the same attitude that a physicist or a chemist would adopt 
towards a substance presented for examination. The latter would 
inquire whether the substance was crystalline or necessarily amor¬ 
phous, whether it was soluble and what were its solvents; what was 
its density and thermal conductivity, and so on. A mathematician 
would inquire whether there were any values of the variable for which 
the function vanished, or alternatively became infinite; in his own 
terminology, he would seek for zeros and poles. He would be interested 
to know whether the function had an addition formula, by which we 
mean an expression for/(a + b) in terms off (a) and /(&), the simpler 
the better and preferably algebraic. Thus the trigonometrical function 
smx has the algebraic addition formula 

sin (a: -f- y) = sinx(l — sin 2 y)» + s in^(l — s in 2 x)*, 
which can also be written in the differential form 


sm(x + y) = sinx — siny -J- siny ~ sinx. 

a V ax 

2S Sr'lT '• ° u!d be erested t0 whether the function 
to&TlT r m ‘ a ’ and m P artio,,lar whether it was asymp-, 
‘° SOme Slm P le ex P ress ,on for large values of the variable just 
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as coshx is asymptotic to \e x . Finally he may decide to tabulate the 
function, though this is usually the last thing he does and many func¬ 
tions have been studied which nobody has yet troubled to tabulate. 

The investigation of these and other points of interest is rarely 
carried out by evaluating the aforementioned integral or by solving 
the differential equation which defines the function. There are other 
methods of approach and we must be prepared to encounter these in 
the following pages. We shall be concerned almost solely with the 
so-called Bessel functions, of which there are several kinds; and it 
may as well be stated at the outset that the graph of a Bessel function 
usually looks like the distorted sine-curve that represents a damped 
oscillation. The amplitude steadily decreases and there is this differ¬ 
ence, that the “ period ” is not quite constant, so that the function 
cannot legitimately be termed periodic. But before we come to the 
Bessel functions there is a certain amount of spadework to be done. 


1-2. A useful limit. 

We begin with a limit for which we shall find almost immediate 
use. If n is a positive number we have 

(1) x n e~ x —> 0, x -> 0, n > 0. 


The arrow notation signifies “ tends to ”, “ approaches ”, or “ has 
the limit ”. The above line signifies that when n is positive and x tends 
to zero, then the expression at the left also tends to zero. The result 
is obviously correct since e~ x tends to unity and x n tends to zero. 
Similarly we have the not quite so obvious result 

(2) x n e~ x -*■ 0, x->cc, n > 0. 


Here x n tends to become indefinitely large and e~ x tends to zero. We 
establish the truth of the statement by writing 



x n x n 

x n e~ x = — =-. 

e 1 1 + x -f- x 2 /2 + . . . 


On dividing numerator and denominator by x n the numerator becomes 
unity. The denominator becomes a limited succession of negative 
powers followed by an infinite succession of positive powers of x. The 
negative powers tend to zero and the positive powers tend to become 
infinite. Thus the fraction tends to zero and the truth of (2) is estab¬ 
lished. 
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EXERCISES 

1. Sketch the graph of x n e~ x when x is positive and n is (i) positive, (ii) zero, 
(iii) negative. Prove that in (i) there is a single maximum; but in (ii) and (iii) the 
value steadily decreases as x increases, whilst in (iii) the initial value ia infinite. 

log# logx 

2. Prove that--> *>, x -> 0. Also -> 0, j: -> ce. 

X * 

3. Deduce that as the number n tends to become very large, the nth root of 
n tends to unity. 

log l log 2 log 3 . , 

4. Prove that the series — -f —^—|-^—j- ... is divergent and deduce 

that the product 1* . 2& . 3*. 4* . . . is infinite. 

5. Prove that as a -> oo, -*■ 0. 

6. Sketch roughly the graph of the function (logx)/x for positive values of 
x. Prove that it has a single maximum value e~ l . 

1-3. The error function, Erf (a:). 

The error function is one of the class of functions defined by an 
integral. It plays a prominent part in the theory of probability, refrac¬ 
tion of light, conduction of heat, and so on. Its definition is 

(1) y = Er((x) = f e-t'dt. 

Jo 

The graph of the integrand is shown in fig. 1. It is symmetrical and 
has a maximum value of unity when t is zero; it is doubly asymptotic 
to the horizontal axis and has two inflections. The variable t is a dummy 
and could equally well be called 8 or any other symbol. The integral 
represents the area included between ordinates at the origin and at 
some variable distance x. Its value will be found tabulated in almost 



x 
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any work on probability and in many sets of tables. The integral can 
be evaluated by expanding the exponential in an infinite series of 
ascending powers of t and integrating term by term; but the result is 
not very informative. 

It not infrequently happens that where an indefinite integral can¬ 
not be evaluated in a succinct form, a neat result can be obtained 
when the integral is definite. We now propose to show that the whole 
area under the curve in the first quadrant is finite and has the value 

( 2 ) f e~ lt (ll = i\/7r. 

•'O 

This last integral is called an infinite integral on account of the upper 
limit; if it exists it is regarded as the limit of the same integral, taken 
between 0 and a, when a is made indefinitely large. Hence we may 
write 

N a ==^ e~ x *dx =j^ e-*'dy, 

N = f e~ x *dx = f e~*'dy\ 

Jo J 0 

N a -> N as a -> oo. 

Here N is a pure number and N„ is a number dependent on a. We 
may write 

N a 2 = re~ z 'dx re-v'dy = f* [%-*•-»'dxdy, 

Jo Jo Jo Jo 

the double integral being taken over a square of side a. On converting 
to polar co-ordinates r, 8 we have e~ x *~ v * = e -r * and the area-element 
dxdy is replaced by rdrdd. The limits for r are 0 and a. and the limits 
for 8 are 0 and ; but as the polar co-ordinates cover only a quadrant 
instead of a square we write 

N a 2 = re-’*rdrdd + R, 

Jo Jo 

where R is known as the remainder. It represents an addition due to 
the area between the quadrant and the square (fig. 2 ), and our first 
step is to show that the value of R is ultimately negligible, or R -> 0 
as a -> oo. Here is the proof. At any point within this area we have 

x 2 + y 2 > a 2 , e-** *' < e 

and since the area over which this is to be integrated is certainly less 
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than half the square, we have R < Ja 2 e a *. Hence R -> 0 as o —> x> . 
We may now let a become infinite and rewrite a previous line as 

N 2 = / ddj re-^dr = , 

whence N = as stated earlier. 

As an example of its application, consider the following intractable- 
looking integral which occurs in the theory of heat conduction. 



Fijj. 2 

The substitution xy = a is equivalent to 

dx __ dy 

xr a 

whence 

^ = 2 a —_ 21 

da a 

The solution of this simple differential equation is I = cc~- a . The 
arbitrary constant c can be determined by putting a = 0. This gives 
c = W* and I = Zv'ttc- 20 , b 
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EXERCISES 


1. Sketch roughly the graph of y = Erf (x) for positive values of x. It close!}' 
resembles the right half of y = tanhx; at what angle does it leave the origin? 

dl 

2. If / = I e~ x ’cos2bxdx, verify that ^ + 2hi = 0; hence deduce that 
1 = W Tre" 6 *. 

J r * t t 

e~ a,x ' cos 2 bxdx = — e“ 6 ^ a ^ 
o ^ 

4. Integrating No. 2 between 0 and fc, deduce that 


i 


.-*• 


, sin 2br 


L 3 


b> 


dx= y/r. 6 - - + 


■0 X l 3 ‘ 5.2! 

6 . Use the substitution x = t — a to establish the result 


f e '.at-t'(i l= JVne 0 *. 
da ^ 

fi. The error function is sometimes defined as Erfc(x) = / e~ r dl. Deduce 
Erf (x) = i C 7 t — Erfc(x). x 

(The function is sometimes called the “ nonnal error ” or “ probability inte- 

2 f* 

gral ” and is tabulated in the form 7= ^ 7 ^/ e~ l ‘dl. This is asymptotic to 
unity for largo values of x and reaches the value 0-99998 when x = 3.) 


1-4. The Gamma function. 

The Gamma function is one of the class defined by an integral. 
It is known that it does not satisfy any differential equation with 
rational coefficients and it is sometimes called the Eulerian integral of 
the second kind. It is one of the earliest transcendental functions that 
one encounters and it is certainly one of the most interesting. 

If n is a positive number, we have as the definition 

(1) T(n) = J l n ~ 1 e~ , dt, n > 0. 

The graph of the integrand for various values of n has already been 
mentioned. The function represents the area to the right of the origin 
for various values of the parameter n. The function T(n) is necessarily 
positive for all positive values of n, for the integrand is positive through¬ 
out the range of integration. The variable t is obviously a dummy. 

Note that the index of l is m — 1 when the argument of the function 
is n. If n is less than unity the integral is of the type known as im- 
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proper since the integrand then tends to infinity at the origin. If n is 
less than, or equal to, zero the integral fails to converge. 

We achieve a first result by equating n to unity. This gives 

(2) T(I) =f*e-‘dt = 1. 

If we integrate by parts we have 

f n-ie~ l dt = - er* + f- e~‘dt. 

J n J n 


On inserting the limits we can write 



t n ~ 1 e~ l dt 



t n e~*dt. 


The first term on the right has already been proved to vanish at both 
limits. We thus have 

(3) nr(n) = r*(n +1). 

This is known as the difference equation for the Gamma function and 
it is one of its most important properties. We proceed to investigate 
some of its consequences. 

(i) Putting n = l gives F(2) = 1T(1) = 1. An application of 
Rolle’s theorem to the two results T(2) = 1 = T(l) indicates that 
F(n) has a minimum value. Actually it occurs when n = 1-4616 

and the corresponding value of r(n) is 0-8856 . 

(ii) By repeated application when n is an integer, e.g. n = 5 we 
get 

T(5) = 4r(4) = 4.3r(3) = 4.3.2r(2) = 4! 

It is easily seen that in general we have T(n) = (n — 1)! when n is an 
integer. As a matter of history, the problem of finding a function of 
a: that should be continuous when a; is positive and take the value x» 
when x is an integer was solved by Euler in 1729. The Gamma function 
can accordingly be regarded as a sort of generalized factorial. 

(in) When n is not an integer, e.g. n = f, we can ultimately reduce 
the argument below unity. Thus T($) = *r(t) = |. *r<i). It follows 
that for purposes of computation the function need only be tabulated 

reduced^ 1 ^ 6 ° *° ^ Ifc ^ appear ,ater that even this range can be 

(iv) By an uncritical use of the difference equation we can carry 
the function over into negative values of the argument. If, for example 
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wo put n = ?. in the relation T(«) = (n — l)(n — 2 )T(n — 2 ) we get 
F(') = (— n)(— 3 )F( — 5) whence T(— £-) is defined as , n 0 T(^). It should 
be noted that the integral definition no longer holds good since the 
integral fails to converge at the lower limit when n = — f. 

With negative integral values we encounter a series of infinities. 
Thus 

T(2) = 1 . T(l) = 1.0. F(0) = 1.0. -1 . T(-l) 

= 1 . 0.-1 . — 2 . T (— 2 ) 

and so on. It follows that T(— n) is infinite if n is a positive integer 
or zero. 

The integral defining the Gamma function can be given various 
forms by changing the dummy variable. The most useful of these is 
the substitution t = x 2 , dt = 2x dx. We then have 

(4) I» = 2 f x 2n ~ 1 e- x 'dx. 

Jo 

With the particular value n = \ we recover the error function, which 
leads to 

Hi) = 2 f"e-**dx. 

Hence the important result 

(5) T(J) = y/ir. 

It is now possible to form an idea of the graph of the function, 
which is indicated in fig. 3. In virtue of the difference equation the 
value tends to infinity with increasing n, and also when n tends to 
zero. The reader should satisfy himself that each branch on the left 
lies wholly on one side of the axis, and the stationary values gradually 
approach indefinitely close to the axis, though the first two are con¬ 
siderably greater in absolute value than the stationary value on the 
right. Incidentally the stationary values do not occur at corresponding 
points; nor are the branches on the left symmetrical, as they are in the 
graph of seex. It is an interesting exercise to prove that the equation 
T(n) = c has an infinite number of roots for all values of c, positive 
and negative. The Gamma function has no real zeros except at nega¬ 
tive infinity. 

The name and notation of the Gamma function are due to Legendre. 
It is as well to point out, since the reader is likely to consult books 
other than this, that an alternative function is sometimes used, known 

(C 150) 
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as Gauss's FI function. It only needs to be remembered that the dif¬ 
ference is one of notation only and that Il(«) = r(« -f- I), so that for 
positive integral values we have ri(n) = nl = T(« -f- 1). 



1-5. The Beta function. 

The Beta function is one of the class that is defined by a definite 
integral. It is sometimes known as the First Eulerian integral and it 
involves two positive parameters which we may take to be p a The 
definition is 

(1) B(p, q) =jfV-i( 1 - ty-'dt, p,q> 0. 


The reader should sketch the possible types of integrand for different 
values of p and q- but a graph of the Beta function itself is not possible 
since it involves three-dimensional representation. The function 
represents the area between the ordinates at t = 0 and t = 1, and the 

X the 1S im T^ " Gither ? ° r * ° r b0th are Iess than 

ToTj nr h Th beC ° meS impr ° per afc one or hoth limits; and if 
f . q, or both are negative or zero the integral fails to converge The 

Integra 1 can be evaluated in finite terms if q is an integer- expansion 

by the binomial theorem permits termwise integration, 0 but the result 

llhed Jf UC ‘p ? Peat latGr tHat the Beta function is too closely 
« Ganima f “ * h - g -t importance 

(O 150) 


2 
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1-6. Change of variable. 

The substitution t = 1 — x, dt = —dx gives 

(1) B(p,q)=[ a^l — xY~ l dx = B(q, p), 

Jo 


so that the function is symmetrical in its two parameters. 
A more fruitful substitution is 

* = sin 2 0, 1— f = cos 2 0, dl = 2 sm 6 cos Odd. 

This gives 

(2) B(p, q) = 2 f sin 2p_1 9 cos 2fl_1 Odd = B(q, p). 

Jo 


This leads to a number of “ reduction formulae ” that one encounters 
when learning the integral calculus. We shall return to these later. 

A third useful substitution is t = y/(l y). This leads to 


B (p, ?) =/ 
Jo 


rC 


y* 




0 (1 + yY+* 


dy. 


In particular, if 0 < (p, q) < 1, p + q = 1 we have the important 
result 



B( 1 - q, q) 



1 -f- y smqn 


This last integral is a well-known result, easily established by contour 
integration. Unfortunately, elementary proofs are not so simple. They 
can usually be found in advanced works on the calculus; in particular, 
see Edwards, Integral Calculus, Vol. 2, p. 61. 


Example .—Consider the integral 

r\* _ r\ n 

1=1 VcotOdO = / coal 0 Bin” * 0efO. 

Jo Jo 

Here we have 2p — 1 = — 2q — 1 = §» whence p= \,q= J- This gives 


1-7. Connexion between B and T functions. 

The two forms 1-4(4) 

CD 

T(p) = 2 f^x^-'e-z'dx; F(q) = 2jf y^-'e-*'dy 

r(n)r(g) = 4 [* [ e- x '~'>'x l1, - l y' lQ - l dxdy. 

Jo Jo 


lead to 
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On converting to polar co-ordinates x = r cos 6, y = r sin 6 we can write 

r(j>)r(?) = 4 [*’f V'V *’**- 1 cos*’- 1 * sin *-'0drd0. 

J 0 Jo 


The reader is left to justify the ignoration of the remainder R by the 
process used when studying the error function. We now have 

r(p)r( 9 ) = 2 f cos 2p-1 0 sin to-'ddd X 2 f e~ r V’ +2 «- 1 dr. 

J o J 0 


The first integral is known to be B(p : q) and the second integral is a 

form of r(p + q). Hence I»r( ? ) = B(p, q)V(p + q), or as it is usually 
written 



B(p, q) = 


I»r(y) 
F(P + <l) 


The result emphasizes the symmetry of B(p, q) in its two parameters. 

It also explains why the Beta function is of little consequence on its 

own account. If we put p + q = 1 the denominator becomes unity 
and we have, using 1-6(3), 


( 2 ) — q) = B(q, 1 — q) 

In particular, if p = J = q we have 


7r 

sinqn 


mm = 


77 


sin l 


7T 


which is our previous result T(^) = yV. 
It appears from the result 


Sin 7177 


oonnerting tho Gamma functions with the circular functions that for 


so that 


mm = ~, 

sin^TT 



table in this range was 



calculated 


■ny/2 

wv 

to twenty figures by Gauss. 



12 


APPLIED BESSEL FUNCTIONS 


1-8. The duplication formula. 

A duplication formula is one which expresses the function of a 
doubled argument in terms of the function of the single argument. As 
a simple illustration, the function cos2j is expressed in terms of cosz 
by the duplication formula cos 2x = 2 cos 2 j — 1. 

We have by definition 

B(n, n) = ['x n ~'( 1 - x)"~'dx = f\x - x-) n -'dx. 

Jo J o 


The substitution 
gives 



b(«. »> = p-L/j 1 - y"-r- ld, J = ^-Jj 1 - y'-'! 


n —1 


dy, 


since the integrand is symmetrical about the origin, 
st it ut ion 



The further sub- 


givcs 


*(».») = - tr-H-'jt = A, B^, i) 

M •'0 “ 


In terms of Gamma functions this is 

!»!» = Tin) H|) 
r(2n) 2 2 "- i r(n + A)’ 

Hence 

(1) y/it • r(2n) = 2 2 »-T(n)r(n + A). 

As it is no part of our programme to develop the Gamma function 
beyond what is requisite for the comprehension of the Bessel functions 
we refrain from pursuing the matter further; but the student to whom 
the subject is new would be well advised to follow it a little further in 
some work on advanced calculus. 
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EXERCISES 


1. J x n ~ 1 e~ tx dx = -pr* & > 

2 . IW(§) = -L nt)r(i) = wi2»im 3*{F(*)>* = -« 2 ir(j). 


3. f '"cos'* e Bin 6 ' 3 e dQ = 

Jo 27Vw 

4. If n i6 a positive integer, 2 n r(n -f- J) = 1.3.5 . . . (2« — 1) Viz. 

5. Provo that f < P-I (a — ty*~ 1 dt = B(p, q)a IH ' 1 ~ 1 . 

Jo 

6. Express the product 2.5.8 . . . (3n — 1) in terms of the Gamma func¬ 
tion and generalize the result. [An3. 3 n r(n + §)/T(§).] 

7. Provo tho duplication formula for the Gamma function by integrating the 
(2n — l)th power of the identity, sin20 = 2sin0 cosO. 

8. Tho sum of two positive numbers is unity. Calculate tho root mean 6quare 
of their product. [Ans. I/V30.] 

9. By means of the substitution t = (x — 6)/(a — b) prove that 



— b)»~\a - x)*~'dx= B(p, q ). (a -by h«-i. 


10. f* cos”0cf0 

Jo 



sin n 040 = 


V7tr(jn + l) 

2 r(|n+ 1 ) ' 


Verify that when « is an integer, odd or even, these correspond to the " reduction 
formulae ” of the calculus. 


II. If n is a positive non-integer, prove that tho sign of T(—?») is ( — ) r 
where r is the integer next greater than n. 



cos 0d0 = 


(2^) s/a 


This result is important since it links tho Gamma function with the Elliptic 
functions; see Whittaker and Watson, Modem Analysis, p. 524. 



CHAPTER II 


Differential Equations 

21. The linear equation. 

The next step in our progress towards the study of the Bessel 
functions is to survey part of the field of ordinary linear differential 
equations. It is presumed that the reader already has some slight 
knowledge of the subject and has at least a nodding acquaintance with 
the standard elementary methods of solution. The object of the present 
chapter, which is mainly theoretical, is to refresh his memory and to 
fill in what are probably a few gaps in his equipment. The justification 
for its inclusion is that we shall ultimately be dealing with differential 
equations. 

The general linear equation may be written 

(1) VV + + ry" + . . . = X, 

where X and the coefficients p, q, r, &c., are in general functions of the 
independent variable x. The order of the equation is the order of the 
highest derivative present. The equation is said to be linear as there 
are no products and powers (other than the first) of the dependent 
variable y and its derivatives. If X be replaced by zero, the resulting 
equation is 

TV + yy' + ry" -f- . . . = 0. 

This is sometimes called the “auxiliary equation”; alternatively, 
“ the reduced equation ”. 

It is a property of the reduced linear equation that if y v y 2 , 2 / 3 * & c -» 
are solutions, so also is Ay x -p + Ctya + • • • "'here A, B, C . . . 

are arbitrary constants. The proof is simple; for if y x is a solution we 

have 

vy\ + yy\ + r yi' + • • • = °- 

py 2 + yyi + r y% + • • • = °> 

vy* + m + r y* + • • • = °> 

14 


Similarly 
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and so on. If we multiply the respective equations by A, B, C . . . and 
add, putting 

Ay\ + By* + Cy 3 = ¥, 

we have pY -+- qY' ~h rY" -f- . . . = 0. 

The result proves that Y is a solution of the reduced equation. The 
arbitrary constants are in practice chosen to suit the assigned con¬ 
ditions of a problem. 

Suppose that y v y 2 , &c., are solutions of the reduced equation. 
They are said to be linearly connected if constants a, b, c, &c., exist, 
not all zero, such that ay Y -f by 2 -f- . . . = 0. If no such constants 
exist, the solutions are said to be linearly independent. It is a pro¬ 
perty of the reduced equation that the number of linearly independent 
solutions is equal to the order. If the order is n and the independent 
solutions are y 1 to y n , then the most general solution is Ay x + By 2 + . . . 
-f By n , where the coefficients are arbitrary constants. This solution 
is known as the complementary function; it contains arbitrary con¬ 
stants, in number equal to the order. 

A solution of (I) other than a form of the complementary function 
is called a particular integral. The full solution of (I) is the sum of the 
complementary function and the particular integral. It is easily 
proved that this sum is a solution; for if Y is the complementary 
function we have 

pY + qY' + rY " + . . . =0, 
and if Z is the particular integral we have 

pZ + qZ’ + rZ" + ... = X. 

Addition gives 

P(Y + Z) + q(Y' + Z') +- X, 

which proves that Y + Z is a solution of (I). 


2 2 . Constant coefficients. 

The reader is probably familiar with the case where the coefficients 

f’J’ r ' • • m are constants. The reduced equation is then 

solved by the substitution y = e™ which leads to 

(p + jm + m*-f . . .)e"« = 0. 

The bracket has a number of zeros in, ?>» or ,j , 

mg solutions are *»,« &c so thTt tbL*’ ’ ’ ! correspond- 

, > oec.y so that the complementary function is 

Y = + Be™** -f. 
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There are rules for dealing with repeated values of m, as also for con¬ 
jugate complex or imaginary values. There are other rules for deter¬ 
mining the particular integral in the commonly occurring cases, such 
as when X is polynomial, exponential or trigonometrical. These tilings 
are not our immediate concern; they can be revised in texts specially 
devoted to the subject. 

2-3. Variable coefficients. 

When the coefficients in 2-1(1) are no longer constants the outlook 
is almost completely changed. There are no longer any rules for finding 
the particular integral, and if one meets an equation which is not 
reduced it is usually a matter of trial or guesswork to get the particular 
integral. Fortunately, most of the equations which one encounters are 
reduced and the need for finding a particular integral does not arise. 

As for the reduced equation, there is no golden rule for solving it. 
Relatively few equations are soluble in finite terms, and in fact the 
vast majority of differential equations remain insoluble. A number of 
methods are available; they include solution in series, solution by 
definite integrals and by contour integration. We shall return to the 
first two methods later. It should be borne in mind that in general a 
differential equation defines a transcendental function, and the study 
of such a function is not altogether precluded by our inability to solve 
the equation in finite terms. In any case a number of quite general 
propositions can be enunciated; they include: 

(i) If a solution of the reduced equation is known, the order can 
be lowered. In practice the known solution is often found by trial or 

inspection. 

(ii) The number of linearly independent solutions of the reduced 

equation cannot exceed the order. 

(iii) Linearly independent solutions of the reduced equation exist 

in number equal to the order. 

Further comment will be confined to equations of the second order. 
2-4. Second order equations. 

As a matter of convenience we take the linear equation of the 
second order in the form 

y" + y'f(x) + yg(x) = h(x). 

We now propose to establish the first proposition, that “ if a solution 
of the reduced equation is known, the order can be lowered ”. The 





DIFFERENTIAL EQUATIONS 


17 


proposition is true for any order. In the case of the second order the 
lowering is to the first order, so that usually the equation is completely 
soluble in these circumstances. 

We accordingly suppose that u(x) is a solution of the reduced 
equation, so that 

u” -f- u'f + ug = 0. 

The substitution y = uv implies 

y ' = u'v + uv\ 
y" = u"v + 2m V + nv\ 

The original equation thus transforms to 

v(.u" + u'f + ug) -f y'(2 u' -f- uf) + uo" = h. 

The first bracket on the left is zero by hypothesis. Accordingly, if we 
put v\x) = t(pc) wc can write 



As this is a first order equation in /, our statement is verified as to 
lowering. The value of t can be found by the use of an integrating 
factor, so that v and y can be found in succession. 

Example .—Consider the equation 

y" + xy’ — tj = ax 2 . 

The reduced equation 

if + xy" — y = 0 

is easily seen to have tho solution y = x, hence we make tho substitution y = xv. 
The original equation now takes the form 



The integrating factor is 

O = exp J(x + <lx = 

This enables us to write , 

v’O + v'Q' = axO = ~ (v'G), 
v' = Q~ l f axGdx + = a(l — ^ 


whence 



18 


APPLIED BESSEL FUNCTIONS 


Here b is an arbitrary constant; and a further integration gives v, so that 

y = xv — afx 2 + 2) -f car -f bx$x~ 2 e~l x, dx. 

The second arbitrary constant is c; the integral in the last term cannot be evalu- 
ated in finite terms, which is in line with our remark that in general a differential 
equation defines a transcendental function. It will be observed that the first 
term on the right corresponds to the particular values b = 0 = c and it is left to 
the reader to verify that it is the particular integral and really is a solution. The 
second term is our original guess at a solution of the reduced equation. 

We now come to our second proposition. It takes the form: “ The 
reduced linear equation of the second order cannot have more than 
two linearly independent solutions ”. It is established by showing that 
a contrary assumption leads to a contradiction. 

We may remind the reader that two simple simultaneous equations 
may be inconsistent. Thus the two equations 

4x + 6y + 7 0 = 6x + 9y + 5 

cannot be simultaneously true, since (6x -f- 9 y) is fifty per cent greater 
than (4x -f Gy). If we remedy the deficiency by replacing 5 by 104, the 
equations become consistent but are effectively one and the same. 
In this case there is a value of y for every arbitrarily assigned value 
of x and the number of solutions is infinite. 

Reverting to our reduced differential equation, 

(1) y" + y'/(*) + yg( x ) = °> 

we assume that there can be three linearly independent solutions, 
which we may denote by y v y 2 and y 3 . We then have 

(2) y/' + yj + ytf = o, 

(3) y 2 " + yj + ytf = 

(4) y 3 " + y 3 '/ +'M 7 = °- 

Multiply the first by A, the second by /x and add all three. It must be 
possible to assign A, /x so that the coefficients of/and g become zero, 

i.e. 

(5) A y x + /xy 2 + y 3 = °> 

(G) Ay/ + y-yi + y/ — °- 

To deny that A, /x can thus be found is to assert that 

y/ _ y% 
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The integration of this small differential equation leads to 

logy! = logy 2 + Iogc; y x = cy 2 . 

This violates the assumption that the solutions are linearly indepen¬ 
dent; hence A and \l can be found. But note that so far there is no 
justification for thinking them to be constants. 

We are also in possession of the further equation 

(7) Ay x " -j- /xy 2 " -f y 3 " = 0, 

which follows automatically from the others. The differentiation of 
equation (5) gives 

A Vi 4- A'yi 4- mi 4- /*'y 2 + yi = 0 , 

so that on subtracting (6) we have 

A'yi 4- ^'y 2 = o. 

Similarly from (6) and (7) we derive 

a V + v-'yi = o. 

In the last two equations the possibility 



'J\ y-i 


is ruled out; hence we must have 

A' = 0 = p’, 

so that A and n are absolute constants. This establishes that the three 
solutions are in fact connected by a linear relation 

Ayi 4- m-i 4- y 3 = o. 

It is one thing to say that there cannot be more than two indepen¬ 
dent solutions, but it is quite a different matter to prove that two 
independent solutions actually exist. A proposition of the latter type 
is known as an “ existence theorem Existence theorems usually 
have three properties; they are long, dull and difficult. It takes a 
special type of mind, even a special type of mathematical mind, to 
revel m existence theorems; and as the reader is presumed not to be 
a mathematical specialist he may be condoned if he accepts the pro¬ 
position on higher authority. We shall accordingly assume that, except 
perhaps in certain artificial and factitious cases, the second order 



20 


APPLIED BESSEL FUNCTIONS 


equation lias two independent solutions which are continuous over 
certain ranges limited by the infinities of the coefficients. 

Parenthetically, let y x , y 2 denote the two independent solutions; 
then the general solution is y = Ay x -f- By 2 where A, B are arbitrary 
constants that can be chosen to fit assigned conditions. This raises the 
question of what two solutions are fundamental. The answer is purely 
a matter of choice or convenience. Consider the oscillation equation 
x -j- oj-x = 0 for simple harmonic motion. It has the solutions 



r*"*, cos tut, sincuJ, 4 sin (ujt -f- 5 cos (ojt — ^ 


and so on. None of these is more fundamental than any other, and a 
linear combination of any two is a solution. With other differential 
equations there is usually some property that makes one form more 
desirable than another, e.g. it may vanish at the origin, or be more 
easily calculable. It results that occasionally mathematicians differ 
as to what are the best forms to employ, and the reader must be pre¬ 
pared to encounter this diversity. 


2-5. Properties of the independent solutions. 

It has already been mentioned that our inability to solve a par¬ 
ticular equation does not wholly deter us from studying the function 
which it defines. The second order equation, even when insoluble, can 
be made to give information of two types. The first concerns relations 
between the two solutions and the second concerns properties of their 

zeros. We begin with the former. 

Suppose that y x , y 2 are two independent solutions of equation 
2-4(1), so that equation 2-4(2) and 2-4(3) hold. On eliminating g we 
have 

( yi”y2 ~ Vivi ) + (yit/i — ya/z )f = 

Note that the first bracket is the derivative of the second, so that if 
we denote the latter by z we can write z' + zf = 0. The integrating 
factor is G = exp jf(x)dx, so that 

L (zG) = 0, z = AG -1 , 
dx 

or 

(1) ViVz ~ VxVz = A ex P {“ 

an extremely useful relation. 



DIFFERENTIAL EQUATIONS 


21 


Example 1.—Consider the equation 

(1 - r - 2x/ + »(n + 1)// = 0. 

It is known as Legendre's equation, and the two transcendental functions which 
are its linearly independent solutions corresponding to our y,, y, arc denoted by 
P n and Q n . Here we have 

— 2x A 

/<*>“ r-^ = rfi lo « (1 ~ x) ‘ 

Hence G = (1 — x 2 ) and the relation becomes 

PnQn - Pn'Qn = {1 

Here A is an absolute constant and its actual value will depend on the forms 
adopted for P n and Q n . 

The relation (1) can be carried a step further; for we have 

y\!k ~ _A_ = d_ /.v, 

Vz Gy* llx \y-i. 

Hence on integration, 

( 2 ) 


y i 


- + A f& 


It follows that if y 2 is known, y x can be determined, as stated 
previously. 


we have 


Example 2.—Using, as before, the reduced equation 

if + xy' — y = 0, 

/(*) = x, $f(x)dx = Jx 2 , 
whence G = el* 1 . As y 2 = x wo derive 

l/i = x(B + >lJar a e“l z Vx), 
which agrees with our previous result. 

2*6. Zeros of the solutions. 

We come now to the second type of information. It concerns the 
zeros of the functions, the values of * for which the function vanishes 
or the places where the graph crosses the x-axis. And lest the reader 

We j 1 *? ? n these zeros > ifc may mollify him to know that 

they usually hold the key to any problem under discussion. 



22 


APPLIED BESSEL FUNCTIONS 


If we take the reduced equation in a form where the coefficients are 
devoid of fractions we can write 

(l) y"p(x) + y'q(z) 4- yr(x) = o. 

The statement that the general solution At/ 1 + By. z can be made to 
fit two conditions has one curious limitation which we shall later find 
useful. It may be stated thus: No solution of (1) can touch the x- 
axis except at a zero of p(x). 

We are assuming that p, q and r are finite and have derivatives of 
all orders. If the solution touches the x-axis we have y = 0 = y', 
whence, if p(x) is not zero, y" = 0 by the equation. Differentiation 
gives 

y'"v + y'\v + q) + y'W + *) + = °> 

so that if p(x) is not zero it follows that y'" = 0 and all the succeeding 
derivatives can be similarly shown to vanish. The application of 
Taylor's theorem to the calculation of the neighbouring ordinate 
gives 

y =/(* + ft) =f(x) + hf'(x) + k - f"(x) + • . . • 

Since every term in the series vanishes by hypothesis, the ordinate 
must vanish identically and the curve degenerates to the x-axis. 

As an example, revert to Legendre's equation, where p(x) = (1 x 2 ). 

No solution can touch the x-axis except at the two points given by 
x 2 = 1. It is a waste of time seeking a solution that touches elsewhere. 

2-7. The normal form. 

In the study of algebraic equations it is often an advantage to 
remove certain terms. It is similarly found advantageous in certain 
theoretical investigations to remove the middle term from the reduced 
second order differential equation. The result is known as the normal 
form and it can be reached by the substitution y = uv. Taking the 
equation in the form 

y" + y'M + yg(x) = 0, 

this gives as before 

v(u " + u 'f + ug) + v\2u' + «/) + = 0. 

The coefficient of v' vanishes if we choose u so that 
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2 u + «/ = 0. ^ + \f = °> l°g « + i ffU) <1* = 0; 

U = exp{— i//(x)<£r}, u' = — £«/, m" = — luf -f ju / 2 
The equation becomes 


rfx 2 


where 


+ vl = 0, 


if 


I(x) = 9 — h ~ U 2 > v = -■ 

ax u 


It will be observed that the normal form is more comprehensive than 
the original w hence it derives. Any number of equations may have the 
same normal form. For a given I we can assign/ arbitrarily and then 
g follows automatically. 


Example .—As an illustration, consider the reduced equation 

y" + * 1 / — y =0. 

This has 

? = /=*. u = exp(— i-r*), 1 = -f G)/-L 

Hence the normal form is 


d 1 v 

fak = M** + 6), V— y exp (I* 2 ). 


It may be added that the normal form is rarely an aid to solution; 

its utility lies in other directions, principally in connexion with the 

zeros. The following is a fair sample of the type of argument. If a and 

6 be two permissible limits of integration, we can integrate the normal 
form 

V + vl = 0 


to g‘ ve [v'j+fjvldx = 0. 

If in particular a and 6 are two consecutive positive zeros of the deri¬ 
vative v we have 1 



in ,mph r tl0n bem ® tlat the “t*g™»d carnot maintain 

Zed aln i“T f f “■ " “ Were further k »own ft* / main 
tamed sign m this interval, we should be driven to the conclusion that 
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v must change sign. In other words, v must in these circumstances 
have a zero between two consecutive zeros of its derivative. All of 
which, of course, depends on the derivative having two zeros; it may 
have none. The converse theorem, known as Rolle's theorem, is en¬ 
countered early in one’s study of the calculus. 

Example. —Consider the equation 

\f -f 2/ sinx -f j/(l -f sin*x) = 0. 

The hope of solving it is rather remote; but it is easily shown to have the normal 
form 

v” + r(l — cos-r) = 0. 

Here I does not change sign in any range. It follows that if any of the infinite 
number of solutions of the equation for v has maxima and minima, these must 
alternately lie on opposite sides of the x-axis. We also have 

u = e.\p[— Ij2 sinxdx] = exp(cosx), 

which vanishes for no value of x. As y = uv, the zeros of y are governed by the 
zeros of v. 

2-8. No common zero. 

Two independent solutions of the normal form cannot have a 
common zero, or their graphs can never cross the x-axis at the same 
place. We establish this interesting and important result by assuming 
that v x and t’ 2 are independent solutions of 

v" + vl = 0. 

From the equations 

i’i" + vj = 0 = v 2 " + rj 

we deduce 

Vy'v 2 ~ v \ v i' = 0 = — ( V 1 V 2 — V l V 2 )» 

whence 

Vi v 2 ~ V l V 2 = A • 

If Vj and r 2 can be simultaneously zero, then the constant A must be 
zero. In that case we could write 

< = 
v 1 v 2 

whence vjv 2 = const., which violates the assumption that v v v 2 are 
independent. We conclude that t’ 2 have no zero in common. 



DIFFERENTIAL EQUATIONS 


2o 


2-9. Interlacing of zeros. 


We can now prove the interesting and important theorem that the 
zeros of t’j and v 2 interlace. Having already established that 


we can write 


whence 



Choose a, b to be two consecutive zeros of v v if such exist. We already 
know that v 2 cannot vanish at either end-point; let us further assume 
that it nowhere vanishes in the range a > x > b. Then v x fv 2 is a con¬ 
tinuous function that vanishes at both end-points, so that the value 
of the above integral is zero as calculated from the left. But this 
establishes a contradiction, for the integrand does not even change 
sign in the range. We conclude that v 1 /v 2 is not a continuous function; 
in other words, v 2 has a zero between the consecutive zeros of v 1 at 
a and b. It is immaterial which of the solutions we call v, and v 2 , so 
that each function must have a zero between two consecutive zeros 
of the other: the zeros of either occur alternately. 

In case the reader finds the matter abstruse, let him consider the 
very ordinary normal form 


It has the two independent solutions sincor, coscox. These evidently 
have no common zero; and their zeros interlace. For the matter of 
that, the same is true of the two independent solutions 5 sin (cox — 

7 cos (cox -{- |7r). 3 


2-10. A comparison theorem. 

Before proceeding further it is necessary to recall a mean-value 
theorem from the integral calculus. It concerns the integral of the 
? r “ d “ C f two hnctums and shows how this might he related to the 

of vSae, T ‘a* Ct ' 0DS “ l0ne - W ° rkin 8 in <M»ite range 
of values a > x > b, we stipulate that the first function f(x) shall be 

one-va ued and continuous. Somewhere in the range is a value of 
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which we may call where f(x) takes the value/(£). With the second 
function <b(i) we are not concerned about its continuity; but we 
insist that it shall be one-valued and not change sign in the range. 
The theorem states that 

f b f(x)<f>(x)dx =/(£)^ <f>(x)dx, a > £>b. 

The proof is simple, for the continuity of f(x) guarantees that it achieves 
its upper bound U and its lower bound L. As <f>(x) does not change 
sign, the value of the integral on the left must lie between 

Uf <j>(x)dx and L,f <f>(x)dx. 

J b J b 

Hence there is a number K between U and L for which we can write 


J a f(x)<f>{x)dx = K.J <f>(x) dx, 


U > K> L. 


The continuity of f(x) further guarantees that it takes every value 
between U and L at least once. Hence there must be a value of x, 
which we have called £, for which/(£) = K. This establishes the result. 

We now revert to our differential equations. The normal form be¬ 
comes most fruitful when comparison is made with a similar and 
simpler form. Suppose that we have two normal forms 

u" + vH = 0 = v" + vK. 


From these we derive 

uv" — u"v = (H — K)uv, 

whence by integration 

[mu' — m'u] (H — K)uvdx. 

Let a, b be two consecutive positive zeros of u; there is no loss of 
generality in assuming that u is positive when a > x > 6. This makes 
v' positive at the left end and negative at the right, as a rough sketch 
will show. We now have 

u(a)v'(a) — u(b)v'(b) = f(H — K)uvdx, 

J b 

and we propose to show that if H > K throughout the range a > 
x > b then m must have a zero between a and b. For since (H K)v 



DIFFERENTIAL EQUATIONS 


27 


is presumed not to change sign, being in fact positive, and u{x) is pre¬ 
sumed to be one-valued and continuous, we can write 

u{a)v'{a) — u(b)v(b) = t/(£)— K)vdx, a > £ > b. 

Inspection shows that u(a), u{b) and u(£) cannot all have the same 
sign, which proves that w(z) must somewhere cross the x-axis and have 
a zero in the range a > x > b. 

As a corollary, if we assume that u and v have a common zero at 
x = b, the relation runs 

u[a)v'(a) = u (£)]\h - K)vdx. 

As the integrand is positive and v'(o) is negative, it follows that u(a) 
and u(£) are of opposite sign. Hence n{x) vanishes a second time before 
v(x) has time to vanish at x = a. 

Example.— Consider the equation + x~'e*>j = 0, of which the solution is 
by no means obvious; in fact, its solution is not easy to obtain. Away from the 
origin, the coefficient x~ l e* is continuous and has derivatives of all orders. Com¬ 
pare this with the equation xf + nhj = 0, the zeros of whose solutions occur at 
intervals of n/n. Whatever the value of n, it must eventually happen for suffi¬ 
ciently large values of x that exceeds n*. so that (ar'e* - »«) is positive. 
It follows that any solution of the given equation must ultimately have a zero 
in every range whose length is w/n. With increasing x we can afford to make » 
larger in the comparison equation, whence we conclude that the interval between 

consecutive zeros of solutions of the given equation probably decreases with 
increasing z. 


2 - 11 . Equations containing a parameter. 

Having expounded aU that we shall need to know of the equation 
in its normal form, we turn to another matter that will figure pro¬ 
minently in the ensuing pages. It sometimes happens that a differ- 

tion a ha qUa i 10n d C ° n i! a,n8 “ parameter »: case of Legendres cqna- 

t.on has already been mentioned as also the oscillation equation 

?n vTl * °j- Th fi J nctl0n of x def ined by the equation then differs 

h JthTn ft CCOrd, T 8 t0 tHe VdUe ° f ”• The symbo1 denoting the function 

As example, the eqnat.on a*" + (1 - x) ,f + ny = P 0 is knQwn 
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as Laguerre’s equation of order n and its solution is denoted by L n (x). 
Its recurrence formula is known to be 

L„ +1 (x) — (2 n + 1 — x)L n (x) + n 2 L„_ 1 (x) = 0, 

and this connects the functions of the three contiguous orders n + 1, 
n and n — 1. In anv given case the recurrence formula can be derived 
only when the solution of the differential equation is known. Con¬ 
versely, the differential equation cannot be derived from the recur¬ 
rence formula alone. Usually there is a connecting differential relation; 
in the present case it happens to be 

L n '(x) — nL' n _ x (x) = —nL n _ x {x). 

212. Orthogonal functions. 

Two functions which merely differ in their order usually possess 
interesting integral properties, one such property in particular being 
known as the orthogonal property. We can illustrate this by Legendre's 
equation. Let P m (x), P n (x ) be the functions that respectively satisfy 
the equations 

(1 - x 2 )P m " - 2 xP m ' + m(m + 1 )P m = 0. 

(1 - x 2 )P n " - 2xP + n(n + 1 )P n = 0. 

On multiplying these respectively by P„, P, n we deduce 

+ l) _ n(n + l)}P„P m = ~ {(1 - ^)(P n 'P m - P„P„')}- 

A bracket on the right suggests limits of integration and we have 

j_Pm^)P n ^) <lx = 0 . m * n - 

This is the orthogonal property, and the functions P v P 2 • • • » re 
said to form an orthogonal family. As a matter of fact the reader has 
for years been familiar with at least two orthogonal families. The one 
is sinx, sin2x, sin3x . . which has the property 

J smmxe'mnxdx = 0, m =t= n. 

The other is obtained on replacing sin by cos. 
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EXERCISES 

1. Solve Laguerre’s equation of order unity, viz. 

+ (i — *)/ + y= o 9 

after first verifying that y = ! — x is a solution. 

2. Prove that any two solutions of the equation 

cosx + y' sinx + y sin 2 x — 0 
are connected by the relation 

y\yi — Vi'yi = «cosr. 

3. Verify that e~ x is a solution of 

jT(cosx — sinx) + 2i/ cosx -f y(cosx + sinx) = 0, 

and hence deduce the other solution. In what circumstances could a solution 
touch the x-axis? 

4. Prove that the last equation has the normal form v* — v * + s * n * x = q 

1 — sin2x 

5. Prove that the equation of damped oscillations, x + 2nx -f- ( a 2 + b 2 )x = 0 
has the same normal form as when the oscillations are undamped. 

6. Reduce the equation ^ (1 — = 0 to normal form. What 

do you conclude from the result ? 

7. Prove that any solution of *y + xxf + (** _ n*)y = 0 must have an 
infinity of zeros. 

8 - Any f olu tion ofxy* + 2/ -f- (1 -f- x)y = 0 has an infinity of zeros whoso 
interval ultimately tends to tt. 

9. Without using the normal form, prove that two independent solutions of 
the reduced equation cannot have a common zero. 

10. If U is greater than K in the two normal forms 

«' + uH = 0 = if + vK, 

r;irr."E;-.2™s.“ ^■*-*-*'-=•■ 

11. Weber's function D„(x) has the recurrence formula 

D n +1 — xD n + 7iD n _, = o 
and satisfies the differential relation 

D n + \*D n — nZ>„_, = 0. 

Deduce that it is a solution of the differential equation 

*lf + y(4n + 2 — X 2 ) = o. 
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12. Deduce the recurrence formula for Laguerre's function from the differ¬ 
ential equation 

ry' H- (1 — x)y' -f ny = 0 
and the differential relation 

7 l ^n-l /== “V 

13. Prove that Laguerre’s function has the orthogonal property 


f e-*L„L m dx = 0. 
JQ 


14. The Tschebyscheff function T n of order n satisfies the equation 

(1 — x~)y" — Ty* + n 2 y = 0. 

Prove that it has the orthogonal property 

n TJT„ . 


/ I 1 T 

* - «• 


15. The Hermitc function //„(x) of order n satisfies the differential relation 

Hn = ZnHn-i 

and has the recurrence formula 

II n+l - 2 xH n + 2 nH n _ x = 0. 

By differentiating the latter, show that H n satisfies the differential equation 

y’ — 2 xy' + 2 ny = 0. 

Deduce the orthogonal property 


f //„ H m r*'dx = 0. 

J — CO 


The Hermite and Laguerre functions occur in the applications of Schrod- 
inger’s wave equation in quantum theory; see H. Weyl, The Theory of Groups 
and Quantum Mechanics , pp. 54-70; B. L. van der Waerden, Die gruppen - 
theoretische Melhode in der Quanlenmtchanik , pp. 12—1G. 



CHAPTER III 


Cylinder Functions 


31. Recurrence formulae. 

A cylinder function C n (x) may be defined as a function of x involving 
a parameter n and satisfying the two recurrence formulae 

<?„-,(*) + C n+1 (x) = — C„(x), 

X 

C-i(*) C+ifc) = 2 C( x )* 

In general x and n are unrestricted; but in this work attention will 
be confined to such values of x as are real and positive, with n real 
but not necessarily integral or positive. 

With a view to simplifying the notation, the argument x will be 
dropped when not in doubt, so that we shall write C n for C„{x ); and 

a prime, or dash, will be used to denote differentiation with respect 
to the argument, so that 


<V = ^<* ); <V<„x) = * CM = 1* <7„(ox). 


^ ' d(ax) nV 7 adx 

The above two formuta are accordingly written 

(1) C„_, + C„„ = — C„, 

X 

(2) C-I - C+1 = 2C n '. 

By addition and subtraction they are equivalent to 

(3) *C-i = nC n + xC n \ 

(4) *C +1 = «C - xC n \ 

mmmmm 
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six different types of Bessel function, each with its own peculiarities 
and bearing the name of some eminent mathematician, yet all pos¬ 
sessing properties in common, there is a certain advantage in retaining 
the term cylinder function for the moment. We can discriminate 
between them later on. 

Much of mathematical physics is dominated by Laplace's equation 
and Laplace’s operator, and when these are translated into cylindrical 
co-ordinates, Bessel’s equation almost inevitably appears, as will be 
found more than once in the ensuing pages. This accounts for the 
name “ cylinder function ” which is common on the Continent and in 
America. It appears to be due to Heine. 

3-2. Bessel’s equation. 

The differentiation of 3-1(3) gives 

(1) xC n " + (» + 1KV = xCn-x + C u . v 

Multiply 3-1(3) by n and subtract from (1) multiplied by x. This gives 

a*7 n " + xC n ' - n 2 C n = x{xC - (n - l)C n _,} = - x*C n , 

the last step being justified by 3-1(4). The result shows that C n satisfies 
the equation 

(2) xty" + xy' + (x* - n 2 )y = 0, 

which is Bessel’s equation of order n. 

The recurrence formula) cannot be directly deduced from Bessel’s 
equation and it is irrational to conclude that any particular Bessel 
function necessarily satisfies the recurrence formulas. This happens 
to be true when the arbitrary constants in the solution of (2) are 
properly chosen as functions of n; but it cannot be substantiated until 
Bessel’s equation has been solved. For the solutions that will be 
adopted, it will prove easy to verify that 3-1(3) is satisfied, and the 
rest will follow\ For, reversing the above argument, a function that 
satisfies (2) and 3-1(3) also satisfies (1) and hence it satisfies 3-1(4). 

3-3. Interlacing of zeros. 

The recurrence formula) rapidly lead to important deductions. 
Thus 3-1(3) is equivalent to 

(1) x n C n .. x = nx”-'C n + x n C n ' = ~ 
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According to Rolle’s theorem, between two consecutive zeros of a con¬ 
tinuous function lies a zero of its derivative. Taking the continuous 
function to be x n C n , between two of its consecutive zeros lies a zero 
of its derivative, i.e. of x n C n _ v If we move away from the origin to 
avoid possible discontinuities or such zeros as may come from x n , we 
conclude that between two consecutive positive zeros of C n lies a zero 
of C„_y 

Similarly 3-1(4) can be written 

(2) x-C n ' - nx--'C„ = -z-C M = L {r-”C„}, 

whence it follows that, except possibly near the origin, a zero of C„ + , 
lies between consecutive zeros of C n . Changing the order from n to 
n — 1, a zero of C n lies between consecutive zeros of C„_ 1 , and taking 
this in conjunction with the similar result obtained from (1) we con¬ 
clude that the zeros of C n and C n _ 1 interlace. Incidentally this is not 
to be confused with the similar proposition, proved in the last chapter, 
concerning two solutions of the same equation. 

Note that (1) can be reversed to give 

(3) J x n C n _ 1 (x) dx = x n C n (x) 
and (2) can be written 

( 4 ) fx~"C n+l (x)dx = —x~"C n (x). 

3-4. Integral orders. 

It seems at first sight reasonable to suppose that functions of 
integral order might be more tractable than others; but they possess 
one awkward drawback. Putting n = 0 in the recurrence formula) 
3-1(3) and 3-1(4) we have 

(!) C 0 ' = C_ x = -C v 

It appears that, of the two functions C x and C_ x of the first order, the 
one is a mere numerical multiple of the other. We proceed to inves¬ 
tigate further. Putting n = 1, —1 i n 3-1(1) we have 

xC 0 + xC 2 = 2 C lt 

u , . *0-2 + XC 0 = -2C_ V . ' 

13y subtraction, 


whence 


*( c 2 - C- z ) = 2(C, + C_,> = 0, 
c -2 = (—l)*C7 a . 
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Putting n = 2, —2 in 3-1(1) we have 

xC, + xC 3 = iC 2 , 
xC_ 3 -f xC_ 1 = — 4C_ 2 . 

By addition 

x(C 3 -f- C_ 3 ) -f- x(C 1 + C-i) = 4(C 2 — C_ 2 ), 

whence 

C-3 = (—1) 3 C' 3 . 

In general, 

± <?_„_!) + x(C n _ x ± C. n+1 ) = 2n{C„ T C„), 

according as n is even or odd, and it is an easy induction, which can 
be left to the reader, to prove that 

(2) C_ n = (—l)"C n . 


Bessel’s ecjuation, being of the second order, certainly possesses 
two linearly independent solutions. The important point is that for 
integral orders the second independent solution is not obtained by 
merely changing the sign of n. This complicates matters considerably, 
and later in the book we shall take up the problem of finding the 
second solution in such cases. 


3-5. The normal form and the zeros. 

The ideas expounded in the previous chapter can now be applied 
to Bessel’s equation. Its leading coefficient is x 2 , whence we conclude 
that no solution can have a repeated zero except possibly at the origin. 
Otherwise expressed, C n can have no repeated positive zero and the 
graph cannot touch the x-axis except possibly at the origin. Com¬ 
paring Bessel’s equation with 

y" + y'A x ) + yg( x ) = °> 

we have, in 2-7, 

f( x ) = -> o( x ) = 1 — \ 

X 


u = exp{— \Jf(x)dx} = x *. 


4x 2 


( 1 ) 


t’" + 


< 


1 + 


1 — 4 n 2 


4x 2 


^ = 0, v = x*C n (x), 


This is the normal form of Bessel’s equation. 



CYLINDER FUNCTIONS 


35 


In the particular case where n = ±2 the normal form reduces to 
the simple oscillation equation v" + v = 0 with the solution v — R 
sin(x + a). This connects Bessel functions with trigonometrical func¬ 
tions and we have the remarkable result that x*Cj and x k C_± are of 
the form a cosx -f- b sinx with an infinite number of zeros at intervals 

Of 77 . 

More generally, I is necessarily positive if (i) 4n 2 <l; or (ii) 
4xr -\-1 > 4 n 2 , which for any value of n must be the case ultimately 
for large x. If n 2 < I > 1 and is monotonic decreasing to unity. 
If n 2 > 7 < 1 and is monotonic increasing to unity. 

We now propose to show that if n 2 < £ the function x*C* n has at 
least one zero in any range 0<a<i<a + s. The proof consists in 
showing that the contrary assumption leads to a contradiction. Assum¬ 
ing then that the function has no such zero, there is no loss of generality 
in taking it to be positive throughout the range. The origin is de¬ 
liberately avoided since we do not know how the function behaves 
there. 

We make comparison with the function to = sin (x — a), which 
has consecutive zeros at a and a + n. Between these values it is posi¬ 
tive and it satisfies the equation w" jo = 0. On the left, jo is posi¬ 
tive at a, whilst it is negative at a -f- tt on the right. From the norma) 
form v" -f- vl = 0 we deduce 

vw" — v"w = (7 — l)vW, 

whence 


+ * f" Tfl + ff r- -ia+f 

J (I — l)i nvdx = \yw' — v'wj = ^yw’j 


= v(a + ir)w'{a -f- tt) — u(a)u?(a). 

In accordance with our assumptions the right side is essentially nega¬ 
tive; but the integrand is nowhere negative in the range. This estab¬ 
lishes the contradiction and we conclude that x*C„ has a zero in the 
given range, and hence C n has an infinite number of positive zeros at 
intervals of not more than it, provided n 2 < 

An upper limit to the interval between the zeros being tt when 

f < we 01111 now establish a lower Umit. Let a be a zero of C n and 
161 / 


Ha) = 1 + 



Then 7(x) < h 2 if x > a. Consider the function w = sin*(x — a) 
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which has consecutive zeros at a, a -f ~/k and satisfies the equation 
w" + kho = 0. Then 

J (k 2 — l)mcdx = J (icv" — w"v)dx = [jet/ — w’v J . 

Take P to be the zero of v next following a and suppose that, if possible, 
it falls short of a + t r/k. This ensures that (ex), ic(x) and k 2 — I(x) 
are all positive (or at least, not negative) between a and P. The integral 
is therefore positive; moreover, 

v((S) = 0 = v(a) — w(a). 

The right side accordingly reduces to ic(P)v’(P) where v'(P) is negative. 
This makes ic(P) negative, thus showing that w has crossed the x-axis 
and contradicting the assumption that P fell short of a -j- 7 r/k. We 
conclude that /? > a -f ir/k, or 



a ~ tt > P > a + 


77a 



Note that with increasing a the interval approaches 77 from below, 
provided n* < J. 

We are now in a position to state that any cylinder function, of 
whatever order, has an infinity of zeros. For choosing any value of n 
such that n 2 < £ and calling it p, we have established that C v has an 
infinity of zeros. These are interlaced by the zeros of C p+1 , and these in 
turn by the zeros of C p+2 , and so on. Working in the other direction, the 
zeros of C v interlace, and are interlaced by, the zeros of C p _j; and so on. 

We shall now establish a theorem concerning the second non- 
negative zero. If the origin happens to be a zero, the proposition 
concerns the magnitude of the first positive zero. Let be the 

first two non-negative zeros of C„, and let /J be the zero of C n between 
them. Obviously £ < 0 ^ and we propose to show’ that n < p. Consider 
the graph of C n between a t and a*. If it passes through aj, from below, 
C n is positive, the curve rises to a maximum and C n " is negative. 
Alternatively, if it passes through a 1 from above, C n is negative, the 
curve sinks to a minimum and C n " is positive. In either case C n and 
C n ‘" are of opposite sign. By hypothesis C n '(P) = 0 and if we put 
x = p in Bessel’s equation we have 

P 2 C n "(P) + (P 2 ~ n z )C„(P) = 0. 

We conclude that (p 2 — n 2 ) is positive and that p is greater than n in 
absolute value. Hence | n | < p < a*. Descriptively this means that 
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as the order rises the graph is more leisurely about reaching the second 
zero. 

It remains to consider the heights and depths to which the graph 
rises and falls between crossings of the x-axis. We shall prove that the 
successive maxima and minima decrease in absolute value. Taking 
the simple case of zero order, Bessel’s equation is 

C 0 " + - C 0 ' + C 0 = 0. 

X 


Multiplying by 2 C„' and integrating between unspecified limits a and 
6 , we have 

M .;+£ <v>*+[c.*];-°. 

If a and b are taken to be positions of stationary value, not necessarily 
consecutive, we have 

2 /‘ l <V* = <V(a) - c„ 2 (6). 

* a •h 

The integral is certainly positive if we are working on the right of the 

origin, and so C 0 (a ) is greater than C 0 (h) in absolute value. 

More generally, we can multiply Bessel’s equation of order n bv 
2C n and write 

2x2 C ""Cn' + -J X . C n ' 2 + 2 C n Cn = 0. 


x 2 — n 2 
As 

dx 

integration by parts gives 
tr 2 


a 2 — n 2 


- / ** \ _ — 2xn* 

& \ x2 — n 2 / (x 2 — n 2 ) 2 ’ 


[A °’i + xfe. + <v*+[>.■].*- o. 

L b ; positions of -*“* — 

f* 2X 3 

J. C "’ 2 ' fx = c » 2 (“> - C» a (4). 

"tJtzk„? :i work ° n **«- 
—*—— - - — g r — z 
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It is an easy deduction from this that there is only one zero of C n ' 
between two consecutive zeros of C n . 

To sum up what we have learned of the function C n (x), we have no 
formula for calculating it; we do not know its value for any particular 
value of x; nor can we say how it behaves at the origin. On the other 
hand it has an infinity of positive zeros which occur at nearly equal 



Fig. 4 .— Match of 3 cylinder functions of consecutive orders 
Cm = 0 when C* + i = Cn-1 (at A and B) 


intervals; between a consecutive pair is a turning value whose abso¬ 
lute magnitude steadily decreases, as also a zero of the next consecutive 
order both above and below. The function of order i has its zeros at 
absolutely equal intervals of tt and is expressible in sines and cosines. 
For the most part the graph looks like a damped vibration, but it can¬ 
not be called periodic; for one thing, there is no constant period, and 
for another thing, the values do not repeat themselves (fig. 4). 


3-6. Transformations. 

It rarely happens in practice that Bessel’s equation makes its 
appearance in the standard form which we have already given, and it 
is desirable that one should be in a position to recognize it under other 
guises. For instance, it is scarcely self-evident that xy" + y = 0 is 
a form of Bessel’s equation. 

One of the most obvious forms is obtained on replacing x by h r, so 
that dx is replaced by kdx. The equation becomes 



dx 2 x dx 



with the solution y = C n (kx). Other exercises are given at the end of 
the chapter and we append an example of the method. 
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Example .—We have, as at 3-3(1), 

j x U n c n ] = * n c n _ t = T[z"-«c n _,]. 

To find the equation satisfied by x n C n we have 

^ [^O n ] = x^C n _ t + 

= * n -'Cn-i + *[x n -*C„_ 2 ] 

= ^~ l C n -i + Cn-1 - C«] 

= ( 2 « - 

The elimination of C n _ x then gives 

y*+ '—^!/ + y=0 

X 

as the equation satisfied by y = *"C n . Note that if n = |, this reduces to the 
oscillation equation for simple harmonic motion and again proves that x*Ci 
can be expressed in terms of sin x and cos x. * 


A comprehensive form that covers most cases occurring in practice 

( 2) r + L= 2 “+ [,/)>*»-.,* + ^=;-vj y _ 0 _ 

whose solution is 

(3 > y = z a c n (/Sx?). 

We establish this result as follows. A solution of 



(Ft , dl 

r z *' hz dz + {f * 2z *~ n2)t = 0 


is i — C n (f3z). If 2 = x y we have, 
tion with respect to x. 


using dashes to denote differentia- 


dz z' ~~ vxy-E 


This can be written 
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whence we get the equivalence of the operators 

d _ 1 d 
2 dz y X dx 

when applied to t. The differential equation can be written 

2 £ ( 2 1 ) + 

whence 

1 i ( x s) + ~ nVl - 

dt 

The substitution t = yx~ a gives x — = y'x 1_a — ayx~ a , 

x L (* £) = y “ x ^ + (l ~ 2a )y' x '-' + «***-• 

After multiplying through by x a " 2 this gives the required result. This 
will be taken as a standard and will be referred to repeatedly in the 
following pages. 

Example .—As an example of its use, let us settle the question whether the 
equation x 2 y" — xy' + (xr x -f l)y = 0 is a form of Bessel's equation; and if so, 
what is the appropriate cylinder function. We recast it as 

»* - £ + (** + = °- 
Comparison with (2) then gives 

1 - 2a= -I, 0y = I, 

2y — 2=2, a *-nV=i; 

whence a = 1, [3=4, y = 2, n = + J. 

We conclude that the equation is of Bessel’s type and has the solutions 

xC^), xC- ,<**•). 

It is worth while making a few observations on equation (2). Note 
that if the middle term is missing, so that the equation is in its normal 
form, we must have a = This enables us to say at sight that if the 
middle term is missing, the solution must have the factor x*. Con¬ 
versely, a function of the form x i C n (fix y ) satisfies an equation in its 
normal form and y' is missing. Note further that y cannot be zero; 
otherwise expressed, the first term in the scpiare bracket cannot be 
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absent. For one thing, the argument of the function ceases to vary 
and degenerates to a mere constant. Moreover the equation ceases to 
be of Bessel’s type. It degenerates to 

*V" + (1 - 2 a)xy' + a 2 y = 0 

and is soluble by elementary means. The reader can verify that the 
results are consistent, a solution being y = x“. Incidentally it is illu¬ 
minating to find the other solution. 

A case that occurs often enough to merit special notice is when 

a 2 = n-y 2 , y = £ so that a = £ n . Its differential coefficients are par¬ 
ticularly simple, for if 

y = x^C^px *) 

the recurrence formula runs 


nC n + /&■»CV = pz*C n _ lt 

where the unwritten argument is px i and the dash denotes differen¬ 
tiation with respect to the argument and not with respect to x. Differ 
entiating a function of a function we have 

% = + */8x*"-*CV, 


(4) 


& { xin °n(M} = 


W rSiflf haVe becn deriVed by Change of variable from pre- 

~ [x"C n (x)] = x*C n _ l (z). 

half tho V T d 1 n differential ^efficient, the degree of the factor * is 

SI if a X at tLe 

C‘ ( ) ItB Vel0C,ty ls and its acceleration 

taST f0rm ’ WhiCh wm *«"" in the applica- 


(6) 


Tx = ~hPx-«"+»c n+1 (p x t). 


(0 160 ) 


4 
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EXERCISES 


1. Show that the stationary values of C 0 are located by the zeros of C,; the 
stationary values of C n are located by the roots of G'„ + , = C n _,. 

2. Prove (i) C 2 - C 0 = 2 C 0 m ; (ii) C 2 = CJ - C 0 '/x; (iii) C s + 3 C 0 ' + 4C 0 " 
= 0; (iv) x~C n ’ — i + ( n ( n “ 1) — 

3. By repeated application of the recurrence formula;, prove 

(i) £zC„_i = n ^n — ( n + 2)C n+2 + (n + 4)C n+4 — . . . . 

(ii) \xC’ n = \nC n - (n + 2)C n+J + (» + 4)C n+4 -- 

4. Prove that 4 C n " = C n _ z - 2 C n + C n+i . Hence prove by induction that 

2r <L r ° n=Cn ~ r ~ rC "- r+:+r(r 2! ~ Cn ~ r+i -••• + <“ 

where r is a positive integer and the coefficients are binomial. 

5. Prove that 

£ [x"C„(ax)] = a^C n . x (ax), 

~ [x-"C n (ax)] = —ax~ n C n+ i(ax), 

^ [X- *"(?„ V (ax)] = - ia»x-R n + , )C n+l V (ax). 

6. Prove that 

Deduce that 

* Sr ^ = 

<- 2 > r m Ax ~" c ^ 

where r is a positive integer. 

2n , 

7. U <? n (x) = £„*(*) prove - Q n +i = Vn • 

8 By methods similar to those used in the text, investigate the properties 
of a function J n (x) which satisfies the two relations 

2 n 

I n -1 + Ai +1 = 5 ^ n " 1 ~~ ^ n+1 ~~ x n ‘ 
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9. Justify the following statements when n is positive ; 

(i) At any positive zero of C n _„ the functions C„ and C n+1 have the same 
sign. 

(ii) At consecutive positive zeros of C„_ t the function C n changes sign. 

(m) Between consecutive positive zeros of C„_ z the function CL*, has an 
odd number of zeros. 

Interchanging U„_, and C„ +1 in the argument, deduce that their zeros inter- 
lace. I his is an extension of the theorem for consecutive orders. 

10. If c„, K„ are two cylinder functions of the same order (solutions of the 
same equation), prove from the recurrence formuhe that 

CnK n — C n K n = C n+l K n — C n K n+1 . 

Deduce from Bessel's equation that the common value is A/x, where A is an 
absolute constant. By considering two consecutivo positive zeros of C deduce 

“ nd '* °' — 

m^SSS**** t ™'°™ Bcs8e, ' s K|u " lk,n “ nd «» 


(1) tT + ;/ + 4 


(*-*>- 


o. y = <?„(**). 


( 2 ) 


** + ^ + *0 ~ V)y = o- y=c n (Vx). 


(3) xiy’ + y = 0. 

(4) y* + xy = 0. 

(5) xy' + y^o. 

(6) y" + - + 1 t/ -f y = o. 


y = ad<7,($rl). 
y = xiC t (2Vz). 

y = ^"C n (x). 


iicu,ar ' -*■ --uh." r„ b : to i n p ar * 




■he 5 0f C " ° < « < P. prove that 

~^Z:^ a 7o^ ^ «»* *- -ae a teto between 

I8 ' ^ ^ " - «V — have a seated poe^ 
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zero, and that C„ and C n ' cannot have a positive zero in common. Hence from 
the recurrence relation deduce that C n and C n+1 cannot have a positive zero in 
common. 

17. Prove from Bessel's equation by considerations of sign that C n ' has only 
one zero between two consecutive zeros of C n . Deduce further that axC n ' + bC n 
has an infinity of positive zeros. 

18. Ii the coefficients a, b are real, prove that axC n ' -f- bC n cannot have a 
repented positive zero greater than n. [Use the function, and its derivative, in 
conjunction with Bessel s equation-] 



CHAPTER IV 


Bessel’s Equation 


4-1. The series solution. 

If we write Bessel s equation of order n in the form 




with the leading coefficient unity, it appears that the only point where 
the other coefficients can become infinite (known as the singularities) 
>s the ongm To the practised eye it is then evident from the theory 
of differentia equations (into which we do not propose to enter) that 
there is a solution in the form of a series of ascending powers of x 
convergent for all values of x except possibly at the origin itself. 

I he formal method of obtaining the series is to assume 

V = a,*' + a x x r+1 + a 2 x r+2 + . . . . 

::V tS uT at r 3 - be ‘ Uhslil,,lcd in «.e differential equation 

newer Of ° , e , t n . tiCaIljr 2 " 0 ' so that «“ ^efficient of every 

p wer of * vanishes. Tins procedure supplies a set of recurrence for 

mu a, whereby all the coefficients can be determined i”s of the 
modified*form ^ ” Shal ‘ take •¥«*» in the 

*y" + (2n + 1)/ + xy = o, 
whose solution is y = x~*C n {x). We now have 

xy” = r(r - 1 ) a ^i + (f + j )m ^ + _ 

(2n -f- l)y' = ( 2n + ljra.x’-i + . . . , 

Xy = ’ 

«o * r+1 + . . . . 


On taking the sum, the coefficient of *r-i must b 

a °‘ be “° (thC Seriea s „„ 1 eudu.rcV::rave S ri a nd- 
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equation r(r -f 2 n) = 0, whence r = 0, —2 n. For preference we shall 
continue the work using the former value. 

We now have 

y = a 0 -F OjX + a. a 2 + ... , 

whence 

xy” = 2 a%t + 3.2a.jX 2 4- 4.3u 4 z 2 4~ • • • » 

xij = 4- Ui-r 2 4~ 4- • • • » 

(2/j 4- !)</' = ('2n + l)«i + 2( )a 2 x 4- 3( )a 3 x 2 4- 4( )a 4 x 3 4- • . • , 


where ( ) is written for (2n 4~ !)• This on addition gives a x = 0 and 
it appears in succession that every coefficient with odd suffix is zero. 
Our series is therefore an even function of the form 

y = a 0 4- a.a 2 4- « 4 x 4 4- . . . . 


Discarding dead matter we now have 

xy" = 2.1a.>x 4- 4.3a 4 ^ 4- 4- ■ 

xy = a^x 4- a** 3 4- a 4 xS + • 

(2n 4- \)y' = 2(2n 4- l)a^ 4- 4(2n 4- 1 )a 4 x* 4- 6(2n 4- IK* 5 + • 


• • > 


Hence 


2(2 n + 2 )a 2 = —a Qy 


a 


a., = — 


o 


2(2 it + 2)’ 


4(2n + 4)a, - -a 2 , a, - + n + 2 ° )(2 „ + 4) > 


6(2 n 4- 6)a 6 = —a 4 . 


a, 


«8 = “ 


2.4.6(2n 4- 2)(2n 4- 4)(2n 4- 6)’ 


and so on. The mode of formation is self-evident and the coefficients 
are uniquely determined so long as n is not a negative integer. \ e 
have 


y = = ""t 1 ~ 2(2« + 2) + 


2 n 4- 2) 2.4(2n 4- 2)(2n 4- **) 

The factor a 0 is arbitrary and it is customary to put 

1 


- • • • 

4) J 


a„ = 


0 2'T(n 4- 1)' 
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On transposing a factor we have 


( 1 ) J*(X) = 


X n 

2 n T(n -f-1) 


(l-___f-_ - _ 

l 2(2 n + 2) 2.4(2w + 2)(2n + 4) 


An alternative form is 



(2) J n (x) = - ( ^ )n (l _ _ . (k) 4 _ \ 

T(n + l)l 1(» + 1) 1 .2(n+ l)(n + 2) 

This J„(x) is known as the Bessel Function of the First Kind, of order 
n In the case of integral orders the T(n -f 1) can he replaced by 
If e is a zero, so that J n (c) = 0, then J n (~c) is obviously zero and no 
particular interest attaches to negative zeros. We proceed to show 
that J n (x ) comes in the category of cylinder functions. 

The general or (r -f- l)th term can be written 

( — ) r (i*)" +2r __ 

f(r+l)r(n + r-f- f)’ r ~ 0 > 1 > 2 > - 


Hence the general term of nj n (x) -f xj n '(x) is 

This is the (r + l)th term of xj n _ x ( x ) y so that 
(3) nJ n + xj n ' = xJ n _ lt 

corresponding to the recurrence formula 3-1 HI i » 

b“d%^, re ^ ProVed fM fUnCti °“ s 


(4) 

(5) 

(6) 


~ */„' = Xj n+1 , 

J n -! ~ J n+1 = 2 J n \ 

J n-i + J„ +l = — «/„, 

a; 


(7) £ {lV " w} = l 

<8) 2 ' d&y 

(9) 
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and so on. In particular, if n = i we have 

r J ( x 2 jfA ) /2 \* • 

< 10 > = 2-InU) l 1 -273 + 2X375 • • •! U 


Since Bessel's equation is unaltered if — n replaces n, we conclude 
that J_ n (x) is equally a solution. The factor outside the bracket 
is enough to show that it is not a mere numerical multiple of </„(x); 
it is accordingly an independent solution and the general solution is 

y = AJ „(x) -f- BJ_ n {x). 

This breaks down when n is an integer, for J n and J_ n are not then 
linearly independent; in fact we now know, in accordance with cylinder 
functions, that J_„ = ( - )"</„• This result is by no means evident 
from the series obtained above and the explanation is as follows. 
Formally we have 


J- n (z) 

fi (**)* i ( ** )4 _ -•••)• 

F( — n -f- 1) 1 1(— n + 1) 1 . 2 (—n+l)(—n + 2) ) 

Just as every denominator after the second term contains the factor 

( n 2), so every denominator after the rth term contains the 

factor (—n + r); and ultimately, every denominator after the nth 
term contains the factor (- n + n). This tendency for the terms to 
become infinite is nullified by the external factor F(- n + 1) which 
is itself infinite. Accordingly, the first n terms are annihilated and the 
reader can convince himself that the rest of the senes behaves like 
J n {x). The troublesome question of the second independent^solution 
for integral orders will be taken up later in the book. 


4-2. Behaviour at the origin. 

We are now in a position to examine the nature of the function at 
the origin. Evidently from the series we have J fl (0) = 0 for all positive 
orders; and the same is true for negative integral orders since J n 
(_)"J For all negative non-integral orders, J_„(0) tends to infini y, 

the sign being that of T(- n + 1), i.e. negative if the greatest integer 
in n is odd; otherwise positive. This frequently bars the use of he 
function of negative non-integral order in physical problems of the 
type where something has to be finite at the origin. The function of 

zero order is unique; for since 

(\*) x 


( 1 ) 


J 0 {X) = 1 — (ix) 2 + 


( 2 !) 


2 


• > 
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we have «/ 0 (0) = 1. Hence the function J„(x) either passes through 
the origin or goes to infinity, with the single exception of the zero 
order, which crosses the axis at a height unity (fig. 5). 



4-3. The zeros. 


Knowing that J n (x) has an infinity of zeros which, for positive 
orders, start with the origin, we can exemplify some of our previous 
work. Thus it was proved that if the order is less than 4, the interval 
between zeros approaches n from below. The following brief table 
gives the approximate values of the first five roots of J 0 (x) and the 
corresponding differences. 


Root 

Difference 


2-4048 5-5201 8-6537 11-7915 14-9309 

3-1153 3-1337 3-1378 3-1394 


Alternatively if the order is greater than 4 the interval approaches 
from above. The corresponding table for </ 3 (x) is: 


Root 

Difference 


6-3802 9-7610 13-0152 16-2235 19-4094 

3-3809 3-2542 3-2083 3-1860 


77 


We have proved that the first positive zero is greater than n; but we 

showing that the first positive zero steadily 
ncreases with the order. For positive n, m put« = x*J v — x*J 
so that as previously by 3-5(1) ~ X 


t/' + t/i + L 
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It is now legitimate to use the origin as the lower limit of integration. 



Assuming that u vanishes at a whilst v is still positive, consideration 
of sign shows that m is greater than n and the function of higher order 
reaches its first zero later. The tables give the first zero for the corre¬ 
sponding order as: 

Order 0 1 2 3 4 5 

Zero 2-405 3-832 5-136 6-380 7-588 8-771 


4-4. Relation between the solutions. 

If y,, y 2 are independent solutions of y" -f y'f(x) yg{x) = 0, 
we have by 2-5(1) 

ViUz ~ V\'ji = A ex P{~ //(*)<**}. 

In the case of Bessel's equation with f{x) = x~ l the relation becomes 

= -• 

X 

We can determine the constant from the series for the four functions 
involved. The form of the answer being known, only the leading terms 
need be considered. This gives 

. _ 1 — n _ n _1_ 

“ 2 n r(n + 1)' 2~"r(- n + lj 2 n r(n + 1)' 2- n r(— n + 1)* 


Since T(n + 1) = nr(n) this becomes by 1-7(2) 

-2 2 . 

-=-sin?i7r. 

T(n)r(l - n) 77 

Thus 

2 . 

(1) J n J-n - Jn'J-n = -- S1I» tin. 


4-5. The orthogonal property. 

Let u = J n (ax), v = J n (px) so that by 3-6(1) 








BESSEL’S EQUATION 


An integration after eliminating n gives 

(a 2 — p 2 )jxuvdx = x{uv' — u'v). 

The right side certainly vanishes at the lower limit zero, provided 
‘In -f 1 > — 1, a stipulation which ensures the integral bein^ con¬ 
vergent. It will vanish at an upper limit of unity if this makes both 
u and v vanish, or 

— 0 = J n (p). 

This means that a and p are two zeros of J n (x). Replacing them by 
(two of the zeros in question) we have 


f Q ZJ n {CrX)J n (CsX)dx = 0 . 


This may be expressed by saying that the family of functions x*J (ex) 

is orthogonal when the c, are the zeros of J n (x). The property is useful 

in expansions, just as the corresponding property for trigonometrical 
functions is useful in Fourier series 


4-6. Lommel integrals. 

It is convenient to treat here certain integrals involving Bessel 

the", 1°^ Q rr atedWith the nan,e ° f E ' C ' J - von Lommel. Using 
the notation of the previous paragraph, but discarding the assumption 

that a, p are zeros of J n (x), we have P ° U 

(1) (a 2 - P 2 )f o X xJ n (ax)J n (p x ) dx 

= *{P J n(<Hr)J n '(px) - aJ n '(ax)J n (p x )}, 

f*..». 

11 ,i “- 

As a particular case we have 

(2) (a‘ xj n (ax)J n (Px) dx = flJ„{a)J 

^ • —Flex 

conjugate complex p — in w hirh Zer °’ ca 14 “• ond let P denote the 
of the last elation is then ‘t T n J* * zer °' The 
contains the ptoduet of two conjugafe 
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necessarily be positive. The contradiction establishes that J„(x) has 
no complex zero. 

A companion result can be obtained by multiplying the equation 
for u by 2 x 2 u'. Thus 

2 xVm" + 2 xu' 2 + 2a 2 xW — 2 n 2 uu' = 0. 

Since 

i. (x 2 u' 2 ) = 2x 2 u'u" + 2xu'\ 
ax 


integration by parts gives 

x 2 u' 2 + a 2 ru 2 — n 2 u 2 = 2 a 2 fxu 2 dx, 

whence 

(3) fxJ„Hax)dx = 


In parallel with the foregoing we have the case where the argu¬ 
ments are the same but the orders different. Let a = J n (ax), v = 
J m (ax) so that 

+ (" 2 - 1 ) " - °- 
X \ x- / 

From these we derive 


uv 


whence 


x(u" v — uv") -f- ( u'v — uv') -f- (nr — n 2 ) — = 0, 

X 


(m 2 — n 2 ) J U — (lx — x(uv' — u'v). 


or, 


(4) (wt 2 — n 2 ) f - J„(ax)J m (ax) dx = ax{J n (ax)J m '(ax) — J n '(ax)J m (ax)\ 

J x 


Here again the recurrence formulae are capable of producing new forms. 

Another fruitful method of integrating products, also due to Loni- 
mel, may be illustrated by examples. Bearing in mind the formula 
4-1(7), 

i = -***•»> 
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we differentiate the identity 

x r J x J 3 = 


and derive 


— (x r J x J 3 ) — (r + 4)z r — x r (J 2 J 3 -f J V J A ), 


Similarly, using the companion formula on the identity 

x'J 2 J t = x'-*(x*J s )(x*J A ) 


we derive 


d 


di (x'JJ,) = (r - 6 + 4 V, + JV,). 


By addition, 

Ic ( Xr(J i J 3 + *W} = ^{f** + 4)./^ + (r - G)«7 2 J 4 }. 

Since r is at our choice, r = 6 gives 

10 i*?j x j 9 dx = x\J,J 3 + J 2 .7 4 ), 

whilst r = —4 gives 

10 Ix-V^dx = + JVi). 

As an example where the orders are equal, consider 

whence ** " 


Similarly 

whence 


* {^} = (, + 

*^3 a = * r - fl (xV3)* 


~ {x r J 3 2} = (r - 6)x-V 3 2 -f- 2x t J 2 J 3 . 

By addition, 

~ {*W + J 3 2 )} = x-{(r + W + (r _ C)V 2 }. 


Hence 


101X s J 2 *dz = X «(J 2 2 _J_ J 3 2 ^ 

10jz~ S J a 2 dx = — X -4(J2 



54 


APPLIED BESSEL FUNCTIONS 


EXERCISES 

1. Prove that the indicia! equation for x 2 ^ -f + (x 2 — n*)y = 0 is r 2 = h* 
and deduce the series for J n (x). 

2. Verify that tho second root of the indicia! equation in the text, r = — 2n, 
leads to the series for </_ n (x). 

3. Verify that the expansion for J n (x) is convergent for all values of x. Prove 
that it is absolutely convergent. 

4. Prove that the solution of Bessel's equation is AJ n + DJ n f ——. 

5. Prove from the series that J n can have no purely imaginary zero. 

6. Prove that J r '(: r) has no complex zero. Prove further that aJ n + t>zJ n 
cannot have a complex zero. [Assume that a, P are conjugate complex zeros 
and eliminate a t b.) 

7. •/_,(*) = □‘cos*, J 3n (x) = (")*( 81 ^- c 081 )* 

Q‘( — 6inx — and evaluate for the orders 2\, — 2J. 

sinwr: 

8 - + l = 2 — — p 


7 TX 


sin MTT 

•^n^-n-1 + J- n Jn +1 = “” 2 ~ . 

t*JL 


9. If u is a positive integer, 

-1- 2)n V(;) 

ar 1 


<*(**)" 


1°. £ {*«"./„(*»)> = 

{x n, «/ n (fcr°)} = 4ox("+ | )«-V n _ 1 (fcr a ). 
ax 


j x {x-"V n (fcc»)} = -iax<>-")“->J n+ 1 (fcr-). 

11. q -— J j 3*/j - 7 - «)•/$ ~~ • • • • 

/ dx 7T ./ _ 

“T5 = - -F 5 - 

x«/ n z 2 sinnTT «/ n 

r * , * / —n 


2 sin H7T 


•v 


13. [ a xJ 0 '(kx)dx = \a*{J*(lca) + 
Jo 
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/ x :n+i 

x~»+'j n *{x)dx = ^- 2 {•/„*(*) + </ S n+I W}. 

fz-**-'J n ^(x)dx = - W(*) + 4.+.*<*». 

IS. - ^.W^.O)). 

The integral is zero if a, p are zeros of aJ„(/) -f btJ„'(i); and when else? 

16. fxJ n *(az)dz = ***{</„*(<**) - ^.(our^o*)}. 

17. If m -f- w > 0, prove that 

J r i j 

Deduce by a limiting process that 

18. By lorana of the formula for d(x’JJ/d x aud its counterpart, prove that 
J 3*J a (x)dx = 2z i J 3 (x) -f- X*(j* — I6)J 4 (z). 

19. Establish the reduction formula 

f*"' +,J - d * + <»* - »*>/*—■'„«fr - + (m _ 

20. Prove that the solution of xy" + \y" + iy = o is 

xHAJ^Vx) + BJ_ h (Vx)}. 

21. Find the general solution of the equation 

*V - 2 zy' + 4(x* - l)y = o. 

22. If c is a zero of J n , prove 

2 f 0 “ -WM = •V.'M 

Deduce that J n has no zero in common with J J i / , 

the origin. n+1 * "-*» J n or oJ n + bxJ n ', except 

. C, ’. C *’ &C ” denote the zeros of J^ whilst d A * . 

of J t , the origin excepted, prove that the iiTt*rv»i a *’ &c ’ denot o the zeros 
r, and that the interval c , ** r c r steadily increases with 

afforded by the J£?&'£&£££ o,TT * *?“ U 
corresponding fi rat five zeros of J t are ° 8 ' Ven “ the tcxt “ d the 

3*8317 7*0155 10*1734 i^.qoq't 

rerffy e byttl™t toxf " 0 ° 6ui,!r 

function. An application is ^‘w^bhap! ^ 10 the **** 

b prove that 261X58 of '.W when n is greater than 

JJ* *"■ » the Ver^mX 

the normal form of the equation satisfied^! w 7 md< f functions - [Consider 
went of y is an increasing function of x.J Y * * prove that the coeffi- 
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Applications 

51. The lengthening pendulum. 

A problem in mathematical physics usually falls into two well- 
defined parts. In the first place there is some principle which holds 
throughout the body of the medium under discussion. It may be that 
the excess of heat-entry over heat-exit for an element is accounted 
for by a rise of temperature; or it may be that the acceleration is the 
ratio of the effective force to the mass. "Whatever it be, the mathe¬ 
matical expression of such a principle almost invariably takes the 
form of a differential equation. 

In the second place there are special conditions at the outskirts 
and these are known as end-conditions, or more often, boundary 
conditions. A rod may be fixed at one end and so be immune from 
displacement; and a moving fluid, in default of creating vacua, can 
have no velocity normal to the walls of its container; and so on. 

The essential problem then is, so to solve the differential equation 
that it contains arbitrary constants in number sufficient to satisfy the 
boundary conditions uniquely. Theory then usually shows that there 
can be no solution other than this. We propose to investigate a number 
of such problems which depend for their solution on Bessel functions 
of the first kind. We begin with a relatively simple dynamical problem 
that requires but little analysis. 

Problem 1 .—Discuss the small oscillations of a simple pendulum 
when the string is paid out from the support at a constant rate. 

At any time t let the string make angle 9 with the downward vertical 
and let its length be r = a -j- bt. The string is thus paid out at a con¬ 
stant rate r = 6 and the initial length when t = 0 is r — a. The trans¬ 
verse acceleration of the mass is rd -f- 2f$, perpendicular to the string 
in the direction of 6 increasing. This acceleration is retarded by the 
component of the weight, the equation of motion being 

»n(rd + 2 rd) = —mg sin 9. 

66 
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For small oscillations we replace sin0 by 0. and on substituting for r 
and r we have 

(a + bt)6 -f 260 +g9 = 0, 

a reduced linear differential equation of the second order with variable 
coefficents. If the leading coefficient be made unity, there is a single 
value of t for which the other coefficients become infinite. This suggests 
that the equation may be of Bessels type; but as the infinity “curs 

when = - o/6 and not when l = Owe change the origin of time by 
the substitution a + bt = bx, dt = dx. We then have X 


( 1 ) 


<fc 2 xdx x 


_ 0 
6 * 


This is not directly comparable with Bessel’s equation and one can 

the oneTnto thToth"^^? 116 that converts 

?, to * h ® other - The readiest method is to compare with the 
general form 3-6(2). We thus deduce P 

1 “ 2 « = 2, = k -, a 2 = 7j2y2 

X f 9 


a = -h 


Y = h » = ±1, 
P = 2 fc. 


Umiffit^solution*" differe "‘ fr ° m ““<* f ° r ™ can take as » 

0 = Ax~*J j(2A*x*). 

thC ° SCi " ationS are merely 
under which this solution is valid'*1 * TUlI t s ' " e can filul conditions 
With I = 0 we have r = U U- - iA y to the initial conditions, 
angle must have been ’ b ' S ° tUat thc lmt,al value of the 




8 n :z::7 or d r asing ^ - 

^--SSS.-SJSsx 


(o 150) 
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EXERCISES 


1. Verify th.it the angular velocity and the angular acceleration are 


Ak 

x 


J.(2kx\), 



Ak- 

x 


J 3 {2kxl), 


and compare the remark on this form in 3-6(4). 

2. It is dynamically evident that between two consecutive positions of tem¬ 
porary angular rest is a single instant of zero displacement; and conversely. 
Interpret this in terms of the Bessel functions. (There is no temporary absolute 
rest since r = b.) 

3. It is known that the graphs of the three functions J v J t , «/ s cross the z- 
axis, from above or below, in that order. It is further known that the zeros of 
any pair of the functions interlace (see 3-6, Ex. 9). Moreover, from the equation 
of motion, if any one of 0, 6, 6 be zero the other two have opposite sign. Prove 
by the aid of a rough sketch that these statements corroborate one another. 

4. In the case of the ordinary simple pendulum it is known that, (i) whether 
the angular displacement be positive or negative, the angular acceleration is 
always directed to the equilibrium position: and (ii) the displacement and accel¬ 
eration vanish simultaneously. In the present case neither statement holds. What 
is the interpretation in terms of the Bessel functions? 

5. Prove that the angular acceleration changes 6ign only once in each swing, 
and that it happens before the vertical position is reached. 

0. Temporary rest is equivalent to maximum displacement; and zero accel¬ 
eration corresponds to maximum velocity. Interpret these. 

7. Eliminate the time from the differential equation of motion and so find 


a solution for 0 in terms of r. 

8. Prove that the time-interval between two consecutive transits through 
the vertical is 6(c r+1 * - c*)/4g, where c r+1 are consecutive zeros of •/,(/)• 
The interval therefore increases indefinitely. 

9. From the principle of angular momentum, that the rate of change of 
angular momentum equals the moment of the impressed forces, we ha\e 



— mgrO, 


whence 


Deduce that 


xJ z (2kx*) = kj xtJ t (2kxi)dx. 


10. It is interesting to speculate on how what one might call the 9 u ®* er 
periods ” compare with each other. Continuing the notation of Ex. 8, let a r+ , 
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be corresponding consecutive zeros of J.jt). The quarter period for an outward 
swing from the vertical is proportional to rf,* _ c r *. The ensuing quarter period 
for the inward swing to the vertical is proportional to c r+J 2 - d r -, and so on. The 
sixth to tenth zeros of «/,(0 taken from the tables are 

19-016 22-700 25-904 29-074 32-190. 

The sixth to ninth zeros of J 2 (t) are 

2MI7 24-270 27-421 30-569. 


These show that the successive quarter periods, beginnin 
raent, are proportional to 


g with an outward move- 


Out: 00-14 71-02 80-89 90-70. 

In: 72-09 81-96 91-83 101-70 




°r. 


r+l » 


V ~' r * <*,+,•-,1 >>,!„,*- 
+ C '+'‘ “r*f *? + i n f < tc w * 


5-2. Vibrations of a taut string. 

direction^of the' ^“beTT ° f •***'* ^ 

s*.s 

titles at B arc X + SX, Y + hY T 7ST Ti° COrre fP ond j n g 
force is SX; and as it is custom ’ * ' T1,6 r<? sultant horizontal 

to be negligible, we have S* ,^1“ x'JZtsl't'r 7T 
so i?ia“ : 8 t °V mM *° bC SeCD "' ith tbe "ye, 

^ to the usual a M „mp t -o™han.he 0 .r P ' by Tbis 

tbe length At riirhf nn 7 T 7 * *. tension 18 consta nt throughout 

sr = «TrAX b! ® the resu,tant force '•p™«7s 

mass of the element is „ Ss andIsTt7’- "°{ necess " ril y ™»stant, the 
‘ion la y. This gives ^equation of mottn ' “ * iU accd ™- 

3 2 v 3 

^pSs = sr = ?. ( rsm^)s,. 
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As sin <// = dy/ds and the distinction between x and s is not worth 
making, we have 


P 


cj/ 

a / 2 



v-l 

dx 2 


So far we have followed the routine procedure. We now propose to 
modify the problem so as to admit of treatment by Bessel functions. 


Problem 2—Discuss the transverse vibrations of a non-uniform taut 
string whose line density at distance x from one end is p( 1 + kx). 

This may be regarded as a first approximation to a wet string in a 
vertical position. The equation of motion is 


3 - 


d'-y 

d: 2 ' 


On the ground that anv periodic variation of the displacement with 
time can be expressed in a Fourier series we write 

y = u sin (pt + <f>), = — P 2 u sin (pt + <f>). 


This introduces a number p whose different values determine the 
frequencies of the various harmonics. The equation of motion is now 

^+ s r <1 + = °- 


A slight modification is required to achieve a tractable form, we put 

1 lex = kz, dx = dz, 


whence 


<^+p*Pkzu = 0 . 

dz 2 ^ 1 T 


The equation is in the normal form, so if it turns out to be soluble in 
terms of Bessel functions they will contain the factor z*. Comparing 
as before with the general standard we have by 3-6(2) 

a = y = o, n = ± b 
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A solution is therefore 

« = Az*Jtft*») + Bz*j_ i (/3z 3 % 

It remains to determine the constants A and B from the fact that 
both ends are fastened. Taking the origin at the left end of the string 
we have as end-conditions 

1 

* • 
k' 


x = 0 = y =-. u, z 


These give 
and 


x=X, y = 0 = ii, 2 =g +A ). 


0 = AJtfk-**) + BJ_^k- 3n ) t 

Putting for brevity 

= c, (1 + kA) w = n, 
the elimination of A, B gives 

The number n is known from the conditions of the problem The 

solution of a t_ den nl T" ' J aVaiIaWa the 
If we write 

f{x) =f{nx) = J^xyj^x) 

the graph of y =f(x) has the zeros of,/ \\ru._ r , 

a zero, the graph switches from + co to thr0M « h 

graph accordingly looks aomewW u- * 4 conve rsely. The 

2fon 0f Js7a{ ( T ) it foLTtl Tt e h Pt tbat ° 

- - a rough 
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minate; but one of them remains arbitrary. The most general solution 
of the partial differential equation for the displacement is 

y = «i(z) sin(p x < + <f>i) + u t {z) sin (p.£ + <f> 2 ) + • •. . 

Such constants as are still present can theoretically be determined to 
conform with initial conditions of shape and velocity. 

It remains to examine the question of nodes. The displacement y, 
and consequently w, is permanently zero if 

JAP**) _ _ B _ Jtfk-™) _ 

A J-dfik-**) J-W 

For a permissible value of c the fraction on the right is a definite num¬ 
ber. The function on the left, as already explained, ranges from -f- co 
to _ oo and therefore takes all intermediate values repeatedly, so that 
there will be nodes when the vibration exceeds the gravest node in 
frequency. 


EXERCISES 

1. Taking the line density as p(l + lcx)\ prove that the solution depends on 

2 ljj((} 2 =), = lpW(%) 

with the corresponding function of negative order. Deduce the equation which 
determines the possible frequencies. 

2. Investigate the problem when the line density is p(l + **)*» showing 
that the solution depends on functions of order §, — -f. 

3. When the line density is p(l + kx)~l the necessary functions are of order 

2 _ 2 

•j ^ • | 

4. Solve the problem when the density is proportional to & and the ends 
arc defined by x = Xj, X 2 . 

5. Prove that the case where the law is p(l + **)"* is not soluble in terms of 
Bessel functions. Also that the law p(I + hc)~ x requires two functions of the 

first order. . 

6 Prove that a law of density can always be found so that a string of given 
length can vibrate in the shape (apart from scale) of any specified loop or series 
of loops of a function x*J n ((3x), where n is assigned and p remains to be foun . 

5-3. Stability of a vertical wire. 

If a straight wire be clamped vertically at its lower end, the posi 
t.ion is stable when the wire is short, A longer wire may find a more 
stable equilibrium in a curved position. We propose to investigate 

this. 
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Problem Z.—A uniform wire of length A and weight w per unit length 

has its lower end clamped vertical. Discuss the condition of insta¬ 
bility. 

Measuring from the lower end 0, consider two adjacent points 
A, B such that OA = x, AB = 8x. Presuming that A is situated a 
distance y from the vertical through 0, the corresponding distance 
for B is y -f- by. At the onset of instability, when y is small, the bend¬ 
ing moment M at A is given with sufficient accuracy by M = Elu" 
where B is Young’s modulus for the material and I is the moment of 
inertia of the cross-section. The weight of wire above B is not sensibly 
different from that above A; but it acts at a shorter leverage. The 

bendmg moment M + 8M at B is accordingly slightly less than that 
at A and we have 


Put 


then 


bU = w (A - x)8y = - Sx. 

(lx 

f-"S— 


A — X = 2 , dx= —dz, J£_ — 9 /.2 


The equation becomes 


El 4 


dfy 9z& 

+ ~ p = 0. 


^ *^1*™ “tt ther n f ° re ** tUInS ^ -'“Wo in 

"ith the general C “ tam the **■ Comparison 

“ = i. v = i » = ±J, /9 = i. 

The solution is therefore 

It remains to determine the constants. 



64 


APPLIED BESSEL FUNCTIONS 


At the upper end, x = A, z = 0 and as there is no applied couple 
we have 



dp 

dz 


In view of previous experience it is possible to write down (4-6, Ex. 10) 

d -l = UkzJ _ x {kz ,ri ) - | BkzJ i (kz w ). 
dz 


It is left to the reader to verify that the series for the former term 
starts with a constant, whereas that for the second terms starts with 
2 2 . Since dp/dz vanishes at z = 0, the value of the constant must be 
zero, whence A is zero. At the foot of the wire, on account of the 
clamping, we have 

x = 0, z = A, p = 0. 


These give 


0 = BX i J_ i (kX sn ). 


The only rational way of satisfying this equation is to take the 
argument *A 3/a as the first zero of J_j(x). This is known to be approxi¬ 
mately 1-87. Id any given case where w, E and I are known, k is deter¬ 
mined and hence A. It remains to add that if the wire is a ribbon or 
of any cross-section other than circular, the instability occurs in the 
direction of greatest flexibility and the appropriate value of I must 

be chosen. . * 

“ I have no satisfaction in formulas,” said Lord Kelvin, unless 

I feel their arithmetical magnitude- at all events when formulas are 
intended for definite dynamical or physical problems.” For a wire of 
radius r made of material of density p we have 

w = Trr 2 gp, l = 

w _ o 1.2 _ *9P k = — l(-\ 

A reasonable value for E is 2.10 3 tonnes/cm. 2 , or 2g. 10 9 dynes/cm. 2 . 
If the diameter is 1 mm. and the density is 7-6 gr./c.c., this gives 

k = 1-64 X 10- 3 and A = 109 cm. 
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EXERCISES 


1. It states in tho text that the only rational way of satisfying the condition 
of stability is to take the smallest zero of J_ j(x). What is the objection, if any, 
to choosing zeros other than the smallest ? 

2. Adopting tho solution for p given in the text, with the accompanying 
value of dpjdz, verify that the normal form is satisfied. 

3. If the wire carries a top weight W, prove that the greatest length for 
stability depends on the first root of the equation 


J' J (na)./ | (a) -f- *7_ J (na)J_ J (a) = 0; 



4. Since 




we have y= — fp dz to determine the form of the wire, 
variable that this is equivalent to dt. The integral 

finite terms. 


Prove by a change of 
cannot be evaluated in 


5. If tho wire is a solid of revolution, the radius r at depth x below tho top 
mg Xx and tho weight above tho section being px", show that the solution 

depends on J (x) = 0, where p = —~ m ~ 1 , 

n — 4m -f- 2 


6 # 4. Instability of the deep cantilever. 

emJe^ 1 ll We )\lr 0 , Wn i th . a J fc * de<?P cantilever tends to lateral instability, 
espemlly if the load rides high. The reader can easily convince him- 

at the Ind* strip of pa P er “"W a paper-clip 

the case of a 1 U ? P<?r l dge - We pr °P° se to investigate this, taking 
the case of end-load on the centre line. b 

Pr0i W i -~ DiSCUSS the of ,he deep cantilever under end- 

„ r d° r IZ'T'lTr ta !; C ‘ v 6 * 8XiS *° thc ri 8 ht > V vertically up- 

(4 6) r dcr ' and «»** acti ”8 «P™rds 

Le bLmt a t f rfl b *“\ b . e “ g * -d end-load IF, the centre- 

point B on it W T- CUIVC Wbch Was ori g inalI y ‘he * axis. Any 
accuracy X v z °^ in * tea *» y ] 2 and lts end-point is with sufficient 

ttrongh anck’s w° rOS r eCUOn at B is P resumed to have turned 
8 angle 0 clockwise when viewed from the origin. 
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There are two component couples at B which we may denote bv 
C x , C t . Their magnitudes and direction cosines are 

C x =W(z-z 0 ); 1, 0, 0. 

C t = JF(A — x)\ 0,0,1. 


It is presumed that a plane section remains plane and normal to the 
centre line; the normal to the section at B thus has direction cosines 

dx dy dz 
~ds ds ds 


where ds is a line-element of the centre line. 




As the twist per unit length is proportional to the normal torque, 

WChaVe C dz + C dx -A dl> 

C -d- s + Cx ds~ A ds ' 

where A is a constant of proportionality depending on the cross-section 
and the material. We can differentiate this with respect to s, whence 


T C 'T + C 

ds ds 


dx\ _ d°-6 

ds) ds 2 ' 


On the grounds that x and s are scarcely distinguishable we put 

dC t 




ds 


= -w. 


dC x _ dz 

ds ds 


= A % = ,r < A - ■» 5? 


whence 
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A line which was originally vertical on the section at B has direc¬ 
tion cosines 0, cos0, sin0; or, 0, 1, 6 approximately, on the grounds 
that B is small. If the couple in this direction be offset against the 
corresponding bending we have 

C,e = JT(A - 1)9 = -El % 

as 2 

Eliminating z from these two equations we reach 

eia + ir 2 (A - xy-e = o. 

This is soluble in terms of Bessel functions if we make some minor 
changes. Put 

X~x — t, dx = —dt, -ill = k- t 


whence 


El A 


d*0 
dt 2 


+ = 0 , 


by 3 h C(2) m the n ° rmaI f0rm * Com P arison our standard gives 

a== i Y= 2, /3 = \k, n = J, 

The solution is therefore 

The arbitrary constants P, 0 remain to be determined. At the free 
ena x A, J = 0 we have no applied couple and 

de = __ 

dx U 

On referring again to 4-6 Ex. 10 we ean write down 

de 

£ = - QU' n J,(Ut*). 

“7^w^e“ C det: S t P istrl 0rmer tem ’ ^ "‘ h “ 

At the clamped end we have 

* = 0> ( = A ' 9 = 0 = <?A‘J_,(UA 2 ). 

«b'es give the fat zero of /_,(*) an 2 0063. so that approximately 

A = 2i-t, or 2 !y/k. 
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An equivalent to the foregoing analysis is usually accepted in texts 
on structures or materials; in practice, the tendency to instability is 
obviated by modifying the section of the beam. An alternative dis¬ 
cussion based on a different set of ideas will be found in Temple and 
Bicklev, Rayleigh's Principle (Oxford). 


EXERCISES 

1. Prove that if the load is uniformly distributed the solution depends on 
functions of order 1 /S, — 1 /6. The smallest zero is given by J n (x) = 0, n = —1/6, 
a- - 2 1423. 

2. If the load varies uniformly from zero at the free end, the solution depends 
on functions of order ±1/8. The smallest zero is given by J„(x) = 0, n = —1/8, 
x — 2-209. 

3. Discuss the stability of the centrally loaded deep beam when both ends 
arc encastree. Prove that it is not fundamentally different from the case dis- 
cussed in the text. 


5-5. Critical load for a variable strut. 

One of the simplest applications of the theory of linear differential 
equations with constant coefficients is to the theory of struts. It is 
presumed that the reader has some acquaintance with this; it can be 
found in numerous texts on structures or materials. 

The problem in its simplest form is best visualized as the equili¬ 
brium position of a flat steel ribbon whose ends are joined by an in- 
extensible string of natural length slightly less than that of the ribbon. 
Taking the line of the string as x axis and one end as origin, let y measure 
the departure from the straight at any point. The couples gi'e as t ie 
approximate equation of equilibrium Ely" Py = 0, where is 
the tension in the string, E is Youngs modulus and 1 is the moment 
of inertia of the section. If these be taken as constants, a suitable 
solution is y = A sinnx where A is arbitrary and n 2 = P/EI. 

The deflection is certainly zero at the end where x is zero. Ta mg 
A as the length of the string, the absence of deflection at the other en 
implies sinnA = 0, whence nX = tt, 2tt, 3tt, . . . . A chosen one o 
these determines n and thence P. The inference is that a possi e 
equilibrium position is any number of loops of a sine curve. A1 t ese 
are very unstable except the single loop. 
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Problem 5 .—Discuss the stability of a triangular strut. 

As a first modification of the analysis sketched above, wo assume 
the ribbon cut to a triangular shape. The moment of inertia / is then 
proportional to the width at any point, and this in turn is propor¬ 
tional to the distance from the vertex. Taking the length of the ribbon 
as x-axis and the vertex as the origin, we can replace I by II.r. The 
equation of equilibrium is then presented to us in the normal form 

EHxy" + Py = 0, 



<Py 
(lx 2 


l- 

+ 4** = °’ 



The solution must contain the factor x*, and a comparison with our 
general standard 3-6(2) gives 

a = h = y, 7i=l, ft = k. 

An acceptable solution is therefore 

y = Ax*J x {kx*). 

This certainly gives zero deflection at the left where x is zero. For 
zero deflection at the right where x = A we have = 0, so that 

• is a zero of ,/,(x). Any particular zero serves to determine k, and 
n “ “ tUr “ detern | ,n es P. The smallest zero corresponds to the stable 
p sition of a single loop. A configuration in several loops is theo- 

uXiroof Z WOM bC highIy UnStaWe; H coms P°" ds to a 

maxtoum df 1^°” in ? singIe loo P there is 11 sin g le point of 
maximum displacement, given by 

^ = 0 = iAkJ 0 (kx*). 

n^t™ :t: heorem that betwecn *•» »» 

fon^'w™ 8 ^, t0 flgUrCS ’ ,et 25 om ' be the length, 1 mm tl.o .ini¬ 
tiated attnTt d' T'‘ hC “* br ° ad e,,d ' This gives // cal- 

«»d the tIZo „3 83 ^ *“ S 3 ’ I0 ° W™' 2 

Attentinn ! x ? 3 ’ 8317 we have k = °- 766 3 and P = 2-45 kg. 
The point of mlximum^ ^ 1^° j • • - 

“fc^v . 8 the 

the more fltbletd “ definitel y dis P' aced from the middle to-.vardl 
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EXERCISES 

1. If the ribbon is blunt-nosed, with / = Hx d, prove that the solution i6 

y = Axh/j(fcrJ), k=*(l L)\ 

2. If the ribbon is sharp-nosed, with / = Hx* ’ J , prove that the solution is 

y = AziJ 2 (kxl), k = 

Deduce that between two consecutive zeros of J 2 is a zero of J x . 

3. Prove that the case / = Hx 2 is not soluble in terms of Bessel functions. 

4. In any case where the ribbon stands in a single loop, what is the physical 
justification for saying there are no inflections? Explain this by the equilibrium 
equation and interpret it in terms of the Bessel functions in any particular case. 



5-6. Railway transition curves, Fresnel’s integrals. 

It is distinctly unusual to find Bessel functions associated with 
geometry ; but the following problem is essentially geometrical. When 
a locomotive rounds a curve there is a definite side-thrust on the track, 
and the reaction supplies the necessary centrifugal force. In practice 
the matter is minimized by giving a super-elevation to the outer rail, 
thus canting the engine inwards. So long as the running is on a straight 
track there is no centrifugal force. It follows that if the straight track 
joins on to a curve of finite radius the force jumps from zero to some 
definite value as the engine passes the join. This is practically an 
impulse and is detrimental both to the permanent way and to the 
comfort of passengers. Moreover, it means reducing speed and losing 
time for the sake of taking things easy. The problem therefore is so 
to design the transition curve as to avoid the necessity for reducing 
speed and to ensure that the centrifugal force rises uniformly from 
zero. 

So long as the engine runs at a constant speed it is immaterial 
whether the uniformity is with regard to time or space. If P be any 
point on the curve, let p be the radius of curvature and s the arcual 
distance, measured from some point to be specified later. For a con¬ 
stant speed the centrifugal force is proportional to p~\ and if it is to 
change at a constant rate we have 

d -Bll = 26 2 . 

ds 
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The constant is given the form 2 b 2 as a matter of convenience. Hence 

P~ l = '2b-s + c, 

and the constant c can be eliminated by measuring the arc from the 
join with the straight track. Since p — ds/dtfi we have 

d 4- — 2b 2 s, 6 = b-s-. 

as 


No integration constant is required. Taking the origin at the join, and 
the x-axis as the straight track continued, we have 


Hence 


^ = cos ip = cos (6s) 2 , ^ = sin^ = sin (6s) 2 . 

as as 


x = f cos (bs) 2 ds, y — f s\n(bs) 2 ds. 

Jo J o 


The change of variable bs = t gives 


bx = f cos t 2 dt, by = f sin l 2 dt. 

JO Jo 


These integrals make their appearance in a subject about as far removed 
from the crudities of rolling stock as well could be; they are known as 
Fresnel’s integrals and they occur in the diffraction of light. If we put 

l z = u, dt = du, 

we have 


Similarly 


2 0' bx =f°Qu) >cosud '‘ =L j - 
<!) h = foiluj Anudu »<«>*«• 


They can be evaluated by the repeated use of the recurrence relation 

w , . = «^n-l “ «/„+!• 

We have m succession 

=j-, - 

= j 3/3 — «/ 7/3 , 

2 J*n = — Jnn * 
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and so on. 
n -> oo , 


so that 


Hence by summation, since J n (x) obviously 

\J- j = + • • • • 



6x — + J s/a + + • 


0 when 


The series is convergent and numerical results can be obtained by the 
use of the tables. The companion result is easily shown to be 



ty = J 


11/2 



• • 


They are both apparently due to Lommel. It need hardly be added 
that no civil engineer could be induced to utilize such results in laying 
out the curve. 



CHAPTER VI 


The Second Solution; Further Applications 

61. The indicial equation. 

If the series y = a 0 x T + a 1 * r+I + ... be substituted in Bessel’s 
equation 

x-y" + xy' + ( x 2 — « 2 )y = 0, 

it follows as in 4-1 that the resulting indicial equation is r- = n-. The 
verification is quite simple; for the coefficient of o 0 on substitution is 
r(r — 1) -f r — n 2 . This agrees with the result stated and the roots of 
the indicial equation are r = n, — n. In general the two solutions J n 
and «/_„ are independent since the one is obviously not a mere num¬ 
erical multiple of the other. Exceptionally, when n is integral, J_ n is 
merely J n except possibly as regards sign. We now have to investigate 
the second independent solution for integral orders. 

Certain general propositions concerning the indicial equation will 
be mentioned here. They are; 

*'1 If t]ie roote of the indicial equation are distinct and do not 
diiter by an integer, each root leads to an independent solution. 

(n) If the roots differ by an integer, the second solution may, or 
ma y not, contain a logarithmic term. 

(iil) In the case of repeated roots, the second solution certainly 
contains a logarithmic term. ^ 

n rn T f heSe .i a i, re P “ rt ° f fte the0ry of diff «™tial equations and no 

wldch VVhe tr bey °? d What appliCS to Bessel ’ s equation, 

r?™," . a ! m0 /‘ pcrfect exc “P lar - They are already substantiated 

f tCL * t r \ r00tS 0f the “ dicial e d uati on differ by 2», and 

ar IL"? hST dd neither ^ “ Whi ° h devolutions 

may or mav not , T*? 1 - u*" “ 8 eneral the elution 

p * fl ? J fc ’ contam a logarithmic term. In the case of Bessel’s 

be ° eS n ° fc; f ° r the funcfcions of °rder half an odd integer can 

has a repeated root and th * \ 1S Z<?ro m dicial equation 
(o ifo? 0t and the0r y ***** tba ^ the second solution must 
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contain a logarithmic term. The suggestion made in 3-6 (for the case 
y = 0) was not devoid of purpose; it can now be disclosed that the 
solution in question is x“ logx. 


6-2. Equations of zero order. 

We now proceed to verify the statement about zero order for 
Bessel’s equation; and we begin with a method that makes the result 
look plausible. From 


x- 


W — 1 -T + • • • 

4 


we have bv reversion of series 

1 


x 2 


Hence 


J„(x)- 1+ 4+ • ' • 

1 = ! + ? + 

—- I I • • • 


:{J 0 (*)} 


x 


From 4-6, Ex. 4, and 2-5(2) it is known that the second solution can be 
expressed as 

y - 

whence we deduce that 


y = H*) { lo g x + ^ + * ' • } ’ 


the first term being followed by a series of ascending powers of x. 

It would be tedious by this method to find the actual form of the 
scries, which is in fact of no great use to us. The important point is 
that the value of y is infinite when x is zero, from which it follows that 
the solution is linearly independent of «7 0 . It is customary to denote it 
by Y 0 . The presence of logx confines us to positive values of x, but 

in practice this is no handicap. 

Having given our statement an air of verisimilitude we turn to a 
more rigorous method, applying it in the first case to the equation of 

zero order 

y" + -y‘ + » = o. 

X 


If we put 

L x 2 _ x* _ 

C y (x) - x^i + 2 ) + 2?4(27+lj(27+T) 
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so that incidentally C 0 (x) = J 0 (x), we know from 4-1(1) that C„ satisfies 
Bessel’s equation 


Differentiating tin's with respect to v we have 

(P id 
dx- 




On equating v to zero, the right side vanishes. We conclude that 




is a solution of Bessel's equation of order zero. It remains to find the 
explicit form of this. Write C v = x'S where S stands for the series 
in brackets. We then have 


£ c, = *« log* + «•• as 

vv 0V 

= C„ log* -f x* 

cv 


On equating v to zero we have 


c ^ x) } „ = ,o e* + {=r) 

\ov /*• = o [dv) ya o 


It is left to the curious to verify by logarithmic differentiation that 

+ + .... 

The result 

T ° {X) = logs: + {<J*)> - ^ (1 + i) + . . . j 

is usuaHy associated with the name of Neumann. The full solution of 
the equation of order zero is 

y = AJ 0 (x) + BY 0 (x) 

and no other solution can have any form but this. We may here insert 
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a remark, of no practical consequence, that in a particular form asso¬ 
ciated with the name of Weber the values 

A = - 1 (log2 - y); B = -? 

7t ■fl¬ 

are taken. The number y is known as Euler's constant: its value is 
0-5772 . . . and it is defined as the limit of 

1 + ^ + 3+ • • • + ~ — 1°S n 

n 

as n approaches infinity. 


6-3. Equation of integral order. 

Having substantiated that the second solution of zero order con¬ 
tains a logarithmic term, we now proceed to render probable the same 
for integral orders. We can write 

J n (x) = 0x n (l -f- a t x z -f n 2 x* -f . . .), 


where fi. a v a 2 , &c., are known constants. We then have an expansion 
of the form 


Hence 


J n (x) 

1 


= -f- b^x- -f- b.gfi +...). 


- (1 -f c x & -f- C 2 X* +...). 


ar{Jn (*)} 2 fPx 2n+l 


The (n + l)th term is a multiple of x _1 and hence on integration we 
get a logarithmic term. If we use the series value of J„(x) in the second 
solution. 



we conclude that the form is AJJx) logs -f- -f- S 2 , where 


A = constant, 

= a finite series of negative powers, 

S 2 = an infinite series of positive powers. 

We do not propose to find the forms of any of these since they would 
serve us no useful purpose. The result is denoted by I „{x) an 1 
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becomes infinite when x is zero (fig. 7). The full solution of Bessel s 
equation of integral order is 

y = AJ n {x) + BY n (x), 

and every solution must be of this type. 

The weakness of the above demonstration lies in our failure to 
establish the presence of the multiple of x~ l . However, as the series <S X 
exists, we are still justified in saying that Y n (x) tends to infinity as x 
tends to zero. Instead of pursuing the matter it is more profitable to 
adopt a new outlook. 



v. 


Let </„(x) be a solution of Bessel’s equation of non-integral order 
By differentiation with respect to v we have 


dx* \du y \ 


+ 


+ ( 1- s) 




Since J- y (x) is also a solution we have similarly 


£_ (!; I , 1 d 

dx*\d v ~ y l ^~xdx 




Multiply the second of these by (-1)", where n is an integer, and sub¬ 
tract In order to avoid piling up formidable masses of symbolism use 
the abbreviation * 
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We then have 
d-F.. 


^ , 1 dj\ 

dx 2 x dx 




We notice that if the non-integral v tends to the integral value n, the 
right side tends to zero. Hence if F n denotes what F v becomes in 
similar circumstances, we conclude that F n is a solution of Bessel’s 
equation, i.e. 

Fn = Lt \^~ v ~ ( —1 )" ^)}> v ~+ n - 


The reason for the foregoing rather abstruse discussion is partly 
the matter of tabulation and partly dissatisfaction with having the 
second solution differently defined according as the order is integral or 
not. The logical thing to do is to adopt a form of the second solution 
that holds in all circumstances. Various writers have done this, and 
the corresponding functions usually go by their names. The only one 
that need concern us is Weber’s, 




cos vtt J v (x) — J ~ v ( x ) 
sin vtt 


For non-integral orders this is evidently a solution since it is a linear 
combination of «/„ and J_„. For integral orders, coshtt = (—1)", 
sin mr = 0, so that the fraction takes the indeterminate form 0/0. 
Resolving this in the usual way by differentiating numerator and 
denominator with respect to v, we have 

9 T 

COS VTT — J v 
OV 

TT COS I TT 




The result shows that as v approaches n the fraction becomes F„/tt 
and is accordingly a solution of Bessel’s equation. 

Much of the foregoing can be laid aside. All that matters in prac¬ 
tice is that there is a second solution Y n (x) that becomes infinite at 
the origin. It may be tabulated under the name of Weber or Neumann 
just as logarithms may be tabulated under the name of Briggs or 
Napier. That is immaterial so long as they are not confounded in the 
same calculation. When we meet a Bessel equation of order n, be the 
order integral or not, we write down the solution as y = AJ n -j- BY n , 
and go on from there. 
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6-4. Recurrence formulae. 

A consequence of Weber’s definition is that Y„ obeys the same 
recurrence formula as J n , which is a great advantage. If we multiply 
the formula 

-y {*■/„} = 

ax ax 

respectively by cot vn, cosec vtt and subtract, we have 

(i) ± {^y.} = 


This establishes that Y„ is a cylinder function and the rest follows 


(2) 

Yy- 1 

+ y v+1 

II 

»l? 

* 

w 

(3) 

5',-i 

~Y v+l 

= 2Y ' 

1-1 

(4) 

*y; 

+ vY y 

ii 

* 

1 

(5) 

*y; 

~vYy 

= -*Y« i, 

and so on. 

The proof is based on the 

assumption 


- - - - - — — X-J- — »•* V- J •• v 1 I v l^OUU'O 

to integral orders. 

Since Y„ is a cylinder function, we conclude that it has an infinity 
of zeros, their interval tending to tt\ and as x*«/„ and x*Y„ are two 
solutions of the same normal form, we see from 2-9 that their zeros 
must interlace. Moreover, there must be a relation of the form 

J.YJ-J r 'Y. = ± 

X 

(see 3*6, Ex. 10). Using Weber’s definition of Y v and its differential 
coefficient, we have 

(6) J V Y, JJY„ = (J V 'J_ V — </„./__/) cosec vtt 

_ 2 

* y 
7TX 

the last step being from 4-4(1). It must be pointed out that the right 

side may take a different value if some second solution other than 

Weber s is used. This has to be borne in mind when consulting other 
books. ® 
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EXERCISES 

1. Establish the following results, which frequently occur in physical prob 
lems: 

(I) J n yn + l-Jn + lY n =- ^ 

I >x 


(2) J n Y n '-J n ’Y n = 




2 / n 2 \ 

(3) J n ‘ Y n ’ - J n " Y„' =- x (l-^> 

2 /3n* \ 

2 * ^{** ny n(fc**)}= ifccKn-Oy,,.^), 

^ {x-i"y n (ixi)} = -§fcr-H"«)y n+I (te»). 

3. ^ (x^y^i^)} = i«*("+‘).-iy n _ I (fa-), 

^ {r-"-y n (fcr')} = -IxrO-")-«y n+1 (Ix-). 

4. y 1 (x)= -J_j(x)= -□‘cosx, 

y_l(*) = •/*(*) = (^ysinx. 


6-5. The lengthening pendulum (resumed). 

Now that we are in possession of a second solution, we can con¬ 
tinue with the applications. We revert to the problem of the length¬ 
ening pendulum. The equation of motion, 5-1(1), 


d 2 d 

dx 2 



has the solution 
whence 

dO 

Jt 


6 = Axr*J 1 (2Axr l ) + Bx* F 1 (2Avr»), 

= d6 = _Ak j, 2kxi) _ Y 2 (2kx'). 

dx x x 


With two constants at our disposal we can satisfy two initial condi¬ 
tions. Suppose we take the origin of time when the angular velocity is 
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temporarily zero. Let the corresponding angular displacement be 
/?. We then have 

t = 0, r = a, 6 = /?, 

dt ’ 6’ \b) ' 

For the determination of A, B \vc have 

/*(£)*= M X W + B l'.(A), A = ? („,)*, 

0 = .<4./ 2 (A) -f- ,(A). 

In virtue of the relation 6-4, Ex. 1, 


«/i(A)J r 2(A)-y 2 (A)r i (A) = - 


we deduce 


77 A’ 


A - -I^(")V 2 (A). B = JttA^fA). 


Hence the angular displacement at any subsequent time is given by 

2 6 (bx \» r „ „ 


EXERCISES 

equation° Ve ^ ** ° f Z ° r ° dis P ,acemi?nt given by the roots of the 

■W _ */,(2fei) 

Y*W r,(2fcr»)’ 

and that the instants of temporary rest arc given by the roots of the equation 

JiW _ Jj(2fcr*) 

W F.(2/b:»)‘ 

What is the significance of the obvious root, X = 2Lr» ? 

deduce the recurrence formula. 1 equation of motion and 

velocity^ 0 ri<i “' “ k “°™ to *» calculate the traverse 
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5. In what circumstances would the solution depend on the function Y with¬ 
out J ? 

6. Prove that the angular displacement and the lineal lateral displacement 
do not reach maxima simultaneously. 

7. Assign two reasons why the successive maximum angular displacements 
decrease. [The case of lineal displacements has to be left open for the moment; 
it depends on the ultimate behaviour of the cylinder function 

8. Prove that the constant A has the physical dimensions Tb 


6-6. Motion of a variable mass. 

Problem 5 .—The rectilineal motion of a variable mass under a variable 
force. 

Suppose a mass moves along OX under an attraction to the origin. 
If the force per unit mass is proportional to the distance and we equate 
it to the rate of change of linear momentum, we have 


i L d *) = - 

dt \ dt / 


c 2 mx, 


where c 2 is the constant of proportionality. Hence 

dt 2 m dt dt 


Suppose further that the mass suffers abrasion or any form of attenua¬ 
tion so that its magnitude at time l is m = (a + bt)~ l . Thus the mass, 
initially a -1 , is asymptotic to zero. Then by logarithmic differentia¬ 


tion, 


1 dm __ b _ 

m dt a -f- bt 


and the equation of motion becomes 


d 2 x 

dt 2 


b dx 
a bt dt 


-f crx = 0. 


Make the substitution 

a + bt = bz, dt = dz, 


and we have 


d 2 x 1 dx 



THE SECOND SOLUTION; FURTHER APPLICATIONS 


S3 


Comparison with our general form 3-6(2) gives 

1 — 2a=—1, ny = a, 2(y — 1) = 0, f$y — c, 


whence 

The solution is 


a = 1 = y = n. /? = c. 

x = AzJ y (cz) -f- BzY^cz), 

and the motion is oscillatory. The velocity at any instant is given by 

^ ^ = AczJ 0 (cz) + BczY 0 (cz). 

To determine the two constants from assigned conditions, let us sup¬ 
pose that when t is zero the mass is at temporary rest at distance h 
from the origin. This gives z 0 = ajb and 

AJ„(X) + Sy o (A) = 0, A = ~ 

b 

AJ X (\) + £y,(A) = h -. 

a 

In virtue of the relation 6-4, Ex. 1, 

- J 0 (X)YM = A = , 

7tA 77 ca 


= 0 gives 


these solve to 

A = Wy 0 (A), B = -hnchJ 0 ( A). 

Suppose we require the time to reach the origin; then z 

J i(&) _ _ B _ J 0 (A) 

^i(«) A Y 0 ( A)' 

The fraction on the right is then a definite number, $ say. Hence 

J^cz) — sY^cz) = 0 . 

Here the left side is a cylinder function and there is therefore an in- 
flnite number of suitable values of cr. Their interval is greater than 

tra^h ^ *° The time interval successive 

transits through the origin tends to the value „/c, as in simple bar 

momc motion. In view of the fact that the foL per uXmass7s 
ment 0r h 10n d t0 a e distance . * would be interesting to see an armi- 

t“e rnofa is° no^ Pr ? Ci , pleS ' e ? laini "g the apparent anomaly that 
motion ,s not simple harmonic. The total mass varies; but the 
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total force varies in the same ratio. Having disposed of that trifle, 
follow it up by explaining why the amplitude continually increases. 
To prove this last statement, multiply the equation of motion 



bx 

a -J- bt 


+ c 2 x = 0 


by 2i and integrate. This gives 



x 2 

a bt 




Taking a, /J as times of temporary rest, not necessarily consecutive, 
we dispose of the first term. The integrand is certaiuly positive and 
hence x is greater at ft than at a. The oscillation increases in amplitude, 
like a vibration picking up on resonance. 

Another curious feature of the motion is that the times for half¬ 
swings inward to the origin continually increase; but the times for 
half-swings away from the origin continually decrease. This is the 
problem referred to in 4-6, Ex. 23. To prove it, let 0, c v c 2 , &c., denote 
the values of cz for zero velocity, and let d v d 2 , &c., denote the values 
for zero displacement. We then have 

d r + 1 d f 77, C r C r _ j < 77, 


whence by subtraction 

(d r+l - c r ) - (d r - c,^) > 0, d r+l -c r > d r - c r . v 


Similarly from 

C r - C r _j <77, d r — d r _ j > 77, 

we deduce 

( c f — d r ) — (c r _i — c? r - 1 ) <0, c r — d r < c r . 1 d r _ v 
A rough sketch shows that this establishes our contention. 


EXERCISES 

1. Derive the acceleration by differentiating the velocity and check by the 
equation of motion. 

2. In dimple harmonic motion, when the displacement is zero the velocity 
is maximum and the acceleration changes sign. Prove that in the present problem 

tbl A^o C ln°si.nple harmonic motion the displacement and the acceleration reach 
their maxima simultaneously. Does this hold in the present case / 
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3. Determine the constants if initially the body is passing through the origin 
with velocity w. 

4. Accor ling to the principles of dynamics, the momentum of t lie body 
when passing through the origin equals the time integral of the force since the 
previous rest position. How does this work out in terms of the Bessel functions? 

5. Examine the problem when the mass at time / is (a -f- bt)~ 3 . 

6. Prove that the kinetic energy of the body when passing the origin is pro¬ 
portional to <x{AJ 0 {a) -f BY 0 (<x)} 2 , where a is a value of c: corresponding to x — 0. 


6-7. Longitudinal vibrations of a bar. 

Problem 6 .—Discuss the longitudinal vibrations in a tapered bar. 

Imagine a long thin metal bar to occupy a finite part of the r-axis 

from the origin. As the result, of a longitudinal blow, the section A 

which was originally at distance x is displaced to distance x -f- u. 

Note that u is not necessarily small since the bar may have moved 

bodily. Similarly the section £ : originally atr + Sr,'is displaced a 

distance u + 8i/. The extension of the element 8x is thus 8m, so that 

t he strain is du/d:r and the stress is E Bu/dr, where E is Young's modulus. 

The convention of sign is fixed by the stress being tensile when 8m is 
positive. 

If a be the area of cross-section A, and F the force acting on it 
we have * ’ 

OX 


Similarly the force acting on the section B is F 4- 8F 
to the right is ^ 



The resultant 


The mass of the element AH is pabx if p « the density and ns its 
acceleration is d*u/Bt* we have the equation of motion 



=S-HSS-£7.=5S 


cl* 3 * 2 ’ 



p 



sc. 
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T'his equation has solutions of the type 



sin 


nx 


cos 

sin 



The choice of sine or cosine is determined by the mode of support. The 
length of the bar then determines n; and the initial conditions deter¬ 
mine A. Thus if both ends of the bar of length A are fixed, so that 
u = 0 both when x = 0 and when x = A, a suitable solution is 

u = Z sin nx(A cos hit -f- B sin lent), 

provided that sinnA = 0. This gives nX = t t, 2t t, 3tt, .... The per¬ 
missible values of n are thus determined, and the constants A, B are 
chosen to suit given initial conditions of motion and displacement. 

The above cursory discussion does not involve Bessel functions and 
fuller treatment on these lines is to be found in texts on sound. If we 
abandon the uniformity of section whilst still treating E and p as 
constants, the equation of motion becomes 

p chi _ 1 3 | du\ 

E dl 2 a ?x \ dx) 

On the grounds that longitudinal vibration is possible we write u = 
X sin (pt + <f>), where A' is independent of t. The equation becomes 

^ + + = o, = 

dxr a dx dx E 


If the bar is slightly tapered we can write 

( x\ da a l da _ 1 

1 + a/ Ji-p 


whence 


1 dX 


dr-x 


dx 2 x + P dx 


+ k 2 X = 0. 


With the slight modification 


we have 


so that 


x -f p = z, dx = dz, 

™ + l d I + i* X = o, 

dz 2 z dz 

X = AJ 0 (h) + BY 0 (kz). 
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As an illustration we suppose the left end to be fixed and the right 
end to be free. At the left we have the conditions 


u — 0 — X = x, z = 

Hence 

0 = AJ 0 (kp) + BY 0 (kB). 


At the free end there can be no stress and we have the conditions 


whence 


x — A, z = A -f ft, 



dX 



0 = AJ x {k{ A + 0)} + BY X {1{\ + fi)}. 


The free end of the bar, a place of zero stress, is a place of maximum 
movement, a loop. The foregoing is in line with the fact that the 
maxima of J 0 and Y 0 are located by the zeros of J x and Y, respec¬ 
tively. The elimination of the ratio A /B gives 

«W) - Ji {*(A + ft} 

y 0 m yah a + My 

The first root of this equation in Jc corresponds to a stress distribution 
which w part of a loop of the curve AJ x (kz) + BY x (kt) t with a zero 
at the right. The higher roots give distributions with one or more 
stress-nodes, there being always one at the right. In each case the 
stress at any point varies sinusoidally with the time. 


EXERCISES 

J *W) + P)> 

y o(*3) Y„{*(X+0)}* 

bt d^T iS “ = * (1 + pn>ve l "“‘ —*» °< tho 

J i(*p) _ + 3)} 

X 0 {I(X + P)}‘ 

Tibmtions h ••**— ■» 

■ number, pusitivl t u^ve = “<> + */W. whom 

Investigate the case where m = 
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6-8. The vibrating membrane. 

Problem 7 .—The vibrations of a stretched membrane. 

A concise history of English literature could be pardoned for omit¬ 
ting the name of John Clare; but the complete omission of Shake¬ 
speare could hardly be condoned. Similarly from a survey of the 
relevant literature it appears that even the most cursory treatment of 
Bessel functions cannot avoid the vibrating membrane. It was first 
discussed nearly two centuries ago by L. Euler in 1764; it has been 
heavily belaboured in innumerable textbooks ever since. 

The reader is aware that a stretched membrane, such as a drum or 
tambourine, is capable of vibrating when tapped. The problem is the 
two-dimensional analogue of the vibrating string and similar assump¬ 
tions are made for effecting a solution. A membrane differs from a 
disc as a string differs from a rod in that flexural rigidity, shear and 
bending are taken as negligible. The vibrations are presumed due to 
the tension applied to the material. 

The equilibrium position of the membrane is taken as the hori¬ 
zontal plane and gravity is ignored. During movement let z be the 
upward displacement of a point whose co-ordinates are otherwise 
r, 9. Consider the element limited by the arcs of radii r and r + 8r, 
and the radii defined by Q and 9 + 8 9. If T be the tension per unit 
length along the edge r80, the force on this edge is F = Tr89 acting 
at an angle 0 to the horizontal. The vertically downward component 
of this is 

02 

Y = F simp, sirn 


At the opposite edge the vertically upward component is Y -f 81 and 
the upward resultant is 



Along the edge Sr we similarly have 

F l = T 1 8r, Y l = F 1 sin$ v 



The upward resultant from the two Sr edges is 

3ft d \T l 
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If p be the surface density, the mass of the element is prdrdd and its 
upward acceleration is d 2 z/ct 2 . As in tlie case of a string, and for tlie 
same reason, it is customary to take the tension as constant; this 
gives the equation of motion 


3 ( ds\ 1 dh = prZh 
3rl3r/ r30 2 T cl 2 ’ 


As z is presumed to undergo periodic variations, expressible in a 
Fourier series, we write 

dH 

z = w exp (ipt), — = —p 2 w exp (ipt), 


where w is independent of I. This gives 


3 3w\ 13 ho 

dr r 8? + ^ + ^ = 0 ’ * 


P-P 
T ' 


We now further assume u> = u cosnd, where u is independent of 6 
and n is an integer. The significance of this assumption will appear 
later. We now have 




of which the solution is 

u = AJ n {kr) + BY n (kr). 
Accordingly, for the displacement, 


cos 


2 = sin pt cosnd { AJ nW + HF n (fe-)}. 


So far nothing has been said about the shape of the membrane; 
but the analysis is obviously adapted to the circular form, and this 
will be employed. If the membrane is complete up to the centre, as 
in a tambourine, the function Y n must be discarded since it makes 
z infinite at the origin. It may be retained for an annular membrane. 

in the simplest case, where n is zero, we can take z ~ A cos ptJJkr) 
As 2 ,s permanently zero at the edge, where r = a, we have JJka) L o' 
The smallest permissible value of ka is 2-4048. For a given membrane, 
this determines k and thence p and the corresponding frequency. 

It is known that J 0 {x) has an infinity of zeros. We denote the suc- 
ession by c 1 , c 2 , . . . , and the corresponding values of Irby k, k 2 

80 (o UO) = C *' A CUriOUS consef l uence of this can be illustrated' by 
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giving s the value 3, say. If 2 is proportional to J 0 {k 3 r), then at the 
boundary where r = a we have 2 proportional to 

= J 0 ( c 3) = 0 . 

But nearer the centre, where r = kjtt/k^ we have 2 proportional to 

•^o(^y) = «/ 0 (M = ^1) = 0. 

The conclusion is that this value of r gives a nodal circle; and the 
same applies if r = k<p/k 3 . In general there are (s — 1) nodal circles 
if J 0 (k/i) = 0. This is the analogue of the possible nodes on a vibrating 
string. 

In the more general case where 2 is proportional to cos ndJ n (kr) 
the displacement is permanently zero if either J n (kr) = 0 or cos nd = 0. 
The former must hold at the periphery and J„(ka ) = 0. This has an 
infinite number of roots, each of which gives a corresponding value 
of k. There are no nodal circles, except the boundary, for the smallest 
root; the higher roots give nodal circles. The alternative cos nd = 0 
gives nodal diameters corresponding to nd = \tt, § -n, .... 

In the above solution, all points simultaneously pass through the 
equilibrium position and hence have the same period. This is known 
as a “ normal mode in general the motion is more complicated 
than this. In a normal mode, with nodal circles and diameters, adja¬ 
cent sectors are in opposite phase. The matter is treated more fully 
in texts on sound. 


EXERCISES 


1. If an annular membrane has inner and outer radii a, b respectively, prove 
that its slowest mode of vibration is given by the smallest root of the equation 

J 0 (ka) _ J 0 {kb) 

Y 0 (ka) ~ Y 0 (kby 


2. Reverting to the problem of the taut string, a non-uniform string of length 
>. has its ends fixed and is kept taut by a tension T. The line-density at distance 
x from one end is p( 1 -f- kx)~ l . With the notation of the text, prove that the 
possible frequencies are given b} 1 


J i(c) _ Ji(nc) 
Y x {e) ~ Yline)* 


n = (1 + *X)», 



3. Investigate the more general case of the vibrating string when the law of 
density is p(l + kx)~ a , a = 2 — to -1 , m an integer. 



THE SECOND SOLUTION; FURTHER APPLICATIONS 


91 


6-9. The oscillating chain. 

Problem 8 .—The small oscillations of a uniform chain suspended at one 
end. 

This is another of the classical problems of our subject. It was 
first discussed in 1732 by Daniel of the Swiss family Bernoulli; later 
on in 1781 it was taken up by Euler. The chain is a flexible filament 
of uniform line-density, devoid of flexural rigidity. Let one end be 
fixed to the origin 0, with OY to the right and OX downwards. If A 
be a point defined by OA = s, the tension T at A makes angle ip with 
OX and has an inward component given by Y = Tsini/;. At an 
adjacent point B defined by s -J- 8s the corresponding outward com¬ 
ponent is Y -|- 8Y, so that the outward resultant on the element 
AB is 



If p be the line density, the mass of the element AB is pBs and its accel¬ 
eration is y; moreover, since the oscillations are small, it is sufficiently 
accurate to take the tension T as the weight of chain below A. Hence 

T = P9i* ~ «). w here A is the length of the chain. This gives as the 
equation of motion 


H dt- ds 



The oscillatory nature of the motion can be allowed for by the sub¬ 
stitution y = u cos { pt -f- tf>) and the equation for u becomes 


(A-s) 


dru 

ds 2 


du p 2 

T 4- — M = o. 
ds g 


The substitution \ — s = z, ds = — dz, which is equivalent to measur- 
ing from the lower end of the chain, gives 


, 1 du 

dz 2 z dz z 



o'solution is obtained by comparison with our general 
3-6(2) and yields 


standard 


Thus 


a = 0, y = \, n = 0, 0 = 2 k. 

u = AJ 0 {2kz*) i- BY 0 {2kz*). 
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For the determination of the constants we examine the free end. Here 
5 = A, 2 = 0 and as the displacement y (and therefore u) is to be small 
we must have B zero. The other constant remains indeterminate since 
the oscillations are merely small. At the fixed end Owe have y = 
0 = u = s, 2 = A so that J 0 {2kX i ) = 0. Thus 2£A* is a zero of J 0 (z); 
any particular zero furnishes a value of k and this in turn determines 
j> and the corresponding frequency of oscillation. 

If we equate t to zero to get the initial shape of the chain we conclude 
that, for the simple type of oscillation we have been discussing, the 
chain must be started in the form of an odd number of half-loops of 
the curve J 0 (2kz-) with a zero at 0- Bv slightly changing the origin 
of time the system can start from rest; but the form will still have to 
be as stated. 


EXERCISES 

1. Examine the ease of the non-uniform chain where the density is propor¬ 

tional to the mth power of the distance from the free end. Prove that the solution 
depends on the function $rg = 4p*(m + 1). It was pointed out by 

Greenhill that a practical approximation would be a large number of parallel 
wires of non-uniform length connected on the principle of the Venetian blind. 

2. If the end of the chain carries a mass, equivalent to a length h, prove that 
the normal modes are given by 

J.(2lhi) J 0 (2LH) 

Y.(2khi) ~ Y0(2111)' 

where / = X -f- h. 

610. Heat conduction in one dimension. 

One of the outstanding characteristics of physical constants is 
that almost inevitably they tend to show variation over an extended 
range. Thus a coefficient of expansion which normally has the con¬ 
stant value a is pretty certain to need modification if it is to be em¬ 
ployed over a wide range of temperatures. The form to be adopted is 
then empirical, and even if theory indicates the nature of the depar¬ 
ture from normal, experimental verification would be needed before 
adopting any specific form. The form most often adopted is a simple 
polynomial a -f /fc + yt 2 , the number of coefficients being a measure 
of one’s fastidiousness, and their values being determined by some 
such method as the principle of least squares. For the most part the 
above method works quite well and it is usually possible to derive a 
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The area a being the same at both ends of the element, we deduce 

d-9 13 k d9 ps d9 

dx 2 k dxdx k dt 


as the equation of heat conduction. 

We attempt to effect a solution by the substitution 9 = X exp (—pi). 
Here X is independent of t\ the negative exponential is preferred 
since a temperature may reasonably be expected to decrease to zero, 
but not to augment indefinitely. The equation becomes 


d 2 X IdKdX 
dx 2 k dx dx 


Wx = o. 


In the simplest case the physical coefficients are all taken to be 
constants. The middle term then drops out and the equation is soluble 
in trigonometrical functions. As a first modification let us assume p 
replaced by fp(x f)~ x . This variable density might be ascribed to 
any of a number of causes such as porosity, humidity, weathering, 
and so on. The coefficients being otherwise constant we now have 


The substitution 
gives 


d 2 X 
dx 2 




P s Pf 

K 


f -\- x = z, dx = dz 


™+*x = o. 

dz 2 Z 


This is soluble by Bessel functions. Comparison with our standard 
3-6(2) gives 

a = h, = 26, y = h n = 1. 

The solution is 

X = Az'J x (2bz') + B 2 »y i (26= i ). 

The temperature is therefore 

9 = Z *e-^‘{A J x (2bz*) + BY 1 (26z*)}. 

The constants A, B would have to be determined to fit pre-assigned 
initial conditions, and no doubt a series of values of p would have to 
be employed. This aspect of the matter will be resumed later. 




THE SECOND SOLUTION; FURTHER APPLICATIONS 


Vo 


EXERCISES 

1. Calculate the temperature gradient. If one face, defined by z = is per¬ 
manently at the temperature of the external medium, deduce that 

AJ U (20)A) + BY 0 {2b\i) = 0. 

2. If the variable density is pP(x + f)~- and the variable specific heat is 
sA -, (x + h), deduce that, with the notation of the text, 

X = zi{AJ n V2bzl) + £F n (26 : »)}, 

where 

b-Kk = psppt 1 + 46 2 c = n\ c = f — h f 

so that the order n is real provided 1 -f 4fc 2 c > 0. 

3. Presuming the density and the specific heat to be constants, let the variable 
conductivity be k(1 -j- x/cz) r . 

(i) Prove that the equation is soluble in terms of J 2i Y 2 if r = 3. 

(ii) If r = 1 the functions are of zero order. 

(iii) If r = 1J the functions arc of the first order. 

(iv) Prove that the equation is soluble by Bessel functions for all values of 
r except r = 2, in which case the solution can be effected by elementary means. 

(v) If r= 2 the equation may still be soluble by Bessel functions if p and 
s vary separately or jointly. Investigate this. 

4. Investigate the flow of heat along a lagged bar of uniform section on the 
assumption that the conductivity, density and specific heat arc proportional to 
*■% x* 1 and & respectively. 


611. The tapered strut. 


Problem 10. —Discuss the stability of the tapered strut. 

An elementary form of this problem has already been considered 
m the previous chapter. We continue with the case of the triangular 
ribbon first removing the tip and thus reducing it to the trapezoidal 

form. This gives I = H(x + h)/h. The number /, defines the fictitious 
vertex of the triangle. The equilibrium equation is 


or, 

d?y , k 2 

dz 2 + 4 2 y ~ °’ z ~ x + h > £ 2 : 
This has previously been solved and we have 


4 Ph 
EH' 


y = *{AJ x (kz') + BYJkz*)}. 
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The end conditions being 

x = 0, 2 = k, 

x = A, z = h + A = 1, 

we derive 

0 = AJ x (kh') + BY x (kk>), 

0 = AJ x {kl*) + BY x (lil k ). 

The ratio A/Ii is determinate, but one of them remains arbitrary. The 
elimination of the ratio gives 

Jtfh^Y^kl*) = J X (W)Y X (W) 

as the transcendental equation which determines k ; this in turn 
gives P. 

As a somewhat less simple illustration we take the left end to be 
clamped horizontal. The origin is at the left and the y axis is positive 
downwards. The right end is under horizontal thrust P, but is pinned 
to prevent lateral displacement. This calls up a lateral force which 
we may call R and presume to act downwards. An equal and opposite 
thrust P acts at the clamp, as also an upthrust R. There is further a 
couple C which is coimterclockwi.se on the assumption that the deflec¬ 
tion is downward. Consideration of the moments about any point 
of the ribbon with co-ordinates x, y gives the equilibrium equation 

Ely " = C — Rx — Py. 

With the above assumption as to the value of I we have 

Eu ^ + ’>) dry^ c _ Rx _p y 
h dx 2 

It is statically evident that C = RA, and hence 

z + h-z)~ Py, z = * + h. 

h dz- 


There is now a particular integral 

y = R(A + h - z)/P. 


The reduced, or auxiliary, equation can be written 



, * 2 n 

H— y = o, 



APh 

EH’ 
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The solution is 

y = + ? (A + h- z), 

whence 

% = %= - 1 
For the end conditions at the left we have 

x = 0, z = h y y = 0 = y'. 

Hence 

0 = AJ,m + BYJkhl) + 

0 = AJ e «h>) + BY 0 (U,>) - 

Pk 

At the right we have 

x = A, z = A h = ?, y = 0. 

Hence 

These three homogeneous equations can determine nothing more than 
the ratios A : B : R even if consistent. It is readily verified that 

_ ^ __ P 7 tR 

2A*y,(AA‘) H- A-Ay 0 (AA*) 2h*J 1 (kh*) + kAJ 0 (kh*) ~ W 

The condition that the three homogeneous equations are consistent is 
obtained by eliminating the ratio A : B. The result is evidently 

+ kAY 0 (kh>) YAW) 

2h*J Akh*) + kAJ 0 (kh*) JAW)' 

This is the transcendental equation that determines k, which in 
turn determines the critical thrust P. 

It is possible to locate the point that tends to maximum deflection 

^• J< ^ Ua *i n ^^n t0 ZCr0; but tbe manipulative algebra involved is more 
te<hous than illuminating. The lateral force R remains indeterminate 

iind ? cur y e taken the nbbon is reminiscent of the profile of the 
underside of a spoon. Consider the function R{X + h — z) — P„ It 

£ zero a fc the right; at the left it is positive but decreasing. It accord- 

g y as a zero somewhere on the ribbon and it is not difficult to see 
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that this locates the inflection. The fact that it vanishes at the right 
is equivalent to the absence of any end-couple. 

Reverting to the case of the ribbon under simple end-thrust, let us 
take 

7 = "fer*)’= "O' = = * + 

The equilibrium equation takes the form 

y " + m *~ r y = / * r "‘ 2 = £i { 

On comparison with our standard 3-6(2) we deduce 

a = i, y = 1 — ir, n = —-—, fi =-. 

2 + r -/ 

It appears from this that the problem is soluble in terms of Bessel 
functions for all values of r except —2. As r falls below —1, the order 
n augments without limit. When r is zero we have a = I = n, so that 
the solutions are of the tvpe i'-J j. But when r is zero, the ribbon is 
uniform and the solution is known to be trigonometrical. We thus 
revert to the theorem that functions of order half an odd integer are 
expressible trigonometrically. 

A more interesting result is derived from letting h become inde- 
finitelv large, so that r also becomes large. Since Br is proportional to 
h~ r , we have /J; 7 ultimately proportional to h , so that the function 
J K (p r 7 ) approximates to J«(h). But as h augments indefinitely, the 
ribbon approaches uniformity, irrespective of the value of n. The 
solution is then trigonometrical. We conclude that, for large values 
of the argument and irrespective of the order, the function «/«(*) 
between consecutive zeros behaves like zr*sinz. The result is one 
aspect of what are known as the asymptotic formulae; it marks an 
advance on our previous knowledge of large values, which was (a) that 
the absolute magnitude of the stationary value continuously decreases, 
and (6) that xJ x (z) increases indefinitely with x. 
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EXERCISES 

1. Taking I = H(x -f- A)/A for a ribbon of length X with both ends built-in 
and subjected to end-thrust P, prove that the condition of stability is 

. , 2 _ 4PA. 

y 0 (AAl) Y 0 {W)’ EH’ + 

2. The ribbon of length X is built-in at the left. At the right is vertical load 
W and horizontal thrust P. Taking 

'X+ A\*/3 




prove that the conditions at the left lead to 

0 = AJ a (p) + BY 3 (p) + 

0 = AJ 2 (p) + BY dp) - 


where 


C 

Phi' 

61FA» 

mP 


m* = 


3CPA 5/3 
Ell 


, m 6 h = p «. 


The end condition at the right leads to 

0 = AJ a (q) + BY a (q), mH =q*, / = X + A. 

The last equation gives the ratio A : B, the other two equations then determine 
A and B in terms of C. Thus C is determinable, and so is the end deflection 8 
given by C = IPX -f PS. 

3. Establish the following results: 


(*) f Ji(x)dx=l. 

Jo 

(») f «/„(*)</* = rJ n+t (.r)dx. 

JO Jo 

J r 00 

3TiJ i (x)dx= (Jtt)*. 

0 

(iv) f X~ l J a {x)dx = 

(v) r^j n+Md x= 


2"r(n + 1)‘ 


By tto ropMkd «° [ tb °— 
hj\{x)dx = J l{x) + J i(x) + J t(x) + m m m m 

Hence by allowing « to tend to infinity the value of the integral is zero. 
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The Modified Functions 


7-1. Function with imaginary argument. 


It frequently happens in physical investigations that we need 
Bessel functions of purely imaginary argument. There need not on 
that account be anything imaginary about the functions, any more 
than there is anything imaginary about cosx, which happens to be 
definable by means of the imaginary exponential exptx. These func¬ 
tions satisfy a differential equation which is a modification of Bessel’s 
equation different from any we have yet considered. In a more ele¬ 
mentary field we have analogously that cosa >t is a solution of the 
equation x + arx = 0, whilst coshwf, or cos iwt, is a solution of the 

allied equation x — ai 2 x = 0. 

We know that a solution of the equation 


^ 4 - 1 ^ 4 - 
dx- x dx 



is 

JnW 


(i*V (j (**)* ■ <**)- _ 

f(» +1)1 l(n + l) 1 1 . 2(» + l)(n + 2) 


Replacing x by ix and dx by idx, we conclude that 


J„(ix) = 


(F*)" (i^) 8 

r(n + 1) l l(n + 1) 


+ • 



is a solution of the equation 



dx 2 x dx 



and so. too, is any numerical multiple of it. Note that only one of the 

four terms differs in sign from Bessel’s equation. 

We remove the undesirable imaginary in defining a new function 

by the relation 

(2) /„(*) = 

100 
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whence 

( 3 ) /„(*) 

This function I n {x) is known as the modified Bessel function of the 
first kind, and it satisfies the equation (I). Certain of its properties 
strike the eye immediately (fig. 8) and show that its behaviour is 
quite different from that of J n (x). For positive values of x, which are 
the only ones that matter in practice, and for positive orders, every 



Fig 8.—March of I H and K n 


_ (W f, , (j*) a 

n / i i > 1 » < 


r(»+i) 


rr + 


(k) 4 


1(»+1) 1.2(» + !)(» + 2) 




brack f is P° s r itive “d increases with x. We conclude that 
in these circumstances /.{*) cannot have a positive zero and the 
external factor (Ji)» shows that the function passes through the origin 
An unmedrate corollary is that no solution of ?he equation (4 can have 

TLTl ? the 6eneral the ° rem the interlacing of 2 er OS 
would then have a zero between them, and this we knfw is false' 
Ihe unique zero order gives 

(4) 4 ( 0 ) = 1 = j 0 ( o). 

!a p ;S z°et:nd h t e he e T rn t al (W ‘ bM “”“ “ wb ™ 
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integral orders we write the relation J „ = (—1 )"</_„ in the form 
(i)-"J n = ( i) n J-„, which proves that 

(5) /„(*) = 

Note that there is now no ambiguity of sign; the functions are identical. 
We shall accordingly be faced with the problem of finding an indepen¬ 
dent solution for integral orders. The series definition of the function 
partially breaks down when n is a negative integer. Some of the deno¬ 
minators become zero and the factor T(n -f- 1) becomes infinite. The 
reconciliation is then the same as was employed in the parallel case of 
the function J- n (x). 


7-2. Recurrence formulae. 

The Bessel functions of order £ are known to be expressible trigo¬ 
nometrically, and on switching over to an imaginary argument it is 
reasonable to expect the modified functions of order £ to be expressible 
in terms of hyperbolic functions. This is the case and the verification 
is left to the reader among the exercises. Similarly it is natural to 
expect recurrence formulae, and in fact we might have investigated 
the functions from this end, as suggested in 3-6, Ex. 8. Replacing 
x by ix in the relation 

x ^ J n {x) = nJ n (x) — xJ n+1 (x), 
dx 


we have 


xx 


d(ix) 


J n (ix) = nJ n (ix) - ixJn+iifr) = x ^ 


Multiply this by (*)-" and it reads 

x I„(x) = nl n (x) -f- xl n + j(x). 
dx 

Note the change of sign in the last term. It is then a simple matter to 
deduce the five relations 


(1) 

Xl' n -f- w/ n — Xl n - j 

(2) 

Xln ~ */„ = *In+ 1 


2n . 

(3) 

I n—1 * n+l ~ * n 

(4) 

l n -l + Ail-1 = » 

(5) 

io = A- 
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It lias already been pointed out. that- for integral orders, every 
terra in the series for /„ is positive and increases with x. The same is 
therefore true for and Hence the grapli rises from the origin 
ever more steeply, rather like the graph of sinh x. Tiie relation 



indicates that 7 n+1 lies lower than 7*^ (fig. 8). For a given x, the value 
of the function decreases as the order rises, so that the higher orders 
are slower in taking off from the x axis. We can illustrate this from the 
tables by considering, say, 7„(4) for different values of n. We find 

w » 012345678 

7„(4), 11-30 9-76 6-42 3-34 1-42 0-50 0-15 004 0-U1 

Two solutions of the same second order equation are necessarily 
connected, and the appropriate modification of the relation 4 - 4 ( 1 ) 

J n J. n ' - J n 'J_ n = _ 

7TX 

gives 




2 sin mr 


f 

7TX 


whence other similar results can be derived. 


7*3. Standard form. 

The functions that appear in the solution of various problems can 
b ; to ‘ake the simple form /.<„ and J aIXg“ 
reqmre a general standard for comparison. If we replace B by ifi in 
our previous standard 3-0(2) we conclude that Jj liBxA or anv 

mUl ‘ ip,e ° f " hich indudcs » a solution of the 


( 1 ) 


&y , 1 - 2a 

dx^~^r~ 


% - [W*- 1 )* + 


2 2 o 

- n *• 


n~y 


X- 


y = o. 


I the^equatio °“ ° f Ieference - •*•** appears that 

equation is in thl ummal forn th n °Ti a ‘ f ° rm| COnveree ^ the 

normal form the solution must contain the factor x* 
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7-4. Lommel integrals. 

We can work out the corresponding Lommel integrals by putting 
a = i, y = 1. Taking two functions of the same order, let 

u = x*/ n (Az), v = £*/„(/*£)' 


The corresponding normal equations are 




where dashes denote differentiations with respect to x. 
of n leads to 


which is 


(A 2 — ix 2 )\uvdx = vu' — v'u, 


The elimination 


(1) (A 2 —/x 2 ) rxI n (Xx)I n (^x)(1 x = x{XI tl ( l xx)I n , (Xx)-ixJ n , ( t ix)I n {Xx)}, 

-'O 

where the dashes now denote differentiation with respect to the argu¬ 
ment. No suitable upper limit exists for making the right side zero; 
the function /„ accordingly has no orthogonal property. 

If we take the two functions to be of different orders with the same 

argument, so that 

u = **/„( Ax), v = x*/ m (Ax), 


we have the normal equations 





From these we deduce 



which is 


(2) (»* - dx = A{7.’(A*)/„(A*) - /.(A*)«**>}• 


X 
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Finally, if we multiply the equation 

»"+j!=( A2+ ?> *= 7 " (A *>- 

by 2 x~y' and integrate, we derive 

ay* = (A 2 x 2 -{- n 2 )*/ 2 — 2A 2 Jry 2 (/r, 


which is 


2 Jx{I n (Xx)Ydx =(* + £){/„(A*)}= - *2{/„'(Ax)}». 


EXERCISES 

1. Check one of the recurrence formulae against the following values 

J 3 = 2-0661 / 4 = 0-8104 / 4 = 0-2651 x = 3*6. 

2. Prove that 

1 »(*) = (~)* 8,nh;r ' f-ifr) = (~)* coshz * 

Evaluate for the orders ± 1|, ± 21. 

Enunciate and prove the general proposition. 

o _ 2 sinmr _ r T r T 

* ^ •*»*-»+1 


4 * U*)^ - * n - r/n - 

ts) r &) = ^ 

^ {ari«J n <p*l)> = iparKn+D/^^p*!). 

6. Prove 1 0 (x) is a solution of the equation — = xt/. 

6. x*// = {„(„ _ I) + *a } / n _ xIn+v 

7- (i) **/«_, = «/ n + {« + 2)/„ +t + („ + 4)/ n+4 + . . . . 

(ii) tf/n-'M 11 * = A.~ / n+2 + / n+< - . . . . 


^ ^n-r + r A,-r+2 + 


da? 
(0 150 ) 


r(r- 1) 

1.2 


n-r+4 


+ • • * + fn-lr* 


8 
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9. Prove from the recurrence formula that / _ n (z) has no positive zero ii 
0 < n < 1. 

10. Obtain the expansion for I n (x) by using the method of Frobenius to solve 
the differential equation. 

11. Show that the second solution for zero order contains the term I 0 {z) logx. 

12. Determine whether the equation xrf = y can be solved in terms of 7 n . 

[Ans. x*/ 1 (2xl).] 

13. If jx is a zero of J n (t) 9 prove that 

(X 2 -f p 2 ) f xI n C>jc)J n ([Lx)dx = (i/ n (X)J n+l (jx). 

Jo 

14. Prove that 


(— = - 

J */„*(*) 


* I -n( x ) 
2 sin nr: I n (x)' 


15. In view of Ex. 2 above, prove that if the order is half an odd integer, 
positive or negative, the function I n (x)er x for large values of x is asymptotic to 
(2wx)“l. An appeal to the principle of continuity makes it reasonable to suppose 
that the proposition holds for all orders. 

16. If y satisfies the equation 7-1(1) and a new variable z is defined by y =«*, 
prove that z satisfies the equation 


d*z ( \\ dz 

d^ + \ 2 + x)di + 


- t ! 2 



7-5. The second solution. 

The complete solution of the equation 7-1(1) is y = AI n {x) + 
BI_ n (x) where A and B are arbitrary constants, and no solution can 
have any form but this. There is effectively only one constant here when 
the order is integral, and in almost all research that involves Bessel 
functions the advent of integral orders is inevitable. The problem of 
finding the second independent solution for integral orders can be 
approached in various ways, each with its own disadvantages and 
compensating advantages. We could employ the method of Frobenius, 
with necessary modifications; alternatively we could define the func¬ 
tion by Y n (ix); or we could use the equation 

y=Ux) f^£k?' 

We fall back on a previous line of argument instead. If the equation 
7 -1(1) be differentiated with respect to the order n, we have 

dx*\dn ) + xdx \3n J \ *V ** 
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Similarly since is a solution we have 

£ PH + - t ( t 1 -"} - ( 1 + "0 %- - 2 ” 7 

ax- [cn I x ax \ cn ) \ x-J cn x z 


If for brevity we put 


we have by subtraction 




When n is integral the right side is zero and we conclude that K is then 
a solution of equation 7-1(1). The procedure is analogous to that for 
obtaining Y„ from J„. 

What is required is a definition that holds for non-integral orders 
and which is equivalent to K for integral orders. A selection, due to 
various writers, is available and we may adopt 

(1) K„(x) = in ! — " ~ 7 " , 

81117 * 77 - 


which is due to Macdonald. This being a linear combination of /„ 
and is certainly a solution for non-integral orders. When n is an 
integer it, takes the indefinite form 0/0; but the usual procedure of the 
differential calculus gives 

d . 

— sin 7177 - = 7 T cos nn = n( — ) n , 


so that 


K n (x) = *( - 



This is the modified Bessel function of the second kind. It is rather a 
tedious business to find its explicit form for integral orders and the 
result is not very informative when acquired; the properties of the 
function are not ascertained from its series. 



7-6. Recurrence formal®. 

The function naturally satisfies recurrence formulae but they are 
different from any yet encountered. This is one of the disadvantages 

Cl r enti r R ed ; , Wh ° reaS J " and 7 " have the same recurrence 
formulae, the modified functions l n and K n as accepted herein do not. 
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It must be borne in mind when consulting other works that slight 
differences may make their appearance, due to the particular defini¬ 
tion adopted for K n . 

If we replace n bv — n in the definition 7-5(1), both numerator and 
denominator change sign, so that 

(1) K n (x) = K_ n {x). 

The companion formula to 

Xln = »I H + Xl n +, 

is xln = — nin 4 xln-V 

Changing the sign of n in the latter we have 

xl-n = nl-n 4 xl_ n . v 

Subtract the first equation and divide by sin nn. In virtue of the rela¬ 
tion sin(n -j- 1)7r = — sinurr we have 

(2) xKn = nK„ — xK ni . v 
Changing the sign of n and using (1) we have 

(3) xK n ' = —nK n — xK n _ v 

Subtraction then gives 

(4) xK n+l = 2nK n + xK n -v 

whilst addition gives 

(5) 2A' n ' = A' n+1 K n -v 

These have been established on the assumption that n is not integral; 
the justification for extending them to integral orders is an appeal to 
their continuity as functions of their order. 

Since I„ and K n are two independent solutions of the same reduced 
equation of the second order they must be related. Using the defini¬ 
tion 7-5(1) we have 

I n K„ — I n 'K n = ^77 COSCC mr{l n (I-n ~ K') ~ Lx^-n — J n)} 

= ^7T cosec N7r{/ n /_ n ' — 

, ( 2sinn7rl 1 

(6) = &7T cosec iiTT !--—j = - 

which can be thrown into various forms by use of the recurrence for¬ 
mula?. 
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When the order is zero we have from the recurrence formulae 

(7) AY = -A\ = ~K_i. 

When the order is i we have from the definition and from the previous 
7-4, Ex. 2 

( 2 V 

(cosh x — sinhx). 


( 8 ) 

Since 


= (£) >e '= K -^- 
K„ = 1 -K i + K. l =(l+^K t , 


it is readily deduced that K n , when n is half an odd integer, is A’ } 
multiplied by a polynomial in x~ l . Consequently all such functions are 
asymptotic to the x axis like x~*c~ x , and an appeal to the principle of 
continuity warrants the belief that the function behaves similarly 
whatever the order. This is another aspect of the asymptotic values 
previously mentioned in the concluding remarks of the last chapter. 


7-7. Graph of K n . 

Coming now to the graph of K n , it has to be frankly admitted that 
satisfactory proofs of an elementary nature are not to be had, especi¬ 
ally of the outstanding proposition that K n (x) has no positive zeros. 
The graph resembles a Boyle’s law curve, or the rectangular hyperbola 

xy = constant in the first quadrant. In default of rigidity the follow¬ 
ing is offered as a second best. 

The relation 7*6(6) can be written 


dx\lj 


xl „ 2 


The right side is essentially negative and steadily decreases in absolute 
value to zero as a: approaches infinity. The function KJI n has no 
stationary values; its graph ultimately becomes horizontal and the 
runction approaches a constant value which we may denote by c 
lhis c must be zero; for since /„ ultimately increases without limit, so 

^ > 15 r th i" g but zero ' But K » certainl 7 does not in- 

ZZnlt ° rdCr iS ha ‘ f an ° dd inte ® er - and 

easonable to think it does not do so in any case. Hence we take c to 
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be zero. And as /„ is always positive, the conclusion seems to be that 
K n lies wholly in the first quadrant. 

Corroboration is afforded by the previous 7-4, Ex. 11. If K 0 comes 
down from positive infinity, we have K 0 ' negative and therefore K x 
positive. The recurrence formula 

xK 2 = 2 K x -f- xK 0 

then shows that K 2 is positive, and so on. The slope is given by 

xK x ' = ~(K l + xK 0 ), 
xK 2 ' = -(2 K 2 + xK x ), 

and so on, so that they are all negative. The relation 

= A.-, + - K n 

X 

shows that, for a given x, K n+1 is greater than K n _ x and the function 
increases in value with the order (fig. 8). The following illustration 
taken almost at random from the tables shows the growth of If„(2-2) 
with the order: 

n, 0 12 34 5678 

K n , 0089 0-108 0-187 0-449 1-411 5-578 26-77 151-6 991 

Looking at the same functions with the doubled argument 4-4, the 
figures are 

K n , 0-007 0-008 0-011 0-018 0-035 0-081 0-219 0-679 2-38, 
which show how very quickly the function approaches the x axis. 

7-8. Lommel integrals. 

The Lommel integrals concerning I can evidently be taken over 
bodily, merely writing K for I. There is, however, one important 
difference. It is no longer permissible to use the origin as the lower 
limit of integration. In compensation, we can use infinity as the upper 
limit. The fact that K n tends to zero as x tends to infinity is one of its 
most valuable assets. 

7-9. Zeros of the modified equation. 

The modified equation has certain peculiarities that are worth 
noticing, though they are rarely mentioned. It has already been 
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mentioned that no solution can have more than one zero, and a reason 
was assigned. If 

y = AI„{x) + BK n {x) 

is a solution, we have 


dx\lj dx \ I„ I xl n 2 

Hence 



Jf a, $ are supposed zeros of y, the left side is zero; but the right side, 
with its positive integrand, cannot integrate to zero, whatever the 
limits. Hence y cannot have more than one zero. 

There is no loss of generality in taking both the arbitrary A, B to 
be positive. It can then be shown that a solution certainly can have 
one zero. We merely adopt the solution y = AI n — BK n which 
ranges from negative infinity at. the left to positive infinity at the 
right, and accordingly crosses OX. This solution is certainly inflected 
and cannot have a stationary value. The proof of both statements 
comes immediately from 

y' = A {'.« +1 a ,} + + n - a :..} . 


The right side is essentially positive and y' cannot be zero, so that. 

there is no stationary value. On the other hand, 1 / ranges from positive 

infinity at the left to positive infinity at the right. It must accordingly 

reach a minimum somewhere between, and a minimum of «' is an 
inflection. 

_I“ P aralle ^ 'rith the above, the alternative form of the solution 
y - AI n -f BK n can have no zero, the function being always positive, 
bince it ranges from infinity at the left to infinity at the right, it must 

ave a minimum somewhere between. Hence it either crosses OX at 
least twice or not at all. Since it cannot cross OX twice, it does not 
cross at all. From the differential coefficient we have 


riJht ra Tff\ fr °? n T tiVG infinity at the left to P° sitive in finity at the 

hi a minirn^rf 01 T haS a Zer °’ Whlch again prOVeS that the funcfc ion 
a minimum. To sum up, a solution must have either a zero or a 
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minimum and cannot have both. No solution can touch OX except 
at the origin, for this would be equivalent to two positive zeros, and 
anyway the matter was disposed of in the theoretical work of Chap. II. 


EXERCISES 

1. Establish the following results: 

(») Yx 

(ii) ^{x-"A' n (x)} = -x-"K n+l (x). 

(iii) K n+1 I n + K n I n+1 = 


(iv) J {K n+l * - K„_ 1 *)xdx = 2nK n * 


(*)• 


2. Prove that 


In = Ai+2 + 


2 n + 1 


n+1 


+ 


w(n — 1) 


n* 


2n - 1 n(n + 1) 

X n " = K n -2 + ~ A n _, + ^ K n . 

Deduce that, for positive integral orders, / n " is positive and steadily increases, 
whilst K n " is positive and steadily decreases. 

3. Deduce from the last exercise that the solution which has a zero has an 
inflection; the solution which has no zero has no inflection. 

4. If a solution has more than one stationary value, it must have an odd 
number of them; prove that it cannot have more than one. 

5. If c, n and m are positive numbers, prove that the equation 

/„(x) = cK m (x) 

has only one solution. 

6. Verify that 

xi{AI n+l (cx ) + BK n+i (cx)} 

is the solution of 


g - ['•+ 


7. Establish the following results: 


(i) £{/o( cx ‘)}= *«-*/,(«!). 

(ii> Tx (• K ’»( cx ‘ ) ) = -l 6 *"**!* 6 * 1 )- 

(iii) ^ {x»/,(cx»)} = ic/ 0 (cxi). 


(iv) {xlA'^cx*)} = — £cA” 0 (cx*). 
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710. Motion of an augmenting mass. 

Having disposed of the necessary theoretical preliminaries, we can 
now turn our attention to some applications of the modified Bessel 
functions. We begin with an example taken from dynamics. 


Problem 11.—A moving body suffers accrefio>i , so that its mass M at 
time t is m (c -f- 0/ c - B 15 repelled from the origin, the force per unit 
mass being proportional to the distance. Initially the body is ejected 
from the origin with velocity u. Discuss the subsequent motion. 

Equating the force to the rate of change of linear momentum, we 
have the equation of motion 

i, (" *) - 


where P is the constant of proportionality. Since M = m(c t)jc 
we have by logarithmic differentiation 

L ( IM = 1 

M dt c -f- t 


The equation of motion becomes 

(Px 1 dx 
dt- c -f- t hi 



Choosing a new independent variable defined by 

c + l = r, dt = dr, 

we have finally pL + - p - k*z = 0. 

dr* r rfr 


This is the modified equation in almost its simplest form, and it has the 
solution 

x = AI 0 (fcr) + BKjJkr). 

From this we derive the velocity 

J t = % = k ( A, ' {kr ^ ~ Sff .(*T)}. 


It remains to determine the arbitrary constants from the initial 
cations 


con- 



dx 




114 


APPLIED BESSEL FUNCTIONS 


We therefore have 

AIq(1cc) + BK 0 (kc), 

AI x {kc) — BK y (kc). 

—IB _ 1 

uK 0 ul 0 I^q+IqKi 

In virtue of the relation 7-9, Ex. 1 (iii) 

I^Koikc) -f l Q (kc)K x (kc) = 

we derive 

A = cuK 0 (kc), B = — cul 0 (kc ), 

and the solution is completely known. It remains to examine some of 
its consequences. Seeing that the coefficients A, B are of opposite sign, 
we can at once say that the solution has a zero and an inflection, but 
no minimum. The dynamical interpretation of this geometrical lan¬ 
guage is that there is an instant of zero displacement and an instant of 
zero acceleration, but no instant of zero velocity or temporary rest. 
We proceed to inspect the matter a little more closely. The full ex¬ 
pression for the displacement is 

x = cu{K 0 (kc)I 0 {kT) — I 0 (kc)K 0 (kr)}. 

Regarding the right side as a function of r, it has the obvious zero 
t = c, which is simply the initial condition of motion. For the velocity 
we have 

d 4 = kcu{K 0 (kc)I x (kr) + I^kcjK^kr)}. 
at 

This certainly has no zero since everything on the right is essentially 
positive; the body accordingly has no rest position. The acceleration 
will be zero if AI X = BK X , or 

K 0 (kc)I^(kr) = Ukc) {fi- 0 (i-T) + i K^kr) J. 

The existence of a positive root for this equation is easily demon¬ 
strated. Differentiation of the series for I x shows that the left side 
increases steadily from a non-zero constant to infinity. Both terms 
in the bracket on the right steadily decrease from infinity to zero. 
Consequently the two graphs must cross just once and there is a unique 


0 = 
u _ 

k ~ 

These solve to 

kA 
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value for r when the velocity is minimum. The body loses velocity on 
leaving the origin, but afterwards it accelerates continuously. It is 
something of a paradox that the acceleration can temporarily vanish 
despite the unremitting activity of the force. 


EXERCISES 

(The notation of the text is used throughout) 

1. If the body is initially at rest at distance h from the origin, prove that 
its distance x at any subsequent time is given by 

= *i(*«)fo(*T) + /,(*c)* 0 (*t). 

2. If the body is projected towards the origin with velocity u from a distance 
h, prove that it fails to reach the origin if 

uK 0 (kc) < khK x (kc). 

In this case its nearest approach to the origin occurs at a time given by 

khl^kc) + ul 0 (kc) 

K t (kx) khK^kc) — uK 0 (kc)' 

3. If the body is ejected from the origin with velocity u and suffere attrition 
in accordance with the law 


M 


= ’"(rib)' 


prove that its distance at any subsequent time is given by 

Deduce that the velocity increases continuously. 

Prove that the momentum at any inBtant is 

mcku{K,{ke)IJJn) + I l (kc)K 0 (kx)). 

r: “ nditi °“ »" d 6 - d «» — 

7-11. The variable transmission line. 

Problem 12 .-Discuss the electrical transmission along a variable line. 
tr»n^ 0meWha o diflerent iUustrati “ « afforded by the theory of line 
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point, not necessarily in the line. If 7 be the current and F be the 
voltage across the loop at A, we have from elementary considerations 
the equations 


-™ = RI + L 

dx 


dl . 
dt' 




= GV + C 


dv 

dt 


In the simple case of a leaky telegraph line it is customary to put 
L =.0 = C. This gives 


d 2 V 
dx 2 


= RGV, 


which is soluble by hyperbolic functions when RG is constant. 
In the less simple case we have 


^ = CL~ + (CR + GL) ~ + RGV. 

dx 2 dt 2 K dt 


In virtue of Fourier analysis one need consider only sinusoidal varia¬ 
tions. We accordingly transform to rotating vector equations on 
replacing V by Fexp(iW). This gives 

~ = (R + iojL)(G + i<oC)V, 


whose solution is by analogy expressible in complex hyperbolic func¬ 
tions of Px , where P is the complex propagation constant defined by 

P = {(7? + iwL)(G + iwC )}*. 


The foregoing can be seen in E. Mallett, Telegraphy and Telephony, 
and in Kennelly’s work and tables of the complex hyperbolic func¬ 
tions. 

It is accordingly evident that if we adopt a series impedance Z = 
R + iioL, and a shunt admittance Y = G + iatC, the analysis can 
be written 




It is customary in practice to boost the inductance by loading the line, 
and for technical reasons that need not concern us the loading is tapered. 
Even the capacity may vary appreciably in a single overhead line with 
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a pronounced sag. We accordingly abandon the constancy of the 
coefficients and derive 


d_ 1 dV 
dx Z dx 




d 2 V 
dx 2 


'""-rzr-o. 

Z dr dx 


It is now apparent that, formally at. least, the problem is closely 
allied to some of our previous work, notably the one-dimensional heat 
conduction in a heterogeneous medium. Among the various possi¬ 
bilities we may take 


These give 


Z = Z 0 r a , Y = Y 0 r\ 


Y 0 Z o 



<PV_adV 

dx 2 x dx 


k 2 Yx^ b = 0 . 


Comparison with our general standard 7-3(1) gives 

1 — 2a = — a, py = k, 2(y — 1) = a -f- 6, 

From these we derive 

a = *(l-f-a), /? = y = 1 + l(« + b), 


a 

n = ~. 
V 


1 + a 


2 -f- q -[- b 


It appears that the case <z -f- 6 = —2 is not soluble by Bessel functions; 
it leads to a well-known type of equation that is soluble by exponen¬ 
tials. Apart from this exception we have the solution 


V = x*{AI n (pxy) + BK n (pxy)}. 

Various writers have investigated simple cases. In general if a■ = b 

the functions are of order £. A particular case of this is a = 0 = 6. 

The line is then uniform and we revert to the simple case of solution 

by hyperbolic functions, which accords with the known values of 
/j and K 

The case -a = 6 was called by Ballantine the Bessel line, and the 
particular case -a = 1 = 6, which has the solution 


V = AI 0 (kx) + BK 0 (kx), 

!f*“ 0 . Wn SS th ° H ? aviside Bessel ,ine ‘ If « = 1, b = 0 the scries imped- 
nce tapers and the shunt admittance is constant. The reverse case 
nas also been investigated, as has the case w-here « ^ -1 b = Q 
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A modification, parallel to some of our previous work, can be made 
by introducing a fictitious origin otherwise by writing Z = Z 0 (x -f p)\p. 
Here p is a complex number such that ZJp is the constant complex 
change of series impedance per unit length. The essential difference 
between this and the previous case is that there is now no correspond¬ 
ing real origin for the series impedance. The equation for V becomes 


d 2 V 

~(h* 


1 dV 

p -f- x dx 


- rz 0 i 


r p + X 




The substitution 


produces 


d 2 V 

dz 2 


p -f x = z, dx = dz, 

\ d J.-khV = 0, lc 2 = ^ 

z dz p 


on the assumption that Y is constant. Comparison with our standard 
7-3(1) gives 

1 — 2 a = —1, py = fc, 2(y — 1) = 1, ny = a, 

whence 

a — 1 , /5 = y = U, n = 5 . 

The solution is therefore 

V = z{AI i (lkz 312 ) + BKtfkz”)}, 

and the current can be derived from 

T __ldV 
Z dz * 

It must be borne in mind that the above is little more than a formal 
solution. Any attempt to utilize the formulae for purposes of computa¬ 
tion must take cognizance of the question of phase arising from the 
complexity of the constant /.-. This technicality lies outside our scope. 

Theoretically there is no reason why both the series impedance 
and shunt admittance should not both vary in accordance with a law 
of the type {x -f- p) a , as was done with the density and specific heat in 
the case of one-dimensional heat conduction in a heterogeneous medium. 
From our point of view these can lead to nothing particularly new or 
interesting, and possibly they are not of sufficient practical importance 
to warrant investigation. Those interested may like to consult a paper 
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by A. T. Starr on the non-uniform transmission line which appeared 
on p. 1052 of Vol. 20 of the Proceedings of the Institute of Radio En¬ 
gineers in 1932. 


7-12. Equilibrium of a non-uniform tie-bar. 

It is well known that in the more elementary theory of beams and 
struts, which usually serves to exemplify differential equations with 
constant coefficients, the problem is generally soluble in terms of 
trigonometrical functions. It is very rarely pointed out that if the 
strut becomes a tie and is laterally loaded, the solution involves 
hyperbolic functions. This being so, it is hardly surprising that certain 
departures of the tie from uniformity lead to a solution by the modified 
Bessel functions. One would expect any such solution to revert to 
hyperbolic functions when the departure from uniformity is made 
indefinitely small. We propose to discuss a problem of this type. 


Problem 13. A horizontal beam of length 2A has uniform load w per unit 
length. It is simply supported at the ends and there is end-pull P. 
The moment of inertia of section, I, at distance xfrom the left end is given 
ty 1 = B(x -f h)/h. Discuss the equilibrium. 

It wall both illustrate the method and serve as a comparison if we 
briefly sketch the solution when the tie is uniform. The y axis is taken 
positive downwards since the deflection is in that direction; the origin 
can be conveniently taken in the middle since the arrangement"is 
symmetrical. The vertical reactions at the ends are each wX, and 

taking moments about any point in the right half-span we have as the 
equation of equilibrium 

Ely" = py -f- $w{x -f A) 2 — wX(x + A). 

Tins evidently has a polynomial particular integral of quadratic type 
We accordingly substitute 

y = a + p(x -f- A) + y (x + A) 2 . 

On comparing coefficients we derive 


w 

tfP’ 





P 

El' 


The full solution can be written 

y = A coshn* + Bsinhwa: + a + p(x + A) + y(x + A) 2 . 
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The two end conditions being y = 0, x = A, —A we have 

0 = A cosh A + B sinhnA + a, 

0 = A coshnA — B sinhnA -f a. 


From these we derive the arbitrary constants as B = 0, A — 
— asechnA. The solution is accordingly 



w I cosh nx 
Pn 2 IcoshnA 


w 


+ * 2 >- 


Turning now to the problem in the proposed form, we have as the 
equation of equilibrium 

^ (x + h + X)y" = Py + iw(x + A) 2 - wX(x + A). 
h 

There is again a polynomial particular integral of quadratic type, and 
carrying out the same procedure as before we derive 


w 


wX 

W 

Pn 2 ’ 

T ' 

P/m 2 ’ 

w 

4 Ph 

= V 2 

= 4/m 2 . 

2 P' 

eh 




The complementary function is derived from 

d 2 y Ph y _ 

dx 2 EH x -j- k + A 

The substitution x + £ + A = z suggests itself, and if we override 
the objection to having two forms of the independent variable in the 
same expression, we have the solution 7*3(1), 

y = z'{AI y (vz') + BK^vz')} + a + /3(x + A) + y(x + A) 2 . 

With a slight reduction on using l = h + 2A, the end conditions left 
and right respectively give 

0 = h>{AI^h*) + BK^h')} - 

w Z 

0 = l^AI^vl') + BK^vl*)} 
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for the determination of the arbitrary constants. They solve to 

iv/Pi, 2 

h{Ilivk^K^vP) - I x ( v P)K x ( v h')} 

_ A 

h'K x (xP) - l'K x (vh') 

_ B 


It is a matter of elementary algebra, which can be left to the reader, 
to show that the particular integral is expressible as 


w 

PiT 2 


*-± 1+1 + £ ( A , _ *>), 


so that the similarity with the previous solution is quite marked. If 
we allow h to become indefinitely large the beam becomes uniform in 
section and the variform solution should go over into the uniform 
solution. The term (x + h -f- A )/h evidently goes over to unity, so that 
three of the four terras in the uniform solution are reproduced. The 
proof for the remaining term requires more than average manipulative 
skill and we shall not pursue it here. Instead, we offer a simpler illus¬ 
tration of the method. 


Problem 14. An investigation on a variform member leads to the result 


where 



?,nif rt kn °f Wn - that 1 . aS 1 k becomes indefi nitely large the member tends to 
uniformity, in which case the solution is known to be 


^ = ^2 [1 — sech n A]. 
It is required to reconcile the results. 


We have 


l =A + A = i( 1 + *), 

" + *'(■+£> 


v(l* - A*) 



(O 160 ) 
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Testing our belief in our unproven asymptotic formula? 7-4, Ex. 15, and 
7-6(8) 


we have 


Similarly 


Hence 


- &*■ 


-> ~ exp - A')} = 


Uvh')KA*i') -* exp 2 ( ^" A) . 

}-►/<* cosh n A, 


and the concordance is established. 

These problems on variable bars furnish physical evidence that the 
asymptotic formulae are independent of the order. By putting 



we can, by a proper choice of r, make the resulting functions to be of 
anv order we like. Whatever the order may be. the bar becomes uni¬ 
form as h increases indefinitely. And there is only one solution of a 
particular problem for a uniform bar. 


EXERCISES 

1. In the problem of the uniformly loaded beam, prove that of the three 
denominators under A, B, &c., none can be zero. This proves that no adjustment 
of the constants will make the solution depend on the one typo of function with¬ 
out the other. This is not usually true of solutions that depend on J n and I 

Prove further that the denominator under B is essentially negative. 

2. A cantilever of length X has end load W and end pull P. The moment of 
inertia of section being given by / = H(x + h)/h, prove that, with the notation 
of the text, the retaining couple C at the wall is given by 

vC{/ 1 (v/*)A’ 0 (v*») + / 0 (v/>‘)A’ 1 (v/‘)} 

= 2IFA*{/ I (v/*)ff 1 (vM) - I^h^K^U)}. 

It is known that, when the section is uniform, the corresponding solution is C = 
(IF/n) tanlinX. Reconcile the two results, the latter being tho limiting form 

h increases indefinitely. 

3. A cantilever of length X has uniform load m per unit length and end pull 
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P . The retaining couple is C and the moment of inertia of section at distance x 
from the wall is I — H{z + h)/h. Verify the following details of the analysis. 
Tho equilibrium equation i6 

Ell if = C + Py + iwx 2 - ,rX*. 

n 

The assumption of a particular integral of the form y = a + fix + yx 2 leads to 



Tho complementary function is 

y = z^/^vz*) + v 2 = 4 n 2 A, z = x + h, 

and the end conditions at the left lead to 


-A = 2pA*/f,(vAI) + av/f 0 (vAl), 
B= 2pAi/,(vAi) — otv/^v/il), 

The right end being free from applied couple, we deduce 


where 


Cn 2 _ [ ] 

w 2 nh{ } 


- 1, l=h+\. 



+ 2(1 - - Xi(vAt)/,(vfi)], 


{ } = + / 0 (vA»)A-,(v/*). 


It can be shown that, in the case of a uniform beam, the corresponding solution 
is 

Cn 2 _ 1 + nX sinhnX _ ^ 

w coshnX 



CHAPTER VIII 


Applications to 
Hydrodynamics and Elasticity 

81. Tidal motion. 

The problems so far discussed in illustration of the theory have 
been of a relatively simple kind, involving but little analysis. We 
come now to a number of problems of somewhat more complicated 
type, beginning with surface waves on a liquid. These are distinct 
from waves of expansion which are propagated throughout the body 
of a medium, like waves of sound or light; and a further distinction 
is drawn according as the liquid is relatively shallow or deep. We 
begin with the former. 

In a current of liquid, the motion of the particles is secular and 
the liquid continuously moves forward with the current. The contrast 
in wave motion is the presumption that a particle makes only a limited 
excursion from a mean position, its displacement being a compound 
of vertical and horizontal oscillations. 

It is first necessary to set up two fundamental equations. One of 
these is a physical equation, known as the equation of continuity and 
based on the incompressibility of the liquid. It states that the mass of 
liquid in any specified space-element is fixed by the volume of the 
element and is invariable. The second equation is purely dynamical. 

Consider a horizontal canal running in the direction OX. Any 
surface phenomenon will be presumed to be uniform across the breadth 
of the canal, so we are working with two dimensions in a vertical plane. 
The x axis is in the surface of the liquid in its rest position and we take 
two fixed vertical planes L, M defined by x and x -f- Bx. The surface 
breadth in L being 6 , we take h to be the mean depth of the possibly 
non-uniform cross-section A, so that A =lh. 

When the surface is agitated by a wave we require two moving 
co-ordinates. We take 77 to be the small surface elevation above the 
rest position, uniform across the breadth in plane L\ £ to be the hori¬ 
zontal forward displacement of any particle in the plane L, so that 
the velocity d£/dt is uniform over the plane L and extends from the 

124 
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surface to the canal bed. The justification for this assumption appears 
in the consistency of the results. 

The rate of liquid entry into the lamina LM from the left is Q = 
bhdg/dt. The corresponding exit across M is Q -f- 8(?. The accession 
—8Q in the lamellar space between L and M causes a surface rise 
dg/dt over an area measuring b by 8x. This gives 



equivalent to 


h + 



m) = 0, 


as the equation of continuity. 

Coming now to the dynamical equation, consider an clement of 
length 8x and small cross-section a, so that its mass is paSx and moves 
with an acceleration d 2 £/ct 2 . This acceleration is caused by the excess 
pressure at the one end as compared with the other, and this in turn 
is due to the difference of head at the planes L and M. We thus have 



We can now eliminate £ and derive 



The particular case of a uniform canal, with b and h both constant is 
not our immediate concern. It evidently leads to 

gh d -h = d - 2 

y dx 2 dt n -' 


and the solution is known to be waves travelling with velocity (gh)* 
ca n nl OSS f 7 g n ng Standing Wavcs - Leavi "g this and coming to 

- rr: & vs ° f 

s{ M S} + y i = °- 

Much of hydrodynamics is so difficult that it is common practice first 
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to solve the necessary equations and afterwards seek the problem 
which the solution fits. One of the simplest hypotheses that one can 
make in the above equation is to treat h as constant and take 6 pro¬ 
portional to x. This gives a crude representation of a V-shaped estuary, 
of uniform depth, that fans out from the origin. The equation becomes 


s+is+^-«. 


CO 


/l 


The solution which gives a finite elevation at the origin is rj = BJ 0 (kx), 
and if we suppose that the distance to the mouth is A where the open 
sea has a tidal rise given by 77 = C cos (cut -{- e), the solution 


1 


° C0S ( wt + 

«/()(* *) 



fits at both ends. In view of the known steady decline in the absolute 
value of the maxima and minima as we move to the right, we get an 
intimation of the increase in wave-height as the tide comes up the 
channel. 


8-2. Canal of finite length. 

The above problem is deceptively simple and there are very few 
like it. We turn instead to the problem of the canal of finite length. 
Reverting to the equation for the uniform canal, we have 

c = 

dx 2 ?l 2 ’ 



showing that the velocity of propagation is “ due to a fall through 
half the depth ”. An infinity of solutions can be built up of the form 

. sin sin 
£ — mx met, 
cos cos 

but if the canal is closed by two vertical barriers at x = 0 and x = A 
the fluid can have no horizontal velocity there. We accordingly have 
as boundary conditions 

|f = 0, x = 0, A. 
ct 

A suitable solution is then $ = B sinfwx sin (met -f- e), provided 

sinmA = 0, mX = it, 2i r, 37 t, ... . 

This gives the possible frequencies of the normal modes of vibration, 
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where every surface-particle goes through its equilibrium position at 
the same instant. 

It is natural to inquire what is the analogue for a non-uniform 
canal. If a simple solution is desired, trouble arises from the equation 
of continuity. Suppose we solve the equation of motion for 77 ; then 
an integration is needed to find £, which must be known in order to 
satisfy the boundary conditions. Very few expressions involving 
Bessel functions arc readily integrable, a difficulty similar to that 
encountered when attempting to evaluate Fourier coefficients. Alter¬ 
natively, if we solve for we have to differentiate to find 17 and thereby 
lose the simplicity of the wave form. 

Taking the equation for 77 we assume with some generality that 
bh oc xP and b oc xP. The equation becomes 


<Pr, 

dx 2 


x dx g 


ox*->’ = 0, 


where a is a constant. 

This, of course, is readily soluble; but our chance of finding an infer¬ 
rable result is based on the relation 


d 

dt 


{C* n+li J n +i wn} = Pyt«"+*>-V „(#>). 


Comparison with our standard 3-6(2) gives 


If 


l — 2a=p, 2y-2 = q — p, ± » = * 

Y 


V = x a J n (fix*), 677 oc x a +ij„(px'), 
the integrabilifcy of 

bh£ = — jbt]dx 

is ensured by taking 

° + ? = y(» + 2) — 1, q = o y _i + _ x> 

P = 1, a = 0 = n. 

As an illustration we may take ? = 2 . The problem and its solution 
then reads: If the variable breadth, depth and cross-section are given 


b = B^l±l 


)■ *-s 


/ 


Bf+z 


_ jf+x 


bh — A 


f ’ 
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the equation for r/ is 


IKr'SK* I'-fT*-”- 


The substitution z = f -f- x gives 


lUtt) 

dz \ dzj 


+ t ££ z2 v = °> 

9fA 


d* 1 ) i 1 dij 9 ~ _ 4 (o 2 B 

d? + z7z + 4^-°’ ^ = 9^T 


The solution is verified in the usual wav to be 

* 

7) = LJ 0 (Pz 312 ) + MY 0 (pz* n ), 

whence we derive 

Af* t _ r.o , 7 _ 


-f *-!+* 


= £-z*»{LJ 1 (f}n + MY 1 (Pz i % 

dp 


If the ends of the canal are defined by x = 0, z =/ and x = A, 
z =/ -f- A = both ends are velocity nodes if 

0 = U x W 3n ) + MYtff). 

0 = LJ i{fil 3r> ) + MY x (pl 317 ). 


The elimination of the constants gives 

Jxifjp _ Ji(pn 

Y l (pr») Y x (pn' 

This transcendental equation determines ft and thence oj and the 
normal modes of vibration. 

Translated into figures, we find from the tables that J X /Y x = 0-954 
when the argument is 12-56 and again when the argument is 15-71. 
Thus 

pl 3r> _ 15-71 
Pf 3p ~ 12-56’ 

whence we derive Iff = 1*161 and / = 6-2A. This means that the 
fictitious origin is distant rather more than six times the length of the 
canal; and the variation in mean depth and cross-section is about 
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sixteen per rent as between one end and the other. The period is 
calculated from 

p., _ 4 cu 2 B _ (12-o(i )- 

P ' 9 jjfA 

Hence fora model tank 1 metre long, the water being 10cm. deep at the 
deep end, we have c u = 9-48 and the period 2 ~ /" = OG6 sees. 

The position of the surface-node is given by yj = U, so that 

Wz”): Y 0 (ftz 3r -) = 0-954. 


Interpolating from the tables we have fiz in = 
we get 



cm. 


14-17. 


As ftf * = 12-56 


This is about what one might expect. 

A good deal more latitude of choice is attained on approaching the 
problem through £. The elimination of tj from the dynamical equation 
and the equation of continuity gives 


d 

dx 






If we make the same tentative assumptions regarding the breadth and 
area, we have, where a is a constant, 


' 4- 2 JL~l3 d A . 

da? x dx + 


f p(P ~ 7~ 1) 
l x 2 



OJ 2 

— or 1 -* 
9 




Provided p - q * 2 this is soluble in terms of Bessel functions and 
comparison with the standard gives 3-6(2) 

l-2a = 2^- ? , a*-n*y* = p(p- q - 1), 2y - 2 = q - p. 

Since 


bhg oc (px y ) 


* C0m pa™tively manageable form is obtained by matin" a + „ 
± ny. We then have ° ry 


n V = (a + 7>) 2 , 

»V-a'= 2 f a+ p !=- f ( J ,_ 5 _ 1)) 

P( 2a + 2p — q — 1) = o. 


Gtdm&ri? Iveseuich Observatory 
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It follows that either p is zero or else the first comparison equation is 
satisfied automatically. The former alternative gives 


a = ^(1 + q), 7 = 1 + \q, 



l±q 

2 + ? 


As a particular case of this we can take q to he also zero. The solution 
is then in terms of x*«/j(x), which is trigonometrical and correct since 
the canal is then uniform. 

The second alternative gives 


a = 1(1 + q) — p, y = 1 — l(p — q), 



1+g 
2 — p + q 


This opens up endless possibilities and we may illustrate by taking 
p = 3, q = 2. The problem and its solution then read: If the variable 
breadth, depth and cross-section are given by 



the equation for £ is 

d r p d ( ( f+ X )• ,n , - 2 e 

dx -+ x ) 2 dx \ P ') J g 

or. 


= 0, 


so that 


The solution is 


£ = + MYJftf)}, 


and correspondingly 

V = ~ 


PA 

2 fB 


z-'{U 2 (Pz>) + MY 2 (Pz »)}. 


As an illustration of the difficulties likely to be encountered in 
translating this into figures, suppose the mean depth shows a 5 per 
cent variation as between the two ends. Then the breadth varies about 
10 per cent and the cross-section varies about 15 per cent. This may 
legitimately be described as gradual. Suppose c is a zero of the cylinder 
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function that defines £; then the next consecutive zero is about 
c -f 77 . The end conditions therefore give 

pP = c_+jr 

Pp c 

Since we are assuming that ///= 21/20 we have (///)* = 41/40 
approximately, and c is about 40 tt. This would exhaust the possi¬ 
bilities of any known tables. 

It is not impossible for the solution to come out in terms of the one 
type of function without the other; compare 7*12, Ex. 1. For example, 
the thirty-fifth zero of Y 3 is 112*273, and if this happens to be the 
value of pf* the first equation of condition would give L — 0. The 
next zero would give pP — 115*416, whence we conclude that / is 
about 17-7A. If A is known, so is /?; the value of to then depends on 
the adopted value of A/B. 

8*3. Surface waves. 

The type of motion hitherto considered, where the horizontal dis¬ 
turbance is presumed uniform throughout the depth of the liquid, and 
the vertical motion of a particle is more or less ignored, is known as 
tidal. We propose now’ to take a less restrictive hypothesis and to 
W’ork in cylindrical co-ordinates. This will have the advantage of 
introducing Laplace’s equation, which plays so large a part in classical 
mathematical physics. 

Let u, v, w be the component velocities of a particle in the direc¬ 
tions of increasing r, 6, z respectively and consider a volume element 
defined by the coaxial cylinders of radii r, r -f- 8 r; by horizontal 
planes at levels z, z-{- 8 z; by axial planes of azimuthal angles 6, 
0 -f- 80. 

The rate of liquid influx across the inner vertical face is 

Q l = ur898z. 

The corresponding efflux across the opposite face is + 8 Q„ and 
the resulting efflux from the volume element is 

SQi = |-(wr808z)Sr. 

cr 

Similarly from the top and bottom faces we have 

SQ 2 — ( wr8$8r)8z. 

oz 
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For the remaining pair of vertical faces we have 


^ (y Sr 82 ) 50. 


The incompressibility of the liquid thus gives the equation of 
tinuity 


a 

dr 



con- 


It is common practice in dynamics and electrostatics to adopt an 
energy potential whose negative gradient gives the force in any par¬ 
ticular direction. A similar idea is adopted in hydrodynamics, a velo¬ 
city potential being assumed whose negative gradient gives the velocity 
in any direction. It corresponds to a type of motion known as irrota- 
tional, in contradistinction to rotational or vortex motion. 

Denoting this velocity potential by 6, we have by definition 




d<f> 

dz 


The equation of continuity becomes 


dr 




and this is equivalent to Laplace’s equation in cylindrical co-ordinates. 

The oidy satisfactory way of attempting to solve a partial differ¬ 
ential equation of this type is to assume initially that <f> has the form 
FH, where F is a function of r and z, whilst // is a function of 0 alone. 
The equation can then be written 



+ r* —= 0. 
II dd 2 dz 2 


Since 0 is a variable independent of r and z, this implies that the middle 
term must be a constant, and with a view to Fourier analysis we put 

I MI 2 

-= — n 2 , 

II dd 


so that II OC cos(n0 -f- e), with n an integer. 
The equation is now 
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There is already an adumbration of Bessels equation, and if we put 
F = RZ, where R is a function of r alone, and Z a function of z alone, 
we have 


1/1 d / n 2 1 1 d-Z 

R\rdr\dr) r- j Z dz- 



As before, we conclude that the last term must be a constant, and in 
our present case it is convenient to take Z oc cosh k(z + h). This leaves 


with the solution 



R = AJ n {Kr) 4- BY n ( K r). 


We now assume that the horizontal plane defined by z = 0 is the 
surface of still liquid of uniform depth h. When the water is mildly 

agitated, a particle at the bottom can have no vertical velocity, so 
that 



when 



This condition is already fulfilled by our choice of Z. 

There is nothing in the analysis so far that specifically connects 
it with wave motion. The connexion is made by introducing a periodic 
time factor into the velocity potential, which thus becomes 

4> = {AJ n ( K r) -f BY n { K r)} cos ( 7 id -f- c ) cosh k(z + h) cos mt. 


Ifji denote the slight surface elevation at any point, its vertical velo¬ 
city is given by 

d<f> . 

1 = — ~ when 2 = 0. 
dz 


As this j s proportional to cos mt, integration gives v oc sin ml, so that 

movinfth wt‘ en 1 W u' Cn v is Zer ° and the surface-particles are 
moving through their equilibrium positions. 

on wZ? 1118 t0 C ° n ? eCt the Various Parameters. It is shown in texts 
on hydrodynamics that m the absence of viscosity and impressed 

forces other than gravity, the plane 2 = 0 is a free surface prodded 


d 2 tf> dd> 


A 
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This evidently leads to 


gx tanli xh = m 2 , 


where incidentally tanli xh is necessarily less than unity and is asymp¬ 
totic to it for large values of xh. 

If the liquid is confined between two retaining walls defined by 
r = a, b, the horizontal velocity u must be zero at each of these, so 
that 


dS A 

U = — -L = 0 

cr 


when r = a, b. 


The arbitrary constants A, B can then be eliminated and we have 

J n '(xa) _ J n '(xb ) 

Y n '(xa) Y n '(xb) 


as the equation to determine k. 

For the sake of arriving at figures we suppose firstly that there is 
no inner bastion. We then discard Y n as giving an infinity at the 
origin. Next, suppose the water is divided quadrantally by two vertical 
and perpendicular barriers. At each of these we need no normal velo¬ 
city, so that 



This is achieved by dropping the phase angle e and putting n — 2. 
We now have 

cf> = AJ 2 (xr) cos 2 9 cosh#c(z + h) cos mt. 

If the outer boundary is defined by r = a, the radial velocity must 
there be zero, and 

u = - d + = 0, J 2 '(m)= 0. 
dr 

The recurrence relation 

2J„ = J „_i J v* i 

then tells us that the required zeros are to be found from 

J^xa) = J 3 (xa). . 

To find, say, the first root of this equation, J 3 is the slower m taking 
off from the axis and its first zero is beyond the first zero of J j, sq 
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the required root is rather less than the first zero of J v From the 
extracted values 

x J l 

30 0-33905 0-30906 

3-1 0-30092 0-32641 

we deduce that the required zero is 3-054 = ko. 

For a tank whose depth equals its diameter we have kTx = 6-1 and 
the distinction between tanli kIi and unity is not worth making. If the 
diameter is a metre, we have 

m = {gK) k = 7-72, 
and the period of oscillation is 

— = 0-81 sec. 
m 


EXERCISES 

1. Discuss the tidal waves in a channel of uniform breadth, the depth shelving 

uniformly from zero down to the open sea where there is a tidal rise rj = 

Ceos (cot -f- e). Prove that the function is of order zero and the argument pro¬ 
portional to xi. 

2. Discuss the normal modes of vibration for surfaco waves in a semicircular 
tank of radius a and depth h. Compute from the values 


X 

Jo 

Jo 

1-8 

0-3400 

0-3061 

1-9 

0-2818 

0 3300 


8-4. Heat conduction. 

The analysis for heat conduction in cylindrical co-ordinates is so 
close to our previous work as almost to amount to a transcript. Wc 
merely call <f> the temperature instead of the velocity potential, and 
pay attention to a few necessary coefficients. Taking the same cylin¬ 
drical space element as before, we let H 1 be the heat per second that 
flows across the inner vertical face. Then H t + 8H, flows across the 
opposite face and the accumulation in the volume element is ~8H , 
The amount H x depends jointly on the area which it crosses, the con- 
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ductivity k , and the negative temperature gradient radially. We thus 
have 


H 1 = 


— ^ xrhdhz, 
dr 



a//, 

dr 


hr = +* 1. 

dr 



Similarly if //., flows upwards through the lower face we have 


-5//, =dz = Kr d ^t- hrhdhz. 

cz dz 2 


Laterally we have 


_W. = »• M - = MM* 

d6 rdff- 


The total heat-gain shows itself in a rise of temperature. Taking density 
p and specific heat s, we have the heat gain per second as 

ps ^ rhrhdhz. 
r dt 


The equation for <f> is therefore 

m 1 5 i d 2 <ft I 1 d 2 <f> _ ps dcf> 

r dr\ dr f dz 2 r 2 d6 2 k dt 

The theory of heat conduction is an extensive study in itself and 
as this is not the place to pursue it, all we can do is to indicate how 
Bessel functions occasionally play a part in the analysis. When the 
temperature has ceased to vary with the time, we have what is known 
as the “ steady state ” and the right side of the last equation falls out, 
thus simplifying the working. 


Problem 15 .—Hot gases are conveyed through a long straight pipe which 
is indifferently lagged; it is required to find the temperature distri¬ 
bution throughout the lagging, V'hose internal and external radii are 
d and c respectively. 

A definite length is under observation and a sufficient number of 
thermo-couples enables the steady temperature distribution at the 
outer surface to be expressed as a Fourier series. It can be assumed 
that the temperature is constant round the circular cross-section of 
any coaxial cylindrical surface. The temperature of the surrounding 
medium can be assumed uniform and will serve as the origin of the 
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temperature scale. It is assumed that tlie heat-loss into the surround¬ 
ing medium is to be accounted for by a constant coefficient of emis- 
sivity h. 

The particular type of Fourier series adopted is a matter for careful 
discrimination since it is very easy unwittingly to introduce extraneous 
hypotheses which are not justifiable. The mean temperature of the 
outer surface is hardly likely to coincide exactly with that of the sur¬ 
rounding medium, so there is pretty certainly a constant term present. 

We accordingly suppose that the temperature of the outer surface is 
given by 

fa = a o + a, cos vz -f- a 2 cos 2 vz -f- . . . 

+ sinvz -f- b 2 sin2nz + . . . . 

The number of terms employed has to be justified by the number of 
thermo-couples; but in any case the coefficients are uniquely known 
from the observational data. 

Under the given assumptions the equation for becomes 

dr- r dr dz 2 


If we suppose the solution is of the form 
function of r alone, we have 


<f> = R exp ioiz, where R is a 


d 2 R . IdR 

M + r Tr ~ W ' R = °' 


This puts us in possession of solutions of the type 

* = {Mo (<*) + BK 0 (wr)} coswz + {C/ 0 (<or) + DK„(wr)} sin wz. 

Sh0 T these ia * h ™-lves are not enough 
he SketofThe^ruZ™ 13 “ ° ft “ h “ - to 


where p and q 

demonstrated, 

solution 


<f> — p -h q logr, 

IZ rff' Th,! existence of ‘his solution is readily 
and by taking vanous values for <o we can build up the 


0 2 {4„/,,(<o„r) ■ B n K 0 (oj„r)} co sco n z 

+ 2{C , , l / 0 (o> n r) -f DnK 0 (a) n r)} sinm „2 
+ V + q logr. 


10 


(O ICO) 
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At the outer boundary, this solution must coincide with (f> e . We con¬ 
clude by comparison, firstly that all the values of to are given by 

oj„ = nv; 

secondly, that the independent terms are related by 

a o = V + 7 

thirdly, that the other coefficients are connected by the relations 

AJ 0 (nvc) + B„K 0 (nvc) = a n , 

CJ 0 {nvc) + D n K 0 (nvc) = b n . 

Each of the last three relations leaves us with effectively one unknown. 

Coming now to the heat emission, we equate the heat emitted by 
any surface element to the heat which it receives. This gives the 
relation 

dd> , , 

—k -L = h<f>, r = c. 
or 


The application of this to each term of our series gives 

— = logc) = ha 0 . 

c 

ha„ 


B^^nvc) — A n I x (nvc) = 
D n K x (nvc) — CJ jCnuc) = 


KHV 

hb n 
" ' • 
KHV 


Taking advantage of the relation 7-9, Ex. 1 (iii) 

I,K X -f- = —» 

twc 

all the constants are now determinable by elementary algebra and the 
solution may be considered complete. 


EXERCISES 

1. Show that tho heat per second that crosses the section defined by z — 0 
can bo expressed as 

— £ — + 2tzk dX{C n I x (nvd) - D n K x (nvd)}. 
v n 

2. Prove that the amount of heat per second that crosses unit area of the 
inner surface averages chajd. 
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8-5. Transverse vibrations of a tapered rod. 

In certain branches of applied mathematics the square of Laplace's 
operator makes its appearance; this is notably the case in elasticity 
and in the motion of a viscous liquid. The fundamental equation is 
then of the fourth order and usually requires more than two different 

types of liessel function for its satisfaction. We begin with a relatively 
simple case. 


Problem 16.—Discuss the transverse vibrations of a tapered rod. 

In the discussion of the transverse vibrations of a taut string, we 
assume the absence of flexural rigidity and we ascribe the vibrations 
to the tension. The essential difference in the parallel discussion of a 
rod is that the flexural rigidity is taken to be paramount, whilst the 
longitudinal stresses arc more or less ignored. 

Imagine a thin rod to occupy part of the x axis and let AB be 8x 
an element of length, the origin for the moment remaining unspecified 
Let« be the cross-section at A; I the moment of inertia of the section*’ 

moment^n i m ^ luS /° r the m »terial; P the density, M the bending 
moment and Q the shearing force. In the approximate theory of bend 

ng the shearing force is the negative gradient of the bending moment 
The latter is given by Euler’s formula and we have * 

= M = El p*, 

dx ar 2 

where y is the deflection. The mass of the element A B x,. i 
accelemion is caused by the element of shearing force SgT This 

"<2.. ‘“‘.o i # _ 

d*y 


-3*—S- 


-£(« 


dx 2 


8x . 


Since oscillatory motion is possible we nut »«_ v • . , 

independent of t. This gives P * ~ * “ P whcrc A * 


IS 


rf2 1 PI** 


} = pauPX. 


x T d * r constants '«“ -l-tion 

valent to «1<* 

various texts on sound, siruetures, elatticity, o^sZ^h of materia,T 
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The most likely departure from uniformity will come from the 
cross-section and its modulus. There is then no standard technique for 
solving the equation; each particular case calls for individual treat¬ 
ment, with the result that considerable ingenuity is often required for 
effecting a solution. If, as an illustration, we assume the rod to be 
uniformly tapered, we can put a = Ax 2 and I = Hx*. The equation 
then becomes 




A poj 2 

~HB' 


At this point, two methods of procedure present themselves; but 
before advancing, it is possibly as well to recall certain features in the 
solution of linear differential equations. 

The equation 



can be written indifferently as 

(D - 2 ){D - 3 )y = 0, (D — 3 )(D - 2)y = 0, D = 


The operators (D — 2) and (D — 3) happen to be permutable. This is 
by no means always the case, as the reader can readily verify from the 
two operators xD and x 2 D. As a result of their permutability, the 
original equation is necessarily satisfied by any solution of either of 
the equations 

(D - 3 )y = 0, (D - 2 )y = 0, 

i.e. by Ae 3x and Be 2x . As only two independent solutions are required 
for a second order equation the given equation is solved. 

Alternatively, certain equations possess a type of homogeneity that 

may be illustrated by 

*?p{ + 3x' l f + 4y = 0. 

dx 2 dx 


The substitution 





leads to an operator $ whose properties may be summarized as 

£. = »(» - 1 )(» - 2 ), = «"(» + «) H - 

dx 3 
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The above equation would then be written 

{«>({> _ 1 ) + 3S> + = 0 

and solved as a linear equation with constant coefficients. 

Applying the former method in the present instance, it is readily 
verified that our equation can be written 

s 

{x 2 D* + 8 xD 2 +12 D-- k*)X = 0. 

The only hope of expressing this in two permutable operators is to 
write 

(xD 2 + aD + k 2 )(xD 2 + aD - k 2 )X = 0. 

It then appears that the two forms are consistent if a = 3 and we have 

(xD 2 + 3 D + k 2 )(xD 2 + 3 D- k 2 )X = 0. 

It follows that a permissible form of X is a solution of either of the 
equations 

F d 2 X 3 dX k 2 

We accordingly have 

X = x~'{AJ£2kx*) + BY 2 (2kx*) + Cl 2 (2kx i ) + DK 2 (2kx*)}. 

Using the second method, which only differs symbolically from the 
former, we should write 

«r*0(& - l){ e »9(S ~1)}X = Pc »X t 

whence 

e' 29 (& + 2)(ff + - l)X = k*X. 

By a slight rearrangement we derive 

+ 2)c^{)({> + 2)X = &X. 


It follows that any solution of either of the two equations 

{e + 2 ) + PpC = 0 = {e~ e &({> + 2 ) — PjAT 

is satisfactory, and it can readily be verified that these are nothing but 
a disguised form of the previous pair, leading to the same result 
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By differentiation we then have 

w 


dX 

dx 


= ~kx~ 3 r ~{AJ z 


By 3 - ci 3 + dk 3 }. 


~ = +Xhr*{AJ i + BY, + Cl , + DK,}, 

the unwritten argument being in every case 2 kx*. The final determina¬ 
tion of and thence a>, is decided by whatever boundary conditions 
we adopt. If, for the sake of illustration, one end is defined by x = 0 
and is free, we discard I’ 4 and K 4 as giving infinity at the origin. If 
the other end is clamped and defined by x = A, so that X = 0 = X', 
we have 

AJ 2 ( 2AA») + C7o(2*A*) = o, 

AJ Z ( 2*A») - C7 3 (2AA*) = 0, 

whence 

J 2 {2kX i ) J 3 (2k\*) _ 

L(2k\*) ^ 7 3 (2/ A*) 

as the transcendental equation for the determination of the normal 
modes of vibration. 


EXERCISES 

1. In the last equation, show that J, and J 3 must have opposite signB and 
deduce that the equation has an infinite number of roots. 

2. If the rod is of circular section, the radius being proportional to z 4/5 , prove 
that the solution involves functions of order 5/6 with arguments proportional 
to x 3/1 . 

3. Investigate the general case where a is proportional to x n and 1 is propor¬ 
tional to x m . Deduce that the problem is soluble when 

3m = n -f 8 or m = n -f 2. 

This does not exhaust the possibilities. 

8 6. Buckling of a circular disc. 

The mathematical theory of elasticity has to be classed as a de¬ 
finitely difficult subject, and' it is not rendered easier by the Cartesian 
method of approach, which happens to be the way it developed. The 
analysis required for the discussion of further problems is too lengthy 
to transcribe here. The reader must accordingly be asked to accept 
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the adopted equations on the authority of the textbooks where he can 
find their justification and master it at leisure, unless he happens to 
be in the fortunate position of being already familiar with it. With 
this reservation we can continue with our applications. 


Problem 17 .—A uniform circular disc, whose middle plane is the z plane, 
has thickness 2 h and radius a. The periphery is subjected to uniform 
radial pressure P per unit area. It is required to investigate the 
critical value of P which causes the disc to buckle. 


The problem is in some respects the two-dimensional analogue of 
the strut but the analysis is more complicated. If we suppose that a 
point whose co-ordinates were r, f, 0 moves to r, f, z, we can discard 
<f> in virtue of the symmetry about the z axis. It is then shown in the 
textbooks that the equilibrium equation has the form 



Adopting the substitution r = e 9 , and writing {> for d/dO as before, the 
equation takes the very simple form 

(e- 2<, $ 2 )(c- 2l, $ 2 + k 2 )z = 0. 

The operators are evidently permutablc, so that any solution of either 
of the equations 

e-v^z = 0, (c-^O 2 + k 2 )z = 0 

will be acceptable. The former gives 

z = A6 -f- B — A Iogr -f- B. 

The latter equation is nothing but Bessel’s equation of order zero, 

<•*»* + = ~ + l ~ ~ + e* = 0 . 

dr 2 r dr 

We are thus in possession of the full solution 

2 = A logr + B + CJ 0 (kr) + DY 0 (kr) 
containing the requisite number of arbitrary constants. Hence also 

~ = At -i - kCJJkr) - kDY^kr). 


If the disc is complete, without a central hole, the values of z and dz fir 
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must not be infinite at the centre. Since Y x goes to infinity logarith¬ 
mically at the origin, we must have both A and D zero, the solution 
reducing to 

z = B+ CJ 0 (kr), d l = -kCJ^kr). 


The arbitrary constants are determined from the boundary conditions, 
there being two standard cases according as the rim is clamped or not. 
We propose to deal with the latter case. 

When the boundary is not clamped there is neither displacement 
nor bending moment at the periphery. Following the textbooks, these 
are covered by 



where a denotes Poisson’s ratio. The former condition is equivalent to 

B -f- CJ 0 {ka) = 0. 

As for the latter condition, the equation for the Bessel function shows 
that it is equivalent to 

kaJ 0 (ka) = (1 — a)J x {ka). 

As Poisson’s ratio is usually about 0-3 and not a subject for meticulous 
accuracy, the equation is quite easy to solve with the aid of the tables. 
It has an infinity of roots; but all except the smallest would correspond 
to highly unstable configurations. The formula for z shows that the 
cross-section of the middle plane looks like the first half-loop of J 0 
and its reflection. The constant B plays no part and might have been 
dropped; it corresponds to a bodily displacement of the disc and has 
no stress effect. The constant C remains indeterminate and gives a 
scale effect. It is of no consequence since if the buckling load is reached 
the damage is usually done. There remains the constant k, a convenient 
abbreviation defined by 


, 3P(1 - a 2 ) 
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EXERCISES 


1. Taking Poisson’s ratio as 0-3, solve the critical equation 

zJ 0 (x) — (l — 

by means of the extracted values 

x J 0 J j 

204 0-2009 0-5738 

2-05 0-1951 0-5730. 


Hence compute for a circular disc of diameter 10 cm. and thickness 1 mm. taking 
Young’s modulus as 2 10® gr./cm a . 


2. In the case of a clamped disc the boundary condition is evidently 

Prove that this leads to J 1 (Aa) = 0 whose smallest root is 3-8317. Hence 
for the disc as above. 



compute 


8-7. Vibrations of a disc. 

Since we have been considering a disc, the natural transition is 
to the discussion of its normal modes of vibration. A disc differs from 
a membrane as a rod differs from a string; the flexural rigidity plays 
the predominant part. In polar co-ordinates r, <f> the vibrations of a 
disc are governed by the equation 

d*z Eh* w i d i an* 

H dfi ^ 3(1 - o*) [dr* + r dr + r 2 2 “ °* 

T^ie significance of the symbols is the same as before, with the addition 
of p which is the density. The analysis leading up to the equation can 
be found in Rayleigh, Theory of Sound , Vol. 1, § 218 and elsewhere. 
Ihe problem is m some respects the two-dimensional analogue of the 
transversely vibrating rod. The boundary conditions for a clamped 
edge are simple. Since there is neither displacement nor slope we 



In the case of a free edge, which implies the absence of shear and 
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bending, the conditions are more complicated. They are usually given 
as 




a idh i dz i dh i i - a a 2 (dz z \ 

dr \ dr 2 ^ rdr~^ r 2 d<f> 2 ] + r 2 d<j? W r) ~ ° 


We assume, as usual, that the normal modes of vibration can be ac¬ 
counted for by the adoption of a periodic time factor exp iwt. For the 
azimuthal angle (f> we assume the presence of a factor exp in<f> with a 
view to Fourier analysis. We accordingly make the substitution 

z = R exp?(to£ -}- n<f>). 


where R is a function of r alone. The equation can then be written 



1 d 


n 


2 \ 2 


+ R = k*R, 

r dr r 



3/>o> 2 (l — CT 2 ) 

M 2 


Alternatively 

{ e -2»(02 _ W 2)}2fl = lAR t 

which is clearly equivalent to the pair of equat ions 

{e-^O 2 - n 2 )}# = + k 2 R, 

and any solution of either of these is acceptable. Reverting to the 
independent variable r we have 



The value of R is thus given by 

R = AJ n (kr) + BY n (kr) + CI n (kr) + DKJkr). 

The fulfilment of the boundary conditions for any stipulated mode of 
support leads to a transcendental equation which ultimately determines 
oj and the frequency of the normal mode of vibration. With so many 
terms present, this may occasionally call for a certain amount of mani¬ 
pulative skill. As an illustration we consider a disc which is complete 
and free, meaning at rest and unsupported in a gravitationless field. 

The disc being complete, the displacement of the centre must be 
finite and the slope there zero, as already explained. This compels us 
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to discard the functions Y„ and K,„ and to work with the more limited 
solution 

R = AJ n (kr) + CI„(kr), 

equivalent to 

2 = (AJ n (kr) -f CI n (kr)} sin n<f> cos tut. 

The first of the two boundary conditions can be rewritten as 


(PR 
dr 2 


+ _ 5? _ - R — ( i-f? /r i- — n 2 tf) = 0. 

r <?r r 2 r 2 \ dr ) 


As far as J„(kr) is concerned, and in virtue of the equation which it 
satisfies, this provides us with 

-R7„(Ar) - {krJ.\kr) - «,({/]}. 


Similarly I„{kr) provides us with 


m n (kr) - {krl„\kr) - nH„(lr)}. 


These can be modified by using the recurrence relations 

3"</ n = M*?n 3cJn + l > %?n = n “f" «+!• 

We conclude that the first boundary condition at r = a is equivalent to 

_ 0 _ (n 2 ~n — m)J n (ka) + kaJ n+1 (ka) _ k 2 a 2 
A (n 2 — n + m)I n (ka) — kal n+1 {ka)' m ~ 1 — d 

The secondary boundary condition can be written 

* I™ -u \ dR _ w * p! « 2 (i - *) n 

drldr 2 rdr r 2 J r 2 \ dr 7 ~ 

The contribution from J n (kr) is 

_z. 2 ^n_ n 2 (l~a)ldJ n J n \ 

dr 7* 17/T 7J* 

By the use of the recurrence formula we have, after a little reduction, 

( n ( n ~m~ n 2 )J n (kr) + (m + n 2 )krJ n+1 (kr)}. 
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Similarly we can show that the contribution from I„(kr) is 

LZI? {„(» + m - n-)l n (kr) + (m - n*)lcrl n+x (kr)}. 

The elimination of the ratio — C/A then gives 

(,< 2 — M — w)./ w (c) + cJ n+1 (c) _ w(n 2 - n + »i)J„(c) - (n 2 + m)cJ n+ ,(c) 
(n 2 n + m)l n (c) — cl n+1 (c) n(n 2 - n - -f ( n 2 - ;>0c/ n+1 (c) ’ 

where the argument is c = ka. This determines k and thence cu and the 
normal modes of vibration for a given value of n. 


EXERCISES 


1. Verify that the last relation can be written 
£ , 2 » 2 (« ~ *> In J n 
C f n +l^n+1 


n* — n* + m 
2m 


{/— + /=-) +»<»- nf/ 5 - - /=-)• 


n+1 ° n+V w n+1 •'fi+l 

2. Deduce that when there are no nodal diameters, the relation becomes 





Establish this result independently. 

3. In connexion with the last equation, prove that the function xI Q (z)/I ,(x) 
has the value 2 when x is zero anl steadily increases with x. Hence prove that 
the equation can have no root smaller than the first zero of «/ 0 (x), but must have 
a root before the first zero of J x (x). With a = 0-3 solve by the aid of the extracted 


values 

x 

Jo 

Ji 

/o 

/. 

2-99 

— 0-25664 

034278 

4-8414 

3-9179 

300 

— 0-20005 

0-33906 

4-8808 

3-9534. 


Hence compute for a disc of radius 2 cm. and thickness I mm., density 7*6 gr. 
per c.e. and E = 2 10 9 gr. per cm. 2 . 

4. When there is a single nodal diameter, derive the relation from No. 1 
above and establish it independently. 

How much of the modified argument in No. 3 can now be applied as an aid 
to solution? 

5. If the disc be complete and clamped at the rim, prove that the discrimi¬ 
nating equation is 

■/„(*») **(*») ’ 

Give an alternative form with orders n and n + 1. 
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8-8. The non-uniform disc. 

Having examined the vibrations of a uniform disc it is natural to 
examine the problem when the hypothesis of uniform thickness is 
abandoned for a certain variability. The law of variation needs to be 
fairly simple, say according to some power of the radius, and if we 
put h = c(a/r)* we can vary the rate of change of thickness, as we move 
along the radius, within wide limits. If 5 is negative the disc has zero 
thickness at the centre; but of course the disc need not be complete. 
Similarly if s is positive the central thickness is infinite. The disc may 
then either be incomplete or be considered as mounted at the end of a 
spindle. 

The curious fact emerges that the solubility of the problem by 
Bessel functions depends on the number of nodal lines, and one par¬ 
ticular case comes out in terms of the elementary functions. The 
fundamental equation is given in the texts as 



which closely resembles the corresponding equation for a non-uniform 
rod. In the present case we have 



so that we can rewrite as 




+ 4. l'« 

W^rdr r*J r dV 

where 


m 2 R j 

r 2 j 


= k*Rr-, 


z = R sin m<f> cosa»L z« — 3/xo 8 (l a 2 ) 

EcW ' 


Further progress is greatly facilitated by the symbolical method, so 
that * 


c __ m 2) e -e(3, + 2){cp» __ 

whence 

e-2*<. + 2){ ( c y _ 3s _ 2)2 _ m 2j ( rp. __ m2)R = ^ 


° f a SOlution lies in expressing the operator as a square; but 
unless there is some re-arrangement of the factors we shall merely be 
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left with the case where s is zero, as the reader can easily verify. We 
accordingly re-arrange the factors in the order 

»> — m — 3s — 2, 0 — m, 5 -f- m — 3s — 2, & -f m. 


We can then write the equation as 

e -*(*+ 2 )(£ _ m _ 2 s)(0 — m + s + 2)c- fi( * f 2) ($ -f m—3s—2)($ -f m)R=k A R. 

The operator is now a square provided 2m = s -f- 2. We accordingly 
remove s and write 

{c -2m9 (0- — 5m -f 4)(5 -f m)} 2 R = k*R, 


which is equivalent to the two equations 

d 2 R 5 — 4m (JR bn(4 — 5m) 
1* + r dr + ( ? 


± J 



These are plainly comparable with our standards 3-6(2) and 7-3(1). 
The comparison gives 

1 — 2a = 5 — 4m, y = m, 

Py = k, a 2 — n 2 y 2 = m(4 — 5m), 


whence we derive 

a = 2m — 2 = s, 




It appears that this gives no solution when m is zero, corresponding to 
a thickness proportional to the square of the radius. The next per¬ 
missible value of m is unity and leads to 

m = 1 = n = y, a = 0 = s, ft = k. 

The disc is of uniform thickness and the solution is given in functions 
of the type J^kr). A certain interest attaches to the value m = 4. 
It leads to y = 4, a = G = s, p = k/4, n = 5/2, so that the functions 
are of order half an odd integer. 



APPLICATIONS TO HYDRODYNAiMICS AND ELASTICITY 151 


EXERCISES 

1. Prove tliat the functions are of integral order only when m is 1 or 2. Prove 
also that the functions are of order half an odd integer only when in is 4. 

2. By writing the factors in the order 

S 4- m — 3s — 2, S -f m, % — m — 3s — 2, S — w, 

prove that the operator can be made a square provided 2 -f 2m -f s = 0. This 
leads to 

(PR 4m -f- 5 dR (w(5m + 4) P 1 
dp r dr ~ I r 2 ^ ^ = 

Discuss the various types of function involved in the solution under different 
conditions, proving incidentally that m cannot be zero and that functions of 
order 3£ may appear; also that only two values of m can lead to functions of 
integral order. Note that the argument of the functions decreases with increas¬ 
ing r. 

3. Prove that the only other re-arrangement of the factors that will lead to 
a square operator is to put (** - m*) first. Show that this gives s the unique 
value — 1 and the thickness is proportional to the radial distance. Verify that 
the equations are 

(PR 3 dR 1 - m* IP 
dr* + r dr r* ± r = 

These are soluble for all values of m and lead to functions of the typo 

r'J im {2kr\). 


CHAPTER IX 


Bessel Coefficients. Integrals 
and Expansions 


91. Bessel coefficients. 

The functions J n {x) of integral order are sometimes called Bessel 
coefficients since they occur as coefficients in a certain expansion 
which we shall now consider. The left side of the identity 

exp {£x(t — r 1 )} = exp (£rf) exp 

can be expressed as a both ways infinite series of ascending and de¬ 
scending powers of t, a so-called Laurent series. We take the product 
of the two expansions 

exp(^) = 1 + (te) A + (te)2 l l -f. . . . + (|*)»£ + .... 

1! 2! n! 


= 1 - «*) £ + ~+... (- 5 *)" • • 

The independent term in this product is 


1 ( I *) 2 ■ ( I ^) 4 

( l !) 2 ~ r ( 2!) 2 



The coefficient of l n is 

(H n _ ih x ) 2 , (H 4 

n\ I 1 (n +1) 1.2(n -f l)(n + 2) 


= Jni?)- 


Similarly the coefficient of t~ n is 

(~^) n h (t x ) 2 ■ (l*) 4 _ 

n! I l(n-fl) 1.2(n+l)(n + 2) 


. J = (—)V n (x) = J_ n (x). 


We accordingly have 

exp{£x(< - <-i)} = 'Z°_J n Jn(x) 

= J Q (x) + - |->) + J 2 (x)(t 2 + «-*) + • • • • 

162 
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This last result is extremely fruitful and we proceed to make deduc¬ 
tions from it. Putting t = exp id we have 

exp ( ix sin 6 ) = J 0 (x) -f- 2 J 2 {x) cos 20 -f- 2 J 4 (x) cos 46 -f- . . . 

+ 2 i{J x (x) sin 9 -f J 3 (x) sin 30 + 

Separation of the real and imaginary parts in 

exp (ix sin 6 ) = cos (x sin 6 ) -f- i sin ( x sin 6 ) 


gives the two series 


( 1 ) cos (x sind) = J 0 (x) -f 2 J 2 (x) cos 2 9 + 2 J 4 (x) cos 40 -f . . . , 

(2) sin (x sin 0) = 2{«/ 1 (x) sin0 -f- J 3 (x) sin 30 -f- . • . }• 


Replacing 0 by its complement i -it — 0 we have correspondinglv 

(3) cos (x cos0) = J 0 (x) - 2 J 2 (x) cos20 -f- 2 J 4 (x) cos40 - . . . , 

(4) sin (a; cos 0) = 2{J 1 (a:) cos0 — J 3 (x) cos 30 


These series are usually associated with the name of Jacobi. Any one 
of the four may be regarded as a Fourier series. Two are in even cosines; 
the others are in odd sines and odd cosines respectively. The law for 
the formation of the coefficients then gives 


(5) 7rJ 2ri (z) ==^ cos(xsin0) cos2n0J0 

= cos (* cos0) cos2u0J0, 

(G) t tJ 2 „ +1 (x) sin (ar sin0) sin(2n -f- 1)9 d9 

= (— ) n f Q sin (a; cos 0) cos(2n -f- 1)0J0. 

In each case the limits can be taken as 0 to if we double the result, 
it 0 be replaced by its supplement tt - 0 in cos (x sin 0) cos h9 there is 
a change of sign if n is odd. Hence we have 


(7) 

similarly 


(0 150) 


f 0 cos (x sin B) cos nBdB = „,/„(*) („ even), 

= 0 (n odd), 

f Q sin (a: sin0) sin nddd = 0 („ even), 

= 7rJ n (x) (n odd). 
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It follows by addition that 

(8) cos {nd — x sin 0 )d0 = ttJ „(x) (n integral). 

These integrals are associated with the name of Bessel. The last of 
them is of historic interest as being the starting point of Bessel’s in¬ 
vestigations. It arose in connexion with an astronomical problem 
concerning what is known as the “ eccentric anomaly As an example 
of the application we have the following. 


9-2. Application to the cycloid. 

Problem 18 .—A circle of unit radius rolls on a straight line. The original 
point of contact describes a cycloid. Express this in a Fourier series. 

With the origin at the initial point of contact, the parametric co¬ 
ordinates of the tracing-point are 

x = 6 — sin#, y = 1 — cos 6 . 

The curve has a base of 2ir and is symmetrical about its middle ordinate; 
the series is therefore cosines only, so that 

y = °o + a i cos x -f- a 2 cos2x 


The limits for x are 0 to 27 t and the same holds for #. Hence 

a n = - f y cos nxdx = — f cos(nd — n sin0)(l — cos 6)-d9. 

7T Jo ' 77 Jo 

It remains to evaluate this integral, which we shall do by using the 
identity 

(1 — cos#) 2 = 2 — 2 cos 9 — sin 2 0. 

If we differentiate Bessel’s last integral with respect to x we have 


f sin (nd — xs'md) sinOdd = 7rJ n '(x). 
Ja 


A second differentiation gives 


f cos (nd — a; sin#) sin 2 6d9= — ttJ , 
Jo 


// 


(x). 
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Moreover, 

2 f cos (tid — x sin 6 ) cos Odd 
Jo 

=J cos{(n -f* 1)0 — x sin $}(19-\-J cos {(«— 1)0 — xsin0j dO 

= i(x) 

277-n . . . 

=-«/„(*)• 

The Fourier coefficients are therefore given by 

27m 


^77a„ = 2 tt«/ n (n)-«/„(«) + 77^ n ''()i) 


71 


or, 


= 2.7 „"(»). 

If we put x = n in Bessel's equation we have 

•'."<») + 1JM + (i - = o, 


whence 


a n = — - J n '(n). 

n 


The solution is completed by verifying that the mean ordinate a 0 is 
3/2. The required Fourier series is therefore 

y = f ~ 2{«7 l '(l) cosx + U 2 '(2) cos2x + £./ 3 '(3) cos3x + . . . }. 


EXERCISES 

1. Obtain the recurrence formulae by differentiating the relation 

oxp{lx(l - !->)} = 

with respect to x, or t, and comparing the coefficients of corresponding powers 

2. By modifying the variables in the above relation, prove that 

oxp{|x(l + <-*)} = 2* co < n / n (x). 

Replace 1 by r» and dcduco 7 n (x) = /_ n (x). 
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3. From the Jacobi series deduce 

1 = J 0 ( x ) + 2J 2 (x) + 2J i {x) + ... , 
cos x = J 0 (x) - 2 J 2 (x) + 2 J,(x) - . . . , 
sinz = 2{J l (x) - J 3 (x) + J,(x) -...}. 

Establish corresponding results for the hyperbolic functions cosh and sinh. Bv 
means of one, two differentiations, prove that 

ix = J x (x) + 3J 3 (x) + 5 J 6 (x) + ... , 
ix sin x= 2V 2 (*) - 4V 4 (x) + 6V s (x) - . . . , 

ix cos x = IV, (x) - 3V 3 (x) + 5*J s (x) - . 

4. From the Fourier series for the cycloid, by taking the highest and lowest 
points, prove that 

i = «/.'(!) + W*l 2) + i-/,'(3) + . . . , 

i = J t '( 1) - *7/(2) + J/ 3 '(3) - . . . . 

Hence 

i = JiV) + W(3) + + ... , 

i = WtW + i«V(4) + K.'(6) + . . . . 

5. A trochoid is described by a point whose distance from the centre is c, 
less than the radius which is unity. If the curve be expressed in a Fourier series, 
the lowest point being on the y axis, prove that the coefficient of cos;tx is 

2c 

a n = “ - J n \nc). 

The parametric co-ordinates are 

x = 0 — c sin 0, y = 1 — c cos 0. 

[The problem on the cycloid introduces functions of equal order and argument. 
Considerable interest attaches to the theory of such functions, especially when 
the order is large. The convergence of the four series mentioned in Question 4 

above is extraordinarily slow. Even when n is as high as 50, the value of -J n '(n) 
is about 0 0006.] n 

2 rl" 1 r* 

C. J 0 (x) — - cos (x sin 0)70 = - / exp (ix cos 0)70, 

Two Two 

Iq(x) = - f* cosh (x sin 0)rf0 = i f exp(xcos0)rf0. 

TWO Two 

7. Prove that the absolute value of J n (x) is less than unity. 

8. Establish the relation 

exp (x cos 0) = I Q (x ) -f- 2/ 1 (x)cos0 -f 2/ 2 (x)cos20 + . • • • 

9. If the anode voltage in a rectifying valve is Ecoscot and the current may be 
taken as / = A exp (IE cos cot), prove that the mean value of the current is AI 0 (bE) 
and that the root mean square as A{I 0 (2bE)}K 
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9-3. The Poisson integral. 

The function J n (x) can be defined by an integral, closely allied to 
the previous forms, which sometimes bears the name of Poisson and is 
otherwise referred to as Bessel’s second integral. It is established as 
follows. The general, or (r -f- l)th, term in the expansion of J„(x) is 

(-m*) n+2r r = 0 ! o 

r(n + r+l)r(r + l)’ ’ ’ ’ ' * * * 


We propose now to multiply the expansion of cos (x cos 6) by sin 2n 0 and 
integrate termwise, n not necessarily being an integer. As groundwork 
for this we have from the Beta functions 1-6(2) 


2 ['"cost'd sin Zn ddd = r( ^+i )r(n + * ) . 
J o r(» + r+l) 


Multiply top and bottom by T(r -J- 1) and use the duplication formula 
1 - 8 ( 1 ) 

T(r + £)I> + 1)2 2 ' = Vn r(2r + 1). 

We then have, if r is an integer, 

fcos 2 '0 sin ^QdQ = 2 V^ r (« + ^)r(2r + 1) 

Jo 2 2 T(w-}-r-f- l)P(r-f 1) 


The (r -j- l)th term in the 


expansion of cos (z cos 0) is 


(~y^ r 

P(2r + 1) 


cos 2r 0, 



Multiply this by sin 2 "# and integrate; wo get 

^ fcos 2 '# sin *"ddd = ± 

P(2r -f l)Jo I> + r + l)r(r +• 1) 

It will be observed that this is the corresponding term of J n (x) mul 
tiplied by 


r ( n + $)V 7T (&)~ n - 

We conclude that 

(1) J n (x) = f cos (x cos 6) sin 2 "#d#. 

This is the formula referred to. The integrand is unchanged if ( be 
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replaced by its supplement tt — 9. and we conclude that the integral 
can be doubled over the half-range, so that 

JJx) oc 2 f* cos (x cos 9) sin 2 n ddO. 

Jo 


By using the complement \tt — 9 we have alternatively 

o 

so that 

n 


JJx) oc 2 f cos (x sin 9) cos 2 n 9d0, 
Jo 


( 2 ) 


j ( x ) = x ^ 

T(n -f- h)\ZrrJo 


f cos (x sin 9) cos 2 n 0d0. 
Jo 


9-4. Application to the circle. 

Students of Fourier analysis, as performed by the integral calculus, 
are usually confined to diagrams consisting of straight lines, with 
occasional excursions into exponentials and parabolas. The circle 
never appears, and here is the explanation. 


Problem 19 .—Express the upper half of the circle x 2 + y z = tt 2 in a 
Fourier series. 

Taking the radius as tt is a mere matter of convenience; by a 
change of scale it can later be made anything we like. The mean 
ordinate is \-rr and since the curve is symmetrical about the y axis, 
the series is cosines only. We accordingly have 

y = £ 7 r 2 + a t cosx + a 2 cos2x + . . . , 

2 r 

where a n = - / y cos nxdx. 


Turning to polar co-ordinates 

x = tt cos d, y = 7T sin 9, 

we have 

9 r 0 

a n = - / (tt sin 9) cos (rnr cos 9)(— tt sin 9d9) 

TT J \' 

\ 


= 2t t f cos (mr cos 9) sin 2 9d9. 
Jo 
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Note the change of limits when 9 replaces r. The integral is equivalent 
to 



nnw* 

(bln) 


Jl(nn) = 


- J l(nn). 
n 


The required series is therefore 

M_. Ji(' 27T ) 


J i(3w) 


77* . (</1 ('ll I , viiufii . «||un/ « , 

= — 4- 7T | ——-—- cos x H--—- cos 2j -f- — cos3r -f 

4 11 



9-5. Modification of the integral. 

A useful modification of the formulae can be made by first trans¬ 
posing the factor x n . Wc then differentiate with respect to x, taking 
advantage of the formula 

'(L = — x ~ nJ »+A x )- 


Afterwards we replace the factor x n . This gives 

f sin (x cos 9) sin 2n 0 cos 6d9 

v - Jo 

_ (£*)" 


(1) - r /„ ( 1 T, , 

A (« + bV 77 


( 2 ) 


r(n + hW* 


f sin (x sin 9) co s 2n 9 sin 9d9. 

Jo 


9-6. The Lissajous figure of eight. 

As an illustration of their use, consider the following variant of a 
familiar problem. 

Problem 20. —A moving lighl-spol simultaneously executes the two per¬ 
pendicular simple harmonic motions 

x = n cos cot, y = c sin2coL 

It is required to find the Fourier series for its path in the time in¬ 
terval 0 to n/co. 

It is well known that the point repeatedly describes a figure of 
eight, one of the so-called Lissajous figures (fig. 9). The time interval 
tt/co is a vertical period, or half a horizontal period. Initially t = 0 = y, 
7r. Finally t = njco, y = 0, x = — n. At half-time we have cot = 
fom = ° = V% The CUrVe is skew and the requisite series is of the 

y = 6 1 sinx-l-6 2 sin2a;+ .... 
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The coefficients are given by 

2 r n 

b n = / y sXnnxdx. 

TT J 0 

Put 

wt = 0, x = 7T cos 6 , y = c sin 20, dx =—wsin Odd. 

The limits for 6 corresponding to x — tt, 0 are 6 = 0, respectively. 
Hence 

2 r° 

b„ = — / (c sin26) sin (mr cos 9)(—7 t sin Odd) 

7 T J i" 

= 2 cf sin(M7r cos0) sin 2 0cos0rf0. 

Jo 

This is patently the upper of the two forms 9-5(1). Comparison then 
gives 

*•-* W- )E K=)-?« n - ) - 

The required series is accordingly 

— = J o(tt) sinx + fa?aft 1 *) sin2x -f- \J 2 {^it) sin3x -f- . . . . 

2c 

The horizontal amplitude is changed from 7r to a if x be replaced by 
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9 7. Transformation and an application to the circle. 

Bessel’s second integral admits of various transformations. They 
mostly amount to substituting, in the first form t for cos0, in the 
second form t for sin 0. In either case we get 

(1) J„{x) = ( j X) / (1 — cos xtdt. 

Here again the integral can be doubled over the half-range 0 to 1. 
Moreover, if we follow the same technique of first transposing the 
factor x n and then differentiating with respect to x, we derive 

(2) ./„ +1 (s) = ( j z >" / «1 - (=)"-» sin xtdt. 

1 (n 2 ) v 77 

As an illustration of the application we offer the following. 


Problem 21 .—It is required to find the Fourier series for the upper 
half of the circle defined in polar co-ordinates by the equation r = 
2 n sin 9. 


The diameter, as a matter of convenience, is 27 t. The origin is on 
the circumference. The Cartesian co-ordinates arc 

x = 2n sin 2 9, y = 2-rr sin 9 cos 9. 

The curve is symmetrical about its middle ordinate and the required 
series hence has the form 

V = + a, cos* -f a 2 cos2x -f . . . . 

The coefficients are given by the rule 

2 

7T J 0 


<*n — — f y cosnxdx 

'TfJ 0 


2 r n l A 

= ~J 0 (*■ si *20) cos (2/i7r sin a 0)(27T sin 2000). 

Further progress is facilitated by writing 

cos (2nrr sin 2 9) = cos(«tt — rnr cos 2 9) 

. = (—) n cos(rj7r cos 20). 

Ihe relation becomes 

r*t 4 

a " = *”J 0 ^" cos ( n7T c °s20) sin 2 20 00. 
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We now make the substitution 

cos 26 = t, sin20 = (1 — < 2 )*, —2 sin2 Odd = dt, 

whence 

a n = —2 ir(—) n J (1 — t 2 )* cosmrtdt. 

Comparison then gives 

(—J-HUh/w 7 , x (-)" 7 / X 
a n = - - ' ~ - ItJAMt) = - 7rJ 1 (M7 t). 

\lX7T n 

This can be checked against the previous circle-problem by changing 
the origin. 


EXERCISES 


1. From the former of the two circle-problems establish the following iden¬ 
tities: 

J = J,(r) + J./,( 3rr) + + • • • • 

= J 1(77) — \J i( 2 t 7) -f- \J i( 377 ) • • • • 

1 — |77 = J i(2t7) + iJ i(477) + $«/ 1(677) + • • • • 

2. An elliptic arch has a span of 40 ft. and rise 18 ft. Taking centre line and 
ground level as axes, prove that the Fourier series, of period 40 ft. for its outline, 
is 

y TZ 7ZX 2 tzx 

— = - + ^i(~) cos 2Q + W i(2tt) cos -go" + • • • • 


3. From the Lissajous figure of eight, deduce the identity 

~ = •/gC 71 ) — J«/ 2 (37 t) + £</ 2 (57t) • • • • 

4. If the position of a moving point is defined at time l by the relations 

x = a cos at, y = c sin 3at, 

prove that its path in the time interval 0 < at < can be expressed as 


- ^ = J 3 ( 77 ) COS^+ J Jj(277) COS + JJ S (3t7) COS 


TZX 

a 


2 tjc 


3 TZX 


+ • 


6. From the last result above, deduce the identities: 

J = J 3 (tt) + §«/ 3 (3-) + J«/ 3 ( 5 ~) + • • • • 
i = J 3 (277) + 3 ( 677 ) + iJ,(1077) 

i = J 3 (4t7) + 1-7,(877) + J^,(1277) + . . . . 
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6. Multiply the identity cos 20 = 1 — 2 sin 2 0 by cos (x cos 0) and integrate from 

2 

0 to thus deducing J 0 + J»= - J x . Treat similarly cos 40 = 1 — 8 sin 2 0 
4- 8sin 4 0. Hence deduce generally that J 2n is expressible in terms of J 0f J l% . . . 
^ n* 

7. Establish the relations 

I Ax) = (- — -) f* cosh {x sin 9)cos 2n 0(/0 

n \T(n+1 )VtzIJ 0 

p\n 

= { )J cosh (x cos 0) sin 5n 0c/0 
= { } \j\ 1 — f 2 )" - * coshterff. 

9-8. Application of Lommel integrals. 

We now transfer our attention to instances where Bessel functions 
figure in the integrand. The cases of Fresnel integrals and Lommel 
integrals were both mentioned earlier in the book; but as no applica¬ 
tion of the latter has yet been studied, it is desirable to fill the lacuna. 
This is best done by an example in heat conduction. 

If <f> be the temperature at any point, we have from 8-4(1) 

1 9 J, \ i d~<f> | 1 3 2 <£ _ ps d<f> 
rdr\ dr I 3z 2 + r 2 d0 2 k dt 

If we can regard <f> as independent of z and 6 , the conduction takes 
place radially and we have 

3 2 <f> , l d<f> _ ps d<f> 
dr 2 r dr k dt 

Assuming for a falling temperature that <f> has the form R exp (—mt), 
where R is a function of r alone, we have solutions of the type 

<f> = {AJ 0 (kr) + By o (tr)}e _m ', Jfc* = £!!?. 

K 

Problem 22. A solid bar of radius c radiates into a medium at zero 
temperature. Assuming that the temperature in the bar is indepen¬ 
dent of the axial distance and of the azimuthal angle , it is required to 
examine the variation in temperature distribution xoith time. 

Discarding the Y function because the bar is solid, we adopt solu¬ 
tions of the type 

<f> = AJ Q (1cr)e~ mt . 
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For the emissivity we have the surface condition 


whence 

or, 




hJ Q (kc) + KkJ 0 '(kc) = 0, 


■vm += °» 



The roots of this equation fix the permissible values of k, and thence 
m. If we build up a series solution 

<f> = 'LA v J 0 {k l r) exp (—m x t), 

the presumption is that initially the temperature distribution is given 
by 

(2) <f> 0 =f(r) = 'ZA 1 J 0 (k 1 r). 

The problem poses two questions at this stage. Firstly, whether an 
arbitrary function f(r) can necessarily be expressed in a series of func¬ 
tions of the type J 0 (k $ r) where the k, are assigned. We assume that it 
can, leaving the pure mathematicians to argue about the conditions of 
its validity. Secondly, how the coefficients A t are to be determined; 
here we assume the legitimacy of termwise integration and we utilize 
Lommel’s integrals. Consider the relation 4-6(1) 

(a 2 — pr)j^ rJ 0 (ar)J 0 (pr)dr = c{pj 0 (ac)J 0 '(Pc) — aJ 0 '(ac)«7 0 (/fc)}. 

The right side is zero if a, are two values of k determined from (1). 
It follows that we have/(r) expressed in a series of orthogonal functions, 
and if we multiply (2) by rJ Q [k u r) and integrate, we have 

AS rJ 0 {k u r)J 0 (k c r)dr = 0, u =1= v. 

Jo 

Turning to the companion relation 4*6(3) 

I" rJ 0 2 (kr)dr = hc?{J 0 2 (kc) + «/ 0 2 (^ c )} 

Jq 

= Jc^l + pWitc), 
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we conclude that 

f Q rf( r )J 0 (k s r)dr = A S J^ rJ 2 (l; s r)dr 

= U s c*(l -f- p 2 )J 0 2 (k x c). 

The method of determining the coefficients is analogous to that em¬ 
ployed in Fourier analysis. Once the equation for k, has been solved, 
A, may be considered as known. Whether one could perform the in¬ 
tegration on the left in any particular case is another matter. Further 
illustrations will be found in H. S. Carslaw, The Conduction of Heat 
Chap. VII. 


9-9. Lipschitz’s integral. 

The number of known integrable expressions involving Bessel 
functions in the integrand is already somewhat extensive and there 
are accretions almost every year. Most of them require complex 
analysis for their establishment; but there are two exceptions which 
we will now give. The first is due to Lipschitz. As groundwork we 
require two elementary integrals. 

(i) Treating fe~ ax cosbxdx symbolically, we have 

jr e-°* cos bx = e~ ax _J— cos bx 
u D — a 


>—ox 


►-ax 


a 2 _j_ &2 & + a) cos&z — ^ (b sin 6* — a cos6x). 

It follows that for the infinite integral we have 


(ii) If we write 


[ e ax cosbxdx — —__. 

jQ a 2 + 6 2 


Ja 2 - j- 6 2 cos 2 * I ' 


<f(tan<) 


(a 2 -f- 6 2 ) -f- a 2 tan 2 * 

1 .( a tan* 


a(a 2 + 6 2 )* 


we have the definite integral 


tan~ J / —fL^ an ^ ) 

l(a 2 + 6 2 )‘j’ 


Jo d 


dt 


7r 


+ 6 a cos 2 < 2a(a 2 + 6 2 )*’ 
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Taking Bessel’s first integral, coupled with a change in the order of 
integration, we have 


f e~ az J 0 (bx)dx = [ e~ 
Jo J o 


ax 


2 f iv 

dx - I cos(6x cos t)dt 

7T J 0 


2 rl* r™ 

= -/ dt e~ ax cos(bx cost)dx 

tt j o J o 


_ 2(1 /**’ 
~~ TT^o a 


dt 


2 + b 2 cos 2 1 


1 

(a 2 + 6 2 )*’ 


EXERCISES 


1. Prove 


r° i r® 

J^ j 0 (bx)dx= J^ J 0 (z) dx = 1 


2. Multiply the expansion of J 0 (bx) by e ^ and integrate termwise, using the 
relation 1-8 Ex. 1, 


j dx = a“ n r(n). 

Jo 


Hence establish Lipschitz's integral. 

3. Differentiate the integral with respect to b and deduce 

T® 6 

J 0 = (a , + pf*- 

Transpose the factor b and differentiate again, utilizing the relation 

{tr*J n m = -t~ nJ n + 1 ( 0 - 


Establish that for n integral 


r~ , . ( 26 ) n r(n + 4) 

J 0 e-“J n (6x)x"dx = VTr(fll + fcl) „ + r 




2a(2b)*r(n + 3/2) 


J. e-‘J n (br)^'dx= VA „, + j.p+w 


Let a tend to zero. Is this legitimate? 
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9-10. Sonine’s integral. 

Sonine’s first finite integral is readily established by the method 
previously employed for Bessel’s first integral. We propose to multiply 
the expansion of J n (x sin 6) by (sin0) ,,+1 (cos0) 2m+1 and integrate term- 
wise. The (r -f- l)th term in the expansion is 


(-) r (i*)" +a ' 


r(» + r+l)r(r + l) 


(sin 0) n+2r , 


so that the essential part of the integral is a Beta function and expres¬ 
sible in Gamma functions. We have 


2 f ‘"(sin0)2«+2r+i (cos0) 2m+I (16 = P(n + r + + *) 

r(#+r + m + 2) 


The reason for taking (n -f- 1) as the index of sin 6 now appears. It pro¬ 
duces the factor T(n -j- r + 1) which cancels with the corresponding 
factor in the denominator. The net result of the integration is 

H r (i*) n+ 2 T(m + i) 

F(r + l)F(n + r + m + 2)‘ 

Inspection shows that if this were multiplied by 

'jHl+l 

2"T(m + 1) 


we should have the (r + l)th term of J n+m+l (x). We conclude that 


(1) *^n+m+l(^) — 


x m + 1 ri* 

2^ f(m 4- l) Jo Jn{x Sin 0)(sin 0)" +1 ( cos 0) tm * l d0. 


EXERCISES 

- 1 ' P !®T° that the validit y of the result is not confined to integral values of m 
and n. Why is it necessary to stipulate that m, n > -1? Try the effect of put- 
tmg n -4 and replacing m by m — J. 

2. Transform the integral by the substitution sin 0 = t and deduce 

” + ” +,(t) = rr (w + 1) 1 J "W 1 - 


f {xfydl = ar l J x (x). 
o 


3. Deduce 
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Hence by differentiation prove that 


f J 1 (xt)( 2 dt= ar 1 J 2 (x). 
Jo 


Generalize by showing that Sonine's integral satisfies the relation 

JZ = -*-V n+1 (z). 


4. Evaluate the integral 


f J 0 (kx)(c 2 — ar) n zdx. 
Jo 


6. Show that the equation 


r l _cosxt = 

Jo l 1 - W 


is satisfied when x is any zero of J 0 . 

911. Weber’s discontinuous integrals. 

If we give an imaginary value to a in Lipschitz’s integral we have 

f e~ ia ?J 0 (bx)dx = ( b 2 — a 2 ) - *. 

•'O 

The left side is complex; but the right side is either quite real or pure 
imaginary, according as b is greater or less than a. Equating the real 
and imaginary parts we have 

/ J 0 {bx) cos axdx = 0, a > b 
Jo 

= (6 2 — a 2 ) -1 , 6 > a. 

f J 0 (bx) sin axdx = (a 2 — 6 2 ) - *, a > b 
Jo 

= 0, b > a. 

This is known as Weber’s discontinuous integral; it is a fair sample of 
a whole hierarchy of discontinuous integrals. 

Integrating the first form with respect to a from a = p to a = q, 
we have 

fj 0 {bx) ^ ~ dx = 0 , (p>b, q > b). 

X 


fubx) S l n - dx(=I, say) 
•'O x 


Hence the integral 
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is independent of a, provided a > b. It remains to see whether / is a 
function of b. Substitute ax = t; then 



This shows that I is a function of c only; but seeing that it is not a 
function of a, it is not a function of b either. We determine its value 
by letting b tend to zero, with J 0 {0) = 1. We then have 

7 = 

Jo t 

Incidentally there are at least fourteen independent, proofs of this last 
relation. We conclude that (if a> b) 

Cubx) dx = Itt. 

•'O X 

Lipschitz’s integral has found applications in the theory of poten¬ 
tial. Sonine’s integral and its near relations appear in some of Ray¬ 
leigh’s work, whilst Weber’s integral was used in an investigation of 
the electrostatic potential of a charged disc. In the main, generaliza¬ 
tions lead to the hypergeometric function, and possibly associated 
Legendre functions. 


9-12. Fourier-Bessel series. 

The type of series which we assumed in 9-8(2) as the expansion of 
an arbitrary function is known as a Fourier-Bessel expansion; it is 
a great stand-by in applied problems where a series expansion is re¬ 
quired. All the functions are of the same order, not necessarily zero 
as in our case, and the arguments are determined from some equation 
of condition that permits the use of Lommel integrals. We take in 
general 

f(x) = AyJ^x) -f AvJ n (k^x) -f . . . . 

Multiply each side by xJ n (k^c) and integrate from 0 to c. We have 
4-6(1) 

(* 2 ~ ^)f‘xJ n (ax)J n (Px)dx = c{pJ n (ac)J n '(pc) - aJ n '(ac)J„(^)}. 

(O 160 ) 


12 
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The right side is zero if a, ft are values of k determined from 

(i) J v (kc) = 0, 

(ii) J n \kc) = 0, 

(iii) PU n '(kc) + QJ n (kc) = 0, 


where P, Q are constants. The first of these was used by Fourier when 
n is zero. The extension to other orders was made by Lommel; and the 
third condition, which we employed in problem 9-8(1), is usually 
associated with the name of Dini. 

The coefficients are then determined from the second of Lommel’s 
integrals 

f o xJ„\ax)ix = |r={(l - + J„' 2 (ac)}, 

which takes different forms according to the type of boundary condi¬ 
tion employed. For example, if the first of the above conditions holds 
we have 

f xf(x)J n (k^t)dx = \A s c-J n '-(kfi). 
j 0 

It has to be admitted that relatively few functions give compactly 
integrable forms. An obvious exception comes from 

r t^+1 

J x n+1 J n (kx)dx = — J n+1 (kx). 


This shows that if f(x) is x n and we write 

x n = A v J n (k x x) + AJJk^c) + . . . 

we have 

^ JWM = J n (kc) = 0. 

k. 

In virtue of the recurrence relation 

kxJ n '{kx) = nJ „(kx) — kxJ n+1 (kx), 

we have 

A 2 c n “ l 

' “ k. J n+1 (M’ 

and the series for x n is thus determined over a range 0 to c. 
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913. The polynomial. 

It is sometimes an advantage to have all arguments the same and 
to allow the order of the function to vary. We accordingly consider the 
possibility of an expansion of the type 

A X J x (x) + AoJ 2 ( x ) + A#J z (x) + • • • • 

In this connexion wc can utilize the Jacobi series to give an expansion 
for x" and hence for a polynomial. For brevity we put 2 cos# = c. 
The trigonometry books then state that 

2 cos 20 = c 2 — 2, 2 cos 30 = c 3 — 3c, 

and in general 

2 cos n$ = c--+ n{n ~ 3 > c—< — «(«--1)(«-5) 

2! 3! ^ t ... . 


We can write the Jacobi series as 9-1(1) 
cos Axe = J 0 (x) - «/ 2 (x)(c 2 — 2) -f- J 4 (x)(c 4 - 4c 2 + 2) — . . . . 
sin Axe = Jj(x)c - J 3 (x)(c 3 - 3c) + J & (x)(<* - 5c 3 + 5c) - . . . . 

If in both cases we expand the left side in powers of xc, we can compare 
the coefficients of corresponding powers of c. In this way we derive, 
for the first three powers. 


1 — Jo( x ) 2J r 2 ( a ') -f- 2./ 4 (x) , 

i x = f ^i( x ) 4" 3</ 3 (x) -j- 5<7 6 (x) 4* . . . , 

£x 2 = J a (x) + 4 J 4 (x) + 9J 0 (x) + . . . . 


It can be left as an exercise for the reader to show that in general 


(Ax) n = 2* ( n + 2r)(n + r — 1)1 . 
0 r! 


n+2 


r(*)- 


9*14. Schlomilch expansion. 

•, * of almost any vibration problem shows the occasional 
necessity for a series of uniform order, the argument being proportional 

form 6 mnk ° f the tGrm ‘ ThiS meanS that WG SGck an ox P ansion of the 

f( X ) = Oq + Oi</q(x) -f- OjJ 0 (2x) -f 


• • • 
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If this were possible, we should have at the origin 

/(0) — «o = «n- 

We should also have by differentiation 

— f '( x ) — a \Ji( x ) + 2 o 2 J i(2x) + .Sffy/j(3x) -{-••• . 

Presuming that we require this to be valid over the range 0 < x < c, 
we substitute x = c sin#; the limits for 6 are then 0 to r. We now 
propose to integrate termwise the expansion 

J x (nc sin 6) = (\nc) sin 6 — (£nc) 3 & -f- ( Inc) 5 s ’ n ^ 

The general, or rth term is, apart from sign, 

_! (sin^-i 


1 . 2 . 2 . 3 


• • • • 


(£nc) : 


r(r)r(r+ 1) 


As part of the integration we have 1-6(2) 

*” Y[r)y/ir 


2 f Q sin »-*0dd = B(r, $) = 


Hr + i) 


The integral of our general term is thus 


r (rW-TT 


r(r)r(r + 1) r(r + i) 

The duplication formula gives 1-8(1) 

2 2 T(r + i)I> + 1) = V^(2r + 1), 

The integral of our general term becomes 

(nc) 2r-1 


Hence 


(2r)! 


f *' r , • m nc (nc ) 3 , (Tic ) 5 

Jo J ' (nc sn ' e)de = 2! 4T + _ 6! 


1 — cos nc 


nc 
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Reverting to the series for f'(r) we have 

I — cosc 


—J Q f'( c si 11 6)d0 = a Y 


+ a 2 


1 — cos2c 


which is better written 

F(c) = {a 0 — /(0)} + a x cosc -f a. 2 cos 2c -f- . . . , 

where 

F(c) = c f f'(c s'mO)dd. 

Jo 

Treating this as a Fourier series, we derive the coefficients as 
°o —/(°) = - f F(c)dc, a n = -f F(c) cos ncdc. 

7TJQ 77 'U 


The series for/(x) may now' be considered as known, at least formally. 
It is evident that in practice the difficulties will arise from the inte¬ 
grations connected with F(c). As an instance of a soluble case we may 
consider the parabolic form f(x) = -This gives f\x) = x, 

F(c) = c I c sin Odd = c 2 , 

Jo 

a o = - r = —, a n = - f cr cos ncdc = (—-l) n —, 
rrJ 0 3 7T JQ ' ' n 3* 

valid in the range 


0 < X < 77. 


CHAPTER X 


Allied Functions 


101. Functions of the third kind. 

It was suggested by Nielsen that in honour of Hankel the symbol 
II should be used to denote the function J„ ± iY n . We accordingly 
have the somewhat cumbrous notation 

(1) Hf'\x) = J n (x) + iY n (x), H n ( '-\x) = J n (x) - iY n (x). 

These being linear combinations of J„ and Y n are necessarily solutions 
of Bessel’s equations of order n. They accordingly satisfy exactly the 
same recurrence formula; as J n and Y n . Moreover, since Y n was de¬ 
fined by means of J n , it follows that any one of «/„, Y n , II „ must be 
expressible as a linear combination of the other two. There is the 
further relation connecting any two solutions of the same equation. 

The // functions find their chief application in the theory of the 
subject. The assessment of their importance in practice is probably 
largely subjective. The standard work by Gray, Mathews and Mac- 
Robert does not even mention them; but they occasionally appear 
in Rayleigh’s work, and they are fairly freely used in Carslaw, Mathe¬ 
matical Theory of the Conduction of Heat. They are readily calculable 
from tables of J n and Y n . and tables of orders 0,1, i are given in Watson, 
Theory of Bessel Functions. 


10-2. Kelvin’s functions. 

There are a round score or so of functions closely allied to Bessel 
functions. For the most part their interest is mathematical; but 
outstanding exceptions are the ber and bei functions and their con¬ 
geners. They play an important part in alternating current theory; 
they have also been used in the two-dimensional motion of a viscous 
liquid. 

The equation for I 0 (l) is 


dfy 

dt- 



ni 
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Put t = x(t)* and the equation becomes 



with the solutions I Q {xy/i) and K 0 (x\/i). The her and bei functions arc 
defined as follows. Since 


/ 0 (0 = i + (40* + 



> 


we have real and imaginary parts in 


I 0 (x\/i) = 



(lx)* (|x)« 

(2!)2 ' r (4!) 2 




+ <'{<**)*- 




= berx -f- i bcix. 



The names are convenient abbreviations for Bessel-real and Bessel- 
imaginary. Both berx and bei x are purely real for real x, and it is 
easily shown that both series are absolutely convergent for all values 
of x. Among the more obvious properties we have 

(3) berO = l, bei 0 = 0, 

r x r x 

(4) / x berxrfx = x bei'x, / x beixr/x = — x ber'x. 

•'o J 0 9 

the latter pair being established from the series. The graphs of both 

functions oscillate (fig. 10). The functions ker and kei arc similarly 
defined from 


( 5 ) K Q (xy/i) = kerx -f i keix, 

and there is yet another pair, her and hei, defined from Hankels 

unction. It can be left to the reader to prove that the four functions 

her, bei, ker and kei are the solutions of the fourth order differential 
equation 

{**(■* ~ S) 2 -f- c«>}y = 0. 

The application of the method of Frobenius then shows that the roots 
ot the indicial equation are r = 0, 0, 2, 2. It follows that there is a 
solution beginning with r°; this is berx. The second solution for the 
ero index certainly contains a logarithmic term; this is kerx. Simi- 
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larly the second solution for the index 2 certainly contains a logarithmic 
term; this is kei x. Finally, in accordance with theory, the first solution 
corresponding to the index 2 may or may not contain a logarithmic 
term. It happens that it does not, and the solution is bei x. All the 
functions have been generalized for order and argument. Those which 
have been mentioned are of zero order, in which case it is customary 
to omit the suffix denoting this. 



10-3. Current distribution in a conductor. 

Probably the simplest application is to the current distribution 
over the cross-section of a conductor. The argument is based on the 
two circuital theorems of Ampere and Faraday. Let H be the mag¬ 
netic intensity at distance r from the axis of a conductor of radius c. 
The work done in carrying unit pole round the circle of radius r is 
JV = 27 ttH. For a circle of radius r + dr we have W + dW. The 
difference dW tallies with 477 times the current through the annulus. 
Hence 

— (27 rrH)dr = 47r(27rr dr) a, 
dr 


where a is the variable current density. Thus 



I - (rH) = 47ra and 
r dr 
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Taking a longitudinal section of unit length, the annulus provides two 
rectangles measuring 1 by dr. The magnetic induction through one of 
these is pHdr. Taking the upper rectangle only, the potential drop 
along the lower edge is E = pa, where p is the resistivity. For the 
upper edge we have similarly E + dE, and for the line integral round 
the rectangle we have 



where p is the permeability. The elimination of H gives 

13 1 3c | 4 rrp da 

r dr \ dr) p dt 

Presuming that the alternating current density is given by 

da 


a = R exp i(col -f- <f>), — = itua, 

dt 


we have 


or, 


1 dldR 


r dr 1 dr 


- n = o, 

P 


%% + - ~ - ik-R = 0, & = 

dr 1 r dr n 


The solution is evidently 

R = C{ber(Ar) -f- i bei(X r)}, 

which remains finite on the axis, C being some arbitrary constant. 

Consideration of the phase of a shows that the current is not neces¬ 
sarily everywhere in the same direction, since the functions ber and 
bei may be of the same or opposite signs, according to the value of 
kr\ in fact, if w and accordingly k is large, the current density vector 
may revolve more than once. 

If a 0 denote the current density at the surface, we have 

a c oc C{ber (£e) + i bei (A-c)}, 

the time periodicity factor being understood. The total current through 
a circle of radius r is 

J 0 27r °rdr = ~~f'{kr ber (At) d(kr) + ikr bei (At) d(kr)} 


2ttC 

A* 2 


[Av{bei' (Icr) — i ber' (At)}]. 
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If we denote this by P, and let Q denote the corresponding quantity 
for the whole cross-section, we have 

P | 2 = /A 2 bei' 2 (At) + ber' 2 (A-r) 

Q | \c) bei' 2 (Ac) -j- ber' 2 (Ac)’ 

If < 7 denote the mean current density we have 

a m = — = ^ {bei'(Ac) — i ber'(Ac)}. 

TTC 1 kc 


It follows that the constant C may be considered known if either a c 
or a m is known. It is customary to write 

OcP = Q(R + *ojL), 


and call R the equivalent or effective resistance, L the internal self- 
inductance, ioL the effective reactance. Whatever they may be called, 
we have 


R -f- iwL = 


ikp ber(Ac) -j- i bei (kc) 
27rc ber' (kc) -f- i bei' (kc) 


The importance of the absolute value of the ratio P/Q lies in its 
ability to account for the “ skin effect ”, whereby at high frequencies 
the current tends to desert the centre of the conductor. As no adequate 
explanation can be given without the use of asymptotic expansions we 
postpone the matter. The various types of function have been fre¬ 
quently calculated and they are usually published in the British 
Association Reports. 


10-4. Asymptotic expansions. 

As these have several times been mentioned in the course of the 
book it is as well that a word should be said on the matter. It is usually 
important to know how a function behaves for large values of the 
argument, especially if tabulation is in view. Tabulation is a tedious 
business anyway, and a formula that will give an approximate result 
is invaluable, especially if at the same time it gives an indication of 
the magnitude of the error. The series expansion of a function may be 
perfectly valid for all values of the argument, and yet useless for 
purposes of computation on account of its slow convergence. 

The obvious suggestion for a way out of the difficulty is to expand 
the function in powers of x~ x , so that as x gets large the terms get 
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small. The answer is that such an expansion is not necessarily valid. 
A function that is defined by a differential equation is very much tied 
in its permissible ways of expression. The reader is probably already 
familiar with asymptotic formulae; thus (1 -j- n -1 )" is asymptotic 
to e. An asymptotic series is something different. In spite of the fact 
that a valid infinite series in descending powers of the variable possibly 
does not exist for some given function, it may still be possible by some 
means to obtain such an expansion formally. The “ means ” is almost 
invariably a contour integral using the complex variable z = x -j- iy, 
and the resulting expansion is divergent. All the same we can always 
say that the function equals the first n terms plus a remainder, the 
remainder being defined at worst as the difference between the function 
and the sum of the first n terms. Everything depends on this remainder, 
whose value is a function of n and the variable. If it can be shown 
that, for a fixed n, R tends to zero as the variable increases, then the 
first n terms constitute the asymptotic expansion of the function, and 
R gives the magnitude of the error. 

An outstandingly good example of an asymptotic series is provided 
by the error function. An integration by parts gives 

= i fz-'e-'dlj?) = - II -ftl J X . 


Making this definite we have 



and more generally 




2x 2r+1 


(2r + 



A repeated application of the last formula gives us 

/V*v* = ^(l_± + + \ 

2x \ 2X 2 (2j 2)2 (2x*j 8 ' ' * / * 

Inspection shows that for any value of x that may be considered large, 
say 10, the magnitude of the terms decreases with unusual rapidity, 
rurther inspection shows that, no matter how large x may be, the 
numerators will ultimately overtake the denominators and the terms 
will then be on their way to increasing beyond all bounds. The series 
i pursued is divergent; but if we examine the end of the sequence on 
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stopping at the rth term we find that, apart from the factor y/ tt (—) r_1 , 
we finish up with 


e ~ xJ T(r — i) 

■27 7 2( '- 1) ’ 



e~ x 'dx 



The integral repays examination. The factor e ~ x * steadily decreases 
throughout the range, showing that the integral must be less than 


r *> 

r(r -f i)e~ xt x~ 2r dx 

J X 


e ~ x> r(r — £) 

2x x 2(r - 1) * 


The integral remainder is thus seen to be less than our rth term. The 
comforting fact emerges that, if we decide to approximate by stopping 
at the rth term, the error is smaller in magnitude than the last term 
retained. Putting it into figures, if x is as low as 3 and we decide to 
retain only two terms, we can bank on the error being below <r 9 /108, 
which is of the order 10 -6 . Comparing this with the alternative pro¬ 
cedure of expanding e -x * and integrating termwise, the latter is evi¬ 
dently out of the running. It only remains to add that, if we calculate 
Erf (r) rather than Erfc ( x) we have 

f e~ xt (lx = f e~ z *dx — f e~ x *dx = \\/tt — f e~ x 'dx. 

J o J o Jx " J x 


10-5. Solution by definite integrals. 


Back in Chapter II we mentioned the matter of solving differential 
equations by means of definite integrals. We shall now give an ele¬ 
mentary exposition of this rather highbrow method by applying it to 
a convenient form of Bessel’s equation and eventually linking up with 
the asymptotic expansions. We propose to examine our chances of 
solving 


A 

dx 2 



2 ^±liL + y = 0, 

x dx 


which is the equation for x~ n C n (x), by means of the integral 

y =J Te ix, dt. 


Here T is a function of t alone, and the integral is eventually to be 
made definite, thus becoming a function of x alone. 

Whenever we attempt to solve an equation by using an integral, 
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much naturally depends on choosing a suitable integrand; here of 
course it is experience that counts. The above integrand may be con¬ 
sidered reasonable on the grounds that the cylinder functions behave 
very much like trigonometrical functions for large values of x. 

This substitution poses two problems, both of which are answered 
in the working. Firstly, what is to be the precise form of T; secondly, 
what are to be the limits of integration. In the meantime we have by 
differentiation and partial integration 



^ = fiTte ixt dl, 
ax J 


tfhi 

</x= 


—J Tt 2 e ixl dt 


t e izC r .ixt 

ti 2 ~ + f (r-r + 2tT) dt. 

txjj IX 


Putting these in the equation, we have with a little rearrangement 

[T(l - - J {r(L - 0) + (2 n - 1 )Tt}c ix ‘dt = 0. 

s 

The first term, in the square bracket, is presumed to be taken at the 
two limits of integration, which are not yet known. The obvious way 
of satisfying the identity is to make both the square bracket and the 
curly bracket vanish separately. To accomplish the latter means solving 
a small differential equation and we get 


T_ 

T 


<^i)< r=(1 



This fixes the form of T and one of the two problems is answered. The 
square bracket obviously vanishes for t = +1; and provided x is a 
positive real, it tends to zero also as t approaches 1 + i°°. This puts 
us m possession of the two solution^ 

Vl = L 1 eiX ' {l ~ t2 ) n ~ idt > y 2 =^ 1+ ’V'(l - 

It remains to interpret these. 
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Concerning the former, it fails to converge at either limit unless 
n > — £. Otherwise, if we substitute t = sin# we have 

r\* 

2 / cos(xsin#) cos 2n 6<W = y, = x~ n C n (x). 

Jo 


Apart from an irrelevant numerical factor, we have evidently reached 
J ri (x) as defined by Bessel’s second integral. 


10-6. The asymptotic solution. 

The second solution is more interesting, if for no reason other than 
its being complex and so providing two solutions by its real and ima¬ 
ginary parts. For positive real x we put 



so that the limits for v are 0 and oo. We then have 

/** .. u’/o i ivW”-* idv 

* = J. exp(.x(^2+ -jj 

or, 

= e "(- + £))"'*■ 


Since 


— i = cis( — \tt) = exp( — \i-rr), 


we can write the multiplier of the integral as 

expi{a: — (n + 

As for the integrand, we can expand by the binomial and integrate 
termwise. The interesting fact presents itself, that if n is half an odd 
integer (a phrase that by this time has a familiar ring about it), the 
binomial expansion terminates and the result is achieved in finite 
terms. Testing this out for n — h, we have 

—2y, = —2x“*Cj(x) = - exp i(x — \tt) f e~ v dv 

'X 

= - (sin* — t cosx). 

X 

Apart from irrelevant numerical factors, these evidently give .7* or 
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If n is not the half of an odd integer, we have the real part of the 
integral as 

/ [i _ [ n v * -|- . . . \dx 

•'» \ 2\(2x) 2 ) 

= r ( n + ^)( 1 -^ri )(,,2 -s ) +...) 

* \ 2!(2.r)* ^ / 

= PI> + i). 

Similarly for the imaginary part we have 


“ 2) „ _ (« ~ *)(« — ?)(« — I) 


V — 


+ 


f e~ v v 

\ 2x 3!(2x) 3 

= r ( « + i>f*zi - - Do* - »(!,»- V) 

“ l 2x 3!(2z) 3 ^ 


= Qr(n + i). 


dv 


The P and Q series are asymptotic. If pursued they ultimately diverge; 
but since they were derived from the binomial theorem, they can always 
be cut short with a remainder. If this was anything more than an 

. • 0 » we should be under the moral obligation of 

investigating that remainder; but we will waive the point. 

If we multiply -2y, by &r)"/vVr(n + i) and pick out the real 
and imaginary parts, we have the asymptotic formulae 


(1) 

- & 

(2) 

- Q‘< 

where 

(3) 

P = x — 


It is readily verified that these give the 
n is half an odd integer. 


usual formulae exactly when 


10*7. The modified functions. 

Seeing that Bessel’s equation was previously modified by writing 
“ f ° r *> we have *~ nI » and x~»K n as solutions of X Dg 

&y , 2n -f-1 dy 



1S4 


APPLIED BESSEL FUNCTIONS 


Hence also the integral 



is a solution, provided that 

[r(l - t 2 )e~*<] + f 'e~ xt {T'(t 2 - 1) - (2n - 1 )Tt}dt = 0. 


The square bracket is zero when t = 1, —1 and tends to zero as t 
tends to infinity. We therefore examine the integral 

y=[ e~ xl (l 2 - 1 )"-*&, t = 1 + 

* 1 X 



= e-.g)” + Wi ) { 


, , 4?i 2 — l 2 , (4n 2 — l 2 )(4n 2 — 3 2 ) , 
l + _ H- — -- ~r • 


l!(8x) 


2!(8x) : 



We know that yx n must be a linear combination of /„ and K n ; but 
as /„ tends to infinity with x we conclude that it is not present here. 
After re-arranging some factors we have the asymptotic formula 




4n 2 - l 2 
l!(8x) 



This can be checked for the value n = h 
The corresponding formula for /„ is 




4n 2 - l 2 
l!(8x) 



The asymptotic formula for the ber and bei functions of general order 
are rather fearsome affairs. Even those for zero order are not very 
alluring. We give them here without commenting further; they can 
be deduced from the asymptotic formula for 7 0 (x). 


berx = 


e* 


\Z(2ttx) 


cos a, beix = 


e* 


V(2 ni) 


sin a, 


X 7T 1 

\/2 8 _ 8y/2x 
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Adverting to tbe skin effect mentioned in 10-3 we have by defi¬ 
nition 

— her A x -f- i bci/.x. 

Differentiating with respect to x, we have in absolute values 

| I^kxi*) j = j bei'/iX — t ber'Avr | 


since the absolute value of i* is unity. Tlie real part of ) ! is 
on using the asymptotic value of we have 


P 

Q 





For any fixed values of r and c the exponent is negative, so that with 
increasing A- the absolute value | P/Q | tends to zero, even when r and c 
are quite close. 


(C150) 
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BIBLIOGRAPHICAL NOTE 


A few words on bibliography and the problems may be of interest. 
For the most part the problems are obviously the common property 
of the mathematical world. A few may be original, but I stake no 
claim to that. For one thing, the claim is scarcely worth making; 
for another, I should risk being non-suited. They came when bidden, 
possibly from the recesses of memory, in which case they left their 
ancestry behind. 

Take, for example, the transverse vibrations of a non-uniform 
rod. Anybody conversant with the literature of the subject knows 
that it was originally tackled by Kirchhoff. It was done again by 
J. \V. Nicholson and it appeared in Proc. Roy. Soc. xciii A (1917). 
Five years later the same journal published another account by 
D. Wrinch in Vol. 101 (1922). 

Similarly Bernoulli's chain problem was modified by Greenhill, as 
stated in the text. The matter was investigated experimentally by 
J. R. Airev who gave his results in Phil. May. (6) xxi, 1911, p. 73G. 
He adds without proof that the uniform chain with end load is soluble 
in terms of./ and Y of orders 0, 1, 2; and he states correctly that the 
problem figures in Routh, Advanced Rigid Dynamics. It also figures 
as No. 32 in the miscellaneous examples at the end of Gray’s treatise 
mentioned below. 

Tiie lengthening pendulum, which I regard as an outstandingly 
good illustration, figures in the last-mentioned source as No. 49 and 
the authors quote Lecornu in the Comptes Rendus of 15th January, 
1894. This I am pretty sure I first met in a record of a conversation 
between Einstein and Max Planck, but I have failed to trace it. The 
upshot is that if anyone feels aggrieved at having their material pressed 
into my service, I hasten to assure them that it was done inadvertently 
and I am prepared to make reasonable amends. And now for the 
books, to all of which I am indebted. 

1. Watson: Theory of Bessel Functions. 

This impeccable and exhaustive treatise is the standard text for 
mathematicians in English. It contains ninety pages of tables, a long 
bibliography, and not a graph in all its eight hundred odd pages. 

2. Gray, Mathews and MacRobert: A Treatise on Bessel Functions 
and their Applications to Physics. 
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This well-known work contains fifty pages of tables, a bibliography, 
a graph and a large collection of miscellaneous examples. The subject 
is treated by contour integrals with more rigour than herein; the 
applications take up half the book and are more advanced than mine. 

3. Riemann: Part idle Different ialgleichungen. 

One of the world’s most delightful mathematical books; some¬ 
times known as Riemann-Weber. To those familiar with the older 
edition, the latest issue looks like vandalism. 

4. E. B. Wilson: Advanced Calculus. 

This miniature encyclopaedia is invaluable to those wishing to 
improve on the calculus of their student years. I am indebted to the 
author for stabilizing my belief that the recurrence formulae made a 
suitable approach and that solution by integrals, leading to the asymp¬ 
totic formulae, was feasible by elementary means. 

5. Jnhnke und Emde: Funktionentafeln mil Fonneln und Kurvm. 

The reputation of this book has grown steadily and deservedly over 

the last thirty years. It is now procurable printed in English and 
German on opposite pages. It contains an unparalleled amount of 
information on tabulated functions, together with excellent graphs. 
It should be in the possession of all who meditate applications of 
higher mathematics. 

6 . Prescott: Applied Elasticity. 

This gives the necessary analysis and the discussion of such of our 
problems as depend on the theory of elasticity. Its form is more easily 
assimilable than the standard work by Love. It also treats the in¬ 
stability of the deep girder. 

7. Case: Strength of Materials. 

This is useful for the ordinary theory of beams and struts, and it 
gives the theory of the instability of the deep cantilever. 

8 . Carslaw: Conduction of Heat. 

The original work was later issued in two volumes, the first dealing 
with Fourier series and the second with heat conduction. The latter 
gives the applications of Bessel functions. 


9. Lamb: Hydrodynamics. 

Tins standard work gives several applications of Bessel functions 
to fluid motion. 


10 . Rayleigh: Theory of Sound. 

leJathere ^ Vibrati °“ P roblcnis will be found treated at greater 
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-zero order, 48, 74. 

— integrals, 154, 157, 101, 182. 
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Clare, John, 88. 

Common zero, 24. 

Comparison thuorom, 25. 

Complomontary function, 15. 

Complex hyperbolic function, 11G. 

— zero, 62. 

Continuity, equation of, 132 . 

Contour integration, 10, 10. 

Critical load, strut, 08. 

Cycloid, 154. 


Definite integrals, solution by, 1.80. 
Difference equation, gamma functions, 7 
Diffraction of light, 71. 

Dini, U., 170. 

Disc, non-uniform, 149. 

— vibration of, 145. 

Duplication formula, 12, 172. 

Edwards calculus, 10. 

Elliptic arch, 102. 

— functions, 13. 

Error function, 3, 179. 

Estuary, 120. 

Filler, L., 8S, 91. 

— constant, 70. 

Eulcrian integral, first kind, 9. 

-second kind, 6. 

Existence theorems, 19. 

Faraday, M., 170. 

Fourier-Bessel scries, 109. 

Frosners integrals, 71. 

Frobcnius method, 100, 175. 

Fundamental solutions, 20. 

Gauss, 0. F., 93. 

— function, 9. 

— table, II. 

Generalized factorial, 7. 

Greenhill, G., 92. 

Half swing, 59, 84. 

Hankcl'functions, 174. 

Heat conduction, 5, 92. 135, 1G3. 
Heaviside-Bessel line, 117. 

Heine, H. E., 32. 

Her, hei functions, 175. 

Hormito functions, 30. 

Heterogeneous medium, 93. 

Humidity, 04. 
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Imaginary argument, 100. 

I mpropter integral, 6. 

Inconsistent equations, 18. 

Independent solutions, 16, 20. 

Indicial equation, 45, 73. 

Infinite integral, 4. 

Integral orders, 33. 

Integrating factor, 17. 

Interlacing of zeros, 25, 32, 101. 
Irrotational motion, 132. 

Jacobi series, 153. 

Kelvin quoted, 64. 

— functions, 174. 

Kennedy tables, 116. 

Kcr, kei functions, 175. 

Laguerre's equation, 28, 29, 30. 

Laplace, P. S., 32, 131. 

Laurent series, 152. 

Leaky telegraph line, 116. 

Least squares, 92. 

Legendre, A. M., 8. 

— equation, 21, 22, 27, 28. 

Lipschitz integral. 165, 168. 

Lissajous figure, 159. 

'Locomotive, 70. 

Lornmel integrals, 51, 72, 104, 110, 163. 

Macdonald, H. M., 107. 

Mallctt, telegraphy, 116. 

Mean value theorem, 25, 29. 

Modified functions, 100. 

-first kind, 101. 

-second kind, 107. 

Monotonic, 35. 

Napier, 78. 

Neumann, K. G., 75, 78. 

Newton s method, 61. 

Nielsen, N., 174. 

Nodes, 62.— - - ■*——----T~ 

Normal erc^L.6; { \ J VI /rtj A 

— form,'22. 

-Bessel’s equation, 34. 

— mode of vibration, 90. 

Order of-function,.27*. . 

OrthogonaLfunetionsr^S.—- 

— property, 50, 104. 

Paper clip, 65. 

Parameter, equations with, 27. 


Particular integral, 15. 

Pendulum, 56, 80. 

Physical constants, 92. 

Poisson's integrals, 157. 

— ratio, 144. 

Porosity, 94. 

Probability, 4. 

— integral, 6. 

Propagation constant, 116. 

Quantum theory, 30. 

Quarter periods, 58. 

Rayleigh's principle, 68. 

— sound, 145, 169. 

Recurrence formula?, 27, 31, 102, 107. 
Reduced equation, 14. 

Reduction formula?, 10, 13, 79. 
Remainder, 4. 

Repeated roots, 73. 

Ribbon, steel, 68, 98. 

-blunt nosed, 70. 

-sharp nosed, 70. 

Rolle's theorem, 7, 24, 33. 

Schlomilch expansion, 171. 
Schrodinger wave equation, 30. 
Second solution, 73, 106. 

Series impedance, 117. 

— solution, 45. 

Shakespeare, 88. 

Shunt admittance, 117. 

Singularities, 45. 

Skin effect, 178, 185. 

Sonine's integral, 167. 

Stability of cantilever, 65. 

— of triangular strut, 68. 

— of wire, 62. 

Standard form, 40, 103. 

Starr, A. T., 119. 

Steady state, 136. 

JStrine, vibrating, 59. 

.Strut! tapered, 95. 

— variable, 68. 

Tambourine, 88. 

Viffwfcd bar. 85. 139. 

•JL'iVlor's theorem, 22. 

-Temple, G., 68. 

Thermocouples, 136. 

Tidal motion, 124. 

Tie bar, 119. 

Tortuous curve, 65. 
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Transverse vibrations, 139. 

TschebychefT function, 30. 

Variable mass, 82. 
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Vertical wire, stability, 63. 

Vibration, disc, 145. 
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— string, 59. 

Viscous liquid, 174. 

Vortex motion, 132. 


Wacnlcn, gruppentheorie, 30. 
Watson, 0. X.. 13. 

Wave equation, 30. 

Waves, 131. 

Weathering, 94. 

Weber, H., 76, 78, 79. 

— discontinuous integral, 163. 

— function, 29. 

Wet string. 60. 

Weyl, quantum mechanics, 30. 
Whittaker, E. T., 13. 

Wire, stability of, 63. 

Zero order, 101. 

Zeros, common, 24. 

— infinity of, 36. 

— interlacing, 25, 32. 

— of the modified equation, 110. 
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